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The  analysis  of  the  sleep  record  is  a required  process 
in  sleep  research.  Usually,  several  EEG  (electroencephalo- 
gram) channels,  and  one  or  two  EOG  (electrooculogram)  chan- 
nels are  recorded  simultaneously  on  the  paper  during  the 
sleep  period.  The  record  of  a single  night's  sleep  data  is 
about  500  meters  long.  In  order  to  provide  a fast  and  an 
objective  method,  many  computerized  sleep  EEG  systems  have 
been  developed.  Though  several  systems  showed  good  perfor- 
mance, there  still  exist  problems  of  a gestalt  perception  of 
the  waveform  and  a data  analyzing  method,  which  is  needed  to 
utilize  the  various  relationships  of  the  data.  This  disser- 
tation describes  a method  of  sleep  EEG  data  analysis.  An  ex- 
pert system  is  designed  to  provide  the  method  which  can 
simulate  the  signal  understanding  procedures  of  the  human 
expert. 
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In  this  system,  the  EEG  signal  of  the  specified 
waveforms  is  regarded  as  a language,  which  has  a set  of 
words  and  a simple  grammar.  Each  waveform  is  related  to  a 
word,  and  channel  relationships  or  waveform  relationships  in 
time  constitute  the  grammar.  Using  these  words  and  grammar, 
the  context  analysis  is  done  on  the  signal  data.  Then  the 
system  synthesizes  all  channel  activities  and  makes  a syn- 
thetic report  using  expert  knowledge  base.  A mixed  tech- 
nique of  a rule-based  system  and  a frame  scheme  is  utilized 
for  the  system  reasoning.  The  frontend  processor  for 
language  token  generation  is  implemented  using  an  A/D  con- 
verter and  a TM-990/101M  microprocomputer  having  32K  bytes 
memory.  The  main  processor  for  the  language  token  processing 
is  implemented  on  a PDP-11/23+  microcomputer  with  the  UNIX 
operating  system.  The  rapid  eye  movement  (REM)  period 
detection  and  the  sleep  staging  which  utilize  multichannel 
information  and  the  expert  knowledge  base  are  shown  as  com- 
ponents of  the  expert  system. 
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CHAPTER  1 
INTRODUCTION 

Analysis  of  the  sleep  EEG  record  is  a productive  method 
in  the  study  of  sleep.  Usually,  several  EEG  (electroen- 
cephalogram) channels,  and  one  or  two  EOG  (electrooculogram) 
channels  are  recorded  simultaneously  on  the  paper  during  a 
sleep  period.  The  record  of  a single  night’s  sleep  data  is 
about  500  meters  long.  In  order  to  provide  a fast  and  an 
objective  method  for  sleep  EEG  analysis,  computerized  sleep 
EEG  systems  [Ma72]  [Fr73]  [Ga73]  [Sm78]  [Ha83a]  have  been 
developed  over  a period  of  years.  Although  several  systems 
performed  well,  today  there  still  exist  significant  problems 
of  visual  perception  and  contextual  understanding  of  the 
waveforms  in  multichannel  sleep  EEG  analysis  using  a comput- 
er. A data  analyzing  method  is  needed,  which  mimics  the  hu- 
man perception  procedures  and  which  includes  more  of  the 
corresponding  expert's  knowledge  base. 

A perceptual  understanding  and  a qualitative  analysis 
of  recorded  signal  data  need  various  kinds  of  knowledge 
bases.  Most  of  the  data  analysis  needs  information  transfor- 
mation which  uses  different  knowledge  bases  depending  on  the 
required  analysis.  In  particular,  when  the  data  analysis  in- 
cludes a visual  perception  process  and  a contextual  under- 
standing simultaneously,  an  expert  knowledge  base  is  neces- 
sarily required  to  understand  and  analyze  the  data. 
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In  other  words,  in  a human  perception  process,  a pic- 
ture can  be  regarded  as  an  integrated  knowledge  source  that 
can  transfer  any  kind  of  information  depending  on  the  situa- 
tion, location  and  the  people  who  are  looking  at  it.  The 
image  of  the  picture  may  stimulate  a memory  about  a similar 
scene  or  related  people,  or  cause  us  to  think  about  implied 
ideas  which  the  picture  might  suggest.  However,  recognition 
processes  of  the  picture,  i.e.,  information  transformation 
from  picture  to  us,  are  performed  using  a knowledge  base 
which  has  been  built  through  our  daily  life.  Depending  on 
the  knowledge  base,  a different  analysis  result  can  be  ob- 
tained with  the  same  picture.  The  necessary  information  may 
be  obtained  directly  from  the  picture  and  indirectly  from 
the  embedded  knowledge  structure  of  the  picture  using  our 
knowledge  bases.  The  understanding  procedures  of  the  sleep 
EEG  record  are  the  same  as  the  picture  understanding  pro- 
cedures except  for  the  fact  that  a sleep  data  record  and  a 
heuristic  knowledge  base  of  human  EEGer  ( electroencephalo- 
grapher)  are  used. 

Let's  take  a close  look  at  the  sleep  data  and  data  pro- 
cessing mechanism  of  a human  expert  to  see  why  most  of  the 
conventional  methods  [Ha83a]  [Ge80]  [La80]  [Sm78]  do  not 
work  properly  and  how  the  data  processing  mechanism  can  be 
realized  appropriately  using  a microcomputer.  Fig.  1*1> 
Fig.  1.2,  Fig.  1.3  and  Fig.  1.4  show  examples  of  sleep  stage 
1 REM  (rapid  eye  movement),  eye  movement  in  sleep  stage  0, 
K-complex  in  sleep  stage  2,  and  delta  activities  in  sleep 
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Fig.  1.1  Sleep  Stage 
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Fig.  1.2  Sleep  Stage 
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Fig.  1.3  Sleep  Stage 


6 


N 

Q_ 

C\J 

N 

Ll_ 

O 

i— ♦ 

CO 

CO 

Lu 

CJ 

o 
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stage  4,  respectively.  Although  slow  waves  which  look  like 
a rapid  eye  movement  appeared  on  each  EOG  channel  (LE-A2), 
only  the  slow  wave  of  the  Fig.  1.1  is  regarded  as  a REM 
period  in  sleep  stage  1 . A big  slow  wave  on  EOG  channel  is 
not  regarded  as  a REM  in  Fig.  1 .2  because  of  the  dominant 
alpha  activity  in  the  occipital  channel  (03-OzPz)  in  the 
previous  seconds  In  Fig.  1 .3>  though  a big  slow  wave  on  EOG 
channel  appeared  without  any  dominant  activities  in  previous 
seconds,  the  simultaneous  slow  activities  followed  by  a fast 
activity  in  the  occipital  channel  makes  it  a K-complex  wave 
instead  of  REM.  In  Fig.  1.4  the  delta  activity  in  all  chan- 
nels shows  that  the  EOG  channel  slow  waves  are  not  REM. 
These  are  good  examples  of  how  the  waves  are  recognized  for 
REM  period  classification. 

The  data  analysis  mechanism  can  be  grossly  described  as 
a two-step  procedure.  The  first  step  is  the  visual  percep- 
tion of  the  waveform  such  as  alpha,  beta,  delta,  sleep  spin- 
dle, etc..  The  second  one  is  the  resolution  of  the  signifi- 
cance of  the  waveforms,  i.e.,  contextual  understanding  of 
the  waveforms  such  as  waveform  priorities  or  total  activity 
time  of  a specified  waveform.  Although  the  waveforms  are 
classified  depending  on  the  frequency  range,  it  is  not  clear 
enough  which  frequency  range  a certain  waveform  belongs  to. 
The  reasons  are  the  following: 

a)  The  term  "frequency"  itself  is  not  clearly  defined 
since  sleep  waveforms  are  usually  composed  of  many  dif- 
ferent frequencies.  The  EEGer  counts  the  number  of 
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peaks  or  the  number  of  zero  crossing  points  in  a fixed 
interval  instead  of  decomposing  the  waveforms  into 
sinusoidal  waves  for  frequency  measurement. 

b)  Depending  on  the  context  of  data,  only  a certain 
frequency  component  is  counted  in  visual  analysis. 

c)  Each  waveform  amplitude  is  subject  dependent,  and 
all  the  waveforms  in  one  subject  are  not  always  big  or 
small.  For  example,  a subject  with  big  alpha  wave  does 
not  always  have  big  beta,  delta,  or  sigma  wave.  There- 
fore, an  amplitude  threshold  value  can  be  changed  for  a 
certain  waveform  detection. 

It  is  difficult  even  for  a human  expert  to  say  which 
component  is  the  dominant  one  in  many  instances.  Also,  the 
boundary  of  a waveform  itself  often  is  not  clearly  defined. 
That  is,  most  of  the  waveforms  have  fuzzy  shapes  and  fuzzy 
boundaries.  These  are  the  principal  reasons  which  account 
for  the  disagreement  of  the  analysis  results  between  dif- 
ferent human  experts  and  why  many  sleep  laboratories  use 
different  sleep  staging  rules. 

However,  depending  on  the  number  of  recognized 
waveforms  or  the  duration  of  a wave  burst  or  wave  activity, 
the  contextual  understanding  of  each  waveform  is  performed. 
Usually,  the  context  understanding  mechanism  for  these  kinds 
of  fuzzy  signals  is  difficult  to  implement  using  convention- 
al programming  techniques.  When  a visual  perception  [ R i 8 3 ] 
[ Wi7  9 3 [To74]  is  required  for  data  analysis,  specific 
methods  or  algorithms  must  be  used  depending  on  the  problem. 
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These  kinds  of  methods  or  algorithms  cannot  be  generalized 
as  the  human  expert  knowledge  cannot  be  generalized.  Gen- 
erally, different  domains  of  expertise  have  different 
knowledge  bases. 

Then  how  can  we  make  a microcomputer  interpret  the 
sleep  data  as  does  a human  expert?  To  make  a computer  inter- 
pret a picture  (data)  means  not  only  giving  all  the  informa- 
tion about  the  picture  (data)  itself  but  also  giving  a meta 
knowledge  which  can  operate  on  a given  knowledge  base  to  ex- 
tract the  necessary  information  from  the  picture  (data). 
Most  of  the  work  in  the  field  of  Artificial  Intelligence  has 
centered  around  the  problem  of  understanding  in  cognitive 
science,  i.e.,  structured  knowledge  representation  for  in- 
formation storage  and  retrieval,  language  understanding  for 
communication  and  knowledge  base  management  for  learning 
i.e.  knowledge  accumulation.  Techniques  already  developed 
to  impart  intelligence  to  a computer  have  been  utilized  in 
many  fields  such  as  CAD/CAM  (computer  aided  design,  and 
manufacturing),  Signal  Processing,  Expert  Systems,  Data  Base 
Management , etc . . 

To  analyze  the  sleep  EEG  data  using  a microcomputer,  an 
expert  system  has  been  constructed  based  on  the  human  visual 
perception  and  contextual  understanding  mechanism  which 
utilize  the  waveform  detection  method  [Sm78]  and  symbolic 
signal  processing  techniques  of  Artificial  Intelligence. 
The  system  performs: 

a)  the  recognition  of  the  specified  waveforms  in  time 
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domain , 

b)  the  examination  of  the  data  back  and  forth  in  time 
and  across  all  channels,  if  necessary  (qualitative 
analysis) , 

c)  the  quantification  of  the  data,  i.e.,  summation  of 
the  total  activity  time  or  the  total  number  of  bursts 
in  each  epoch  (quantitative  analysis), 

d)  the  rule  selection,  which  is  dependent  on  the  re- 
quired analysis  (system  flexibility), 

e)  the  rule  management,  i.e.,  knowledge  integration 
(learning  mechanism), 

f)  an  explanation  of  the  reasoning  mechanism  (user 
friendliness) . 

At  present,  the  user  friendly  interface  (such  as  natur- 
al language  communication  or  data  display)  and  knowledge 
base  management  system  have  not  been  fully  considered.  Only 
the  basic  frame  of  the  sleep  EEG  analyzing  expert  system  has 
been  designed  to  demonstrate  how  it  can  simulate  a human  ex- 
pert. In  the  present  system,  a data  analysis  report  more 
detailed  than  the  one  provided  by  the  human  expert  is  gen- 
erated using  the  multichannel  EEG/EOG  signal  inputs.  The 
symbolic  signal  processing  is  performed  using  a simple  sig- 
nal language  understanding  mechanism  and  the  heuristic  rules 
of  the  human  expert. 

The  functional  structure  of  the  data  analyzing  system 
is  divided  into  two  blocks,  a token  generator  and  a token 
processor.  In  the  token  generator,  the  multichannel  EEG 
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signals  are  converted  to  a token  stream  which  contains  the 
information  about  all  channel  waveform  activities.  From  the 
token  stream,  a sequence  of  words  is  generated,  which 
describes  a channel's  waveform  activity.  This  sequence  of 
words  is  interpreted  in  the  token  processor.  The  token 
processor's  interpreting  mechanism  is  very  similar  to  that 
of  high  level  computer  language  compiler.  A set  of  words  is 
regarded  as  a signal  interpretation  language.  Each  waveform 
is  related  to  a word,  and  the  channel  relationships  or  the 
time  domain  waveform  relationships  constitute  the  grammar  of 
the  signal  interpretation  language.  Using  the  selected  words 
and  grammar,  the  context  understanding  of  the  signal  is  per- 
formed. Although  all  the  signal  characteristics  are  not 
represented  in  the  grammar,  basic  signal  patterns  and  rela- 
tionships are  recognized  and  analyzed,  using  a reasonably 
simple  grammar.  Once  a channel  description  list  is  ob- 
tained, an  epoch  description  is  obtained  using  a knowledge 
base  of  the  heuristic  rules  of  the  human  expert.  A hybrid 
technique  of  a rule  based  system  and  a frame  scheme  is  util- 
ized for  the  rule  implementation  and  the  system  reasoning. 

The  token  generator,  used  as  a frontend  processor,  is 
implemented  using  TI-9900  microprocessor  and  an  A/D  convert- 
er in  which  four  single  end  input  channels  are  utilized.  The 
token  processor,  used  as  a main  processor,  is  currently  im- 
plemented on  a GOULD  CONCEPT  32/8280  minicomputer  using  the 
UNIX  operating  system.  (In  the  near  future  it  will  be  im- 
plemented on  PDP-1 1/23+) . The  TI  assembly  language,  and  the 
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high  level  language  "C",  and  "LISP"  are  used  in  the  token 
generator  and  processor.  The  standard  software  of  the  UNIX 
operating  system,  the  LEX  (lexical  analyzer)  [Le75l  and  the 
YACC  (yet  another  compiler  for  compiler)  [Jo75]  program,  are 
utilized  for  generating  the  signal  language  compiler. 

In  Chapter  2,  most  of  EEG  analysis  methodologies  which 
are  related  to  sleep  EEG  analysis  are  reviewed.  In  Chapter  3 
a more  elaborate  approach  to  a sleep  EEG  analysis  expert 
system  is  discussed.  In  Chapters  4 and  5,  the  design  con- 
cept, the  detailed  structure  of  the  system  and  the  system 
status  are  described.  The  system  performance  and  the  con- 
clusions are  discussed  in  Chapter  6. 


CHAPTER  2 

OVERVIEW  OF  EEG  SIGNAL  ANALYSIS 
In  this  chapter  the  EEG  signal  source,  model  and  com- 
ponents are  briefly  reviewed  and  the  basic  EEG  analysis 
methods  are  discussed  to  present  a better  picture  of  the  ad- 
vantages and  limitations  of  each  method  for  the  sleep  EEG 
analysis  application. 

EES  Signal  Source , Model  and  Components 
It  is  generally  known  that  the  scalp  EEG  derives  from 
graded  synaptic  potentials  generated  by  pyramidal  cells  in 
the  cerebral  cortex,  which  are  triggered  by  rhythmic 
discharges  from  the  thalamic  nuclei.  A pacemaker  system  si- 
tuated in  the  thalamus  and  in  the  reticular  formation  prob- 
ably regulates  the  synchrony  of  the  cortical  signals  [An68] 
[Ge75].  The  relatively  slow  time  course  of  EPSPs  (excitato- 
ry post-synaptic  potentials)  and  IPSPs  (inhibitory  post- 
synaptic  potentials)  is  comparable  with  the  EEG  and  the  sum- 
mation is  facilitated  by  the  columnar  structure  of  the  neu- 
rons reaching  from  upper  to  lower  layers  of  the  cortex 
[C08O]  . 

Several  models  have  been  devised  to  locate  the  signal 
source,  to  characterize  the  wave  propagation,  and  to  analyze 
the  signal  characteristics.  One  of  them  is  the  current  di- 
pole layers  model  [Nu8la]  which  assumes  that  the  scalp  po- 
tentials are  due  to  current  dipole  layers  occupying  various 
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surface  areas  of  the  cortex.  Approximately  2/3  to  3/4  of 
cortical  neurons  are  oriented  perpendicular  to  the  cortical 
surface.  The  large  number  of  interconnections  makes  synchro- 
nous activity  possible.  Using  this  current  dipole  model,  a 
relationship  between  cortical  current  density  and  resulting 
scalp  potential  was  obtained  by  means  of  arguments  about 
volume  conduction  through  brain,  skull,  and  scalp.  And  a 
theoretical  model  of  neural  interaction  at  the  macroscopic 
level  of  brain  hierarchy  was  developed  to  predict  cortical 
current  density  and  the  resulting  spatial- temporal  proper- 
ties of  scalp  EEG . It  was  noted  that  the  EEG  frequency  and 
amplitude  is  dependent  on  the  threshold  value  which  is  the 
average  difference  between  the  number  of  excitatory  and  in- 
hibitory synapses  on  a neuron,  weighted  by  the  appropriate 
threshold  [Nu8la]. 

A neuronal  population  model  [Ch73],  which  is  based  on 
the  hypothesis  that  the  surface  potentials  are  a combination 
of  ESPSs  and  IPSPs  which  occur  both  at  different  depth  and 
different  latency,  was  used  to  quantify  the  electrogenesis 
of  the  evoked  surface  potentials  and  to  examine  the  rela- 
tionship between  scalp  potentials  evoked  by  visual  stimula- 
tion and  neuronal  functioning.  This  model  is  used  to  deter- 
mine the  location  (depth)  of  the  signal  sources  and  the  ap- 
proximate size  of  the  neuronal  population  responsible  for 
the  potential.  Using  an  array  of  scalp  electrodes  over  the 
occipital  cortex  and  computer  simulated  experimental  data, 
it  was  shown  that  the  sequence  of  neuronal  excitation,  the 
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magnitude  of  excitatory  current,  the  depth  of  the  neuronal 
population,  and  the  period  of  activity  of  the  EPSP  and  IPSP 
could  be  determined  [Ha75]  [Ch77]. 

For  the  quantitative  analysis  of  the  EEG  signal,  a 
hierarchical  model  [Sa80]  which  utilizes  a recursively  es- 
timated autoregressive  (AR)  model,  a piecewise  stationary 
model,  composite  source  model  and  character  string  and  syn- 
tactic model,  was  developed  considering  that  EEG  signals  may 
appear  to  be  stationary  over  short  periods  of  observation, 
but  exhibit  markedly  nonstationary  behavior  when  observed 
for  a long  period.  In  this  model  a discrete  time  signal 
model  of  a piecewise  stationary  stochastic  sequence  is  used 
for  data  segmentation.  Each  AR  model  for  a data  segment 
represents  a type  of  waveform.  The  transition  of  the  ac- 
tivity is  detected  using  a one-step  prediction  error.  Using 
this  segmented  or  piecewise  stationary  model,  the  decomposi- 
tion of  the  original  EEG  signal  is  done  to  produce  a compo- 
site EEG  signal  source.  This  composite  source  generates  the 
output  of  character  strings  by  random  switching  among  the 
subsources.  The  output  of  the  composite  source  is  then 
analyzed  using  a stochastic  grammar  in  which  the  detected 
substrings  are  regarded  as  a nonterminals.  To  remove  the 
subject  dependency  of  the  model,  a global  subsource  set  for 
the  given  population  was  developed.  With  a ninety  minutes 
data  record,  about  800  segments  and  46  identified  source 
models  were  obtained  at  a speed  of  five  to  ten  times  real 
time  processing  of  a single  EEG  channel. 
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Although  certain  models  for  EEG  signal  sources  could 
locate  the  spatial  origin  and  could  explain  wave  propagation 
characteristics  and  wave  dispersion  relationships  in  a 
nonhomogeneous  volume  conduction  such  as  the  brain,  cerebral 
cortex,  skull  and  scalp,  there  is  not  sufficient  knowledge 
of  the  EEG  signal  origin  and  the  significance  of  EEG  activi- 
ties such  as  the  relationships  between  physiological  and 
clinical  states.  However,  the  analysis  of  the  EEG  record  is 
a useful  aid  in  the  research  of  the  clinical  diagnosis,  the 
evoked  response  potential  and  the  sleep  analysis  etc.  [Ge75] 
[Re77]  [Ba79]  [Ge80]  [Dz84]  [B08O]  [Az82], 

Two  kind  of  signals  are  of  main  concern  in  brain  wave 
research.  The  first  signal  is  the  event  related  potentials 
(ERP).  External  inputs  to  the  human  brain  such  as  sensory 
inputs  generate  these  potentials.  For  instance,  the  N100 
component  [C08O],  which  has  about  100  ms  latency,  can  easily 
be  evoked  by  auditory,  visual  and  tactile  stimuli.  It  can 
be  seen  in  a central  channel  EEG  if  a simple  stimulus,  such 
as  a hand  clap,  is  given.  The  amplitude  of  this  component 
is  dependent  on  a variety  factors  including  intensity  of  the 
stimulus,  its  expectedness  and  attention  paid  to  it.  Usual- 
ly, event  related  potentials  are  masked  by  the  background 
noise  which  normally  exists  in  the  low  frequency  band.  To 
extract  the  ERP  from  the  ongoing  activity  which  seems  like 
background  EEG  or  noise,  special  averaging  techniques  such 
as  cross-correlation,  median  representative,  synchronous  and 
latency-corrected  averaging  have  been  used  [We72]  [Mc77l 
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[ Vi77] [Au8l 3 . The  second  signal  is  the  clinical  EEG  signal. 
The  clinical  EEG  signal  can  be  distinguished  by  physically 
looking  at  the  record  because  the  signal  amplitude  is  gen- 
erally much  larger  than  the  background  activities.  Physical 
interpretation  and  clinical  interpretation  are  possible  by 
visual  analysis  of  the  clinical  EEG.  For  physical  interpre- 
tation it  is  necessary  to  recognize  and  classify  specified 
EEG  waveforms  which  have  an  individual  identity  and  temporal 
patterns  which  result  when  they  occur  in  a sequence.  For 
clinical  interpretation  the  information  is  obtained  from  the 
comparison  of  empirical  observations  between  the  EEG  and 
clinical  observations.  The  sleep  EEG  is  one  of  the  clinical 
EEG  signals.  The  computerized  physical  interpretation  of 
the  clinical  EEG  data  is  the  main  concern  in  this  study. 

The  signal  used  in  sleep  research  is  in  the  • frequency 
range  of  0.5  Hz  to  120  Hz.  The  specific  frequency  bands 
used  in  the  present  system  include  : rslow  (0. 3-6.0),  oslow, 
pslow,  qslow  (0. 3-3.0),  delta  (0. 5-2.0),  theta  (3.0-7. 0), 
alpha  (8.0-12.0),  sigma  spindle  (11.75-16.0),  beta  (16.0— 
32.0),  and  artifact  (32.0-120).  Although  waveforms  are 
classified  in  terms  of  frequencies,  the  number  of  local 
peaks  in  a fixed  time  interval  or  the  half  wave  duration  is 
utilized  instead  of  sinusoidal  wave  number  for  frequency 
measurement  in  visual  sleep  EEG  analysis.  There  are  slight 
differences  in  the  range  of  each  frequency  band  depending  on 
the  laboratory  standards. 
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A variety  of  methods  have  been  tried  to  analyze  EEG 
data  [Ge75]  [Re77]  CBa79l.  They  can  be  roughly  divided  into 
two  groups.  One  is  the  frequency  domain  analysis  in  which 
the  spectral  parameters  obtained  using  various  spectral  es- 
timation methods  are  utilized  for  classical  pattern  recogni- 
tion methods.  The  other  is  time  domain  analysis  in  which 
the  waveform  descriptors  such  as  wave  period,  peak  ampli- 
tude, derivatives,  statistical  characteristics,  and  wave 
criteria  are  utilized  to  describe  each  waveform  activity 
directly  or  indirectly.  In  the  following  section,  various 
spectral  estimation  methods,  frequency  domain  methods,  time 
domain  methods,  and  pattern  recognition  methods  for  clinical 
EEG  classification  are  reviewed. 

£££  Sppctrai  Eptj.rgjiipp 

An  appropriate  spectral  estimation  is  the  key  to  the 
spectral  analysis.  In  general,  for  a consistent  spectral 
estimation  it  is  assumed  that  the  signal  is  stationary  or  at 
least  piecewise  stationary  [Op75].  Though  piecewise  sta- 
tionarity  of  the  EEG  [Co77]  [Ka73J  and  Gaussian  distribution 
behavior  [ El 67 ] [Mc75]  were  shown  in  some  cases,  it  is  also 
shown  that  the  signal  characteristics  such  as  stationarity 
and  Gaussianity  are  changing  [Ge75]  depending  upon  the 
behavioral  state  of  the  subject.  Therefore,  the  change  of 
the  state,  i.e.,  the  duration  of  the  same  state  and  the 
boundary  of  the  state  change  need  to  be  considered  for  a 
proper  spectral  estimation. 
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The  power  spectrum  can  be  estimated  using  the  indirect 
method  [B158],  direct  method  [Du73]  [Ma73]  or  data  model 
[ Fe7 1 ] [Ze69].  The  correlogram  and  periodogram  have  been 
used  for  indirect  and  direct  method,  respectively.  AR  (au- 
toregressive), MA  (moving  average)  and  ARMA  (autoregressive 
and  moving  average)  techniques  have  been  used  to  generate 
data  models.  The  direct  method  has  been  widely  used  with 
the  FFT  (fast  Fourier  transform)  or  FWT  (fast  Walsh 
transform) . Data  model  has  been  used  as  a minimum  variance 
and  unbiased  estimator  in  modern  spectral  analysis  to 
represent  the  time  domain  data  structure.  The  direct  method 
and  data  models  related  mainly  with  the  clinical  EEG 
analysis  are  discussed  in  the  following  section. 

Pj.rec£  Metlap^ 

The  power  spectrum  is  usually  obtained  by  taking  the 
Fourier  transform  of  the  properly  windowed  data  and  squaring 
its  absolute  value  [De73l  [Ma733 . The  characteristics  of 
the  spectrum  obtained  using  the  direct  method  are  equivalent 
with  the  ones  obtained  from  the  periodogram.  Although  one 
can  easily  select  an  unbiased  estimate  of  autocorrelation 
function  for  the  calculation  of  the  periodogram,  the  finite 
limits  of  the  data  cause  a biased  estimate  of  the  power 
spectrum,  and  the  variance  of  the  estimated  power  spectrum 
does  not  approach  zero  even  though  the  number  of  the  data 
points  used  for  the  periodogram  estimation  approaches  infin- 
ity [Op75],  Therefore,  a consistent  estimate  of  the  power 
spectrum  can  be  obtained  from  the  modified  periodogram.  To 
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reduce  the  variance  of  estimates  the  average  of  the  indepen- 
dent periodograms  has  been  used.  To  increase  the  spectral 
resolution,  an  appropriate  spectrum  window  is  applied  to  the 
periodogram  [Ba53J.  But  the  averaging  process  increases  the 
bias  of  the  estimates  and  decreases  the  spectral  resolution 
as  the  segmentation  number  increases  in  a fixed  data  record 
length.  Another  approach  to  this  problem  is  to  smooth  the 
spectrum  by  applying  an  appropriate  window  function  directly 
to  the  data  sequence  before  computing  the  periodograms 
[We70].  Overlapping  segmentation  has  also  been  utilized  to 
minimize  the  decrease  of  the  resolution  caused  by  segmenta- 
tion of  the  fixed  data  in  the  averaging  process,  where  the 
modified  periodograms  are  not  independent.  That  is,  one  can 
only  improve  the  variance  of  estimates  at  the  expense  of  the 
spectrum  resolution.  Therefore,  when  a certain  spectral  es- 
timation procedure  is  applied  to  the  unknown  signal,  the  in- 
terpretation of  the  obtained  spectrum  has  to  be  done  consid- 
ering the  spectral  estimation  procedures,  i.e.,  windowing 
effects  and  sampling  effects  simultaneously. 

Since  the  earliest  EEG  spectrum  was  obtained  by  manual 
calculation  (see  [Sm78]),  FFT  and  FWT  have  been  widely  util- 
ized for  EEG  spectral  estimation  [ GE7 5 1 [Ge80]  [Bo77]  [Ba79l 
[Ye72]  [B174]  [La80].  The  FFT  [C067],  whose  algorithms  are 
based  upon  that  of  decomposing  the  computation  of  the 
discrete  Fourier  transform  of  a sequence  of  length  1 N * into 
successively  smaller  discrete  Fourier  transforms  in  time  or 
in  frequency,  appreciably  decreases  the  computing  time  for 
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spectral  calculation.  At  present,  in  many  of  the  labora- 
tories [Ba79]  [Re77],  multichannel  FFT’s  are  used  routinely 
with  the  aid  of  hardware  FFT  processors  which  are  dedicated 
to  a micro  or  minicomputer. 

To  decrease  the  time-consuming  computational  procedure 
of  the  FFT,  other  orthogonal  functions  have  been  tried. 
Walsh  functions  [Ye72]  [La76]  is  one  of  them.  The  fast  Walsh 
transform  (FWT)  employs  Walsh  functions  in  the  series  ap- 
proximation of  the  original  signal.  Walsh  functions  [B174] 
are  a set  of  periodic  orthogonal  functions  which  are  binary 
in  nature  and  appear  to  be  a squared-off  version  of  the 
sinusoidal  functions.  These  binary  waveforms  make  the  fast 
Walsh  transform  very  attractive  for  many  kinds  of  signal 
processing  applications.  Because  only  two  function  values, 
+1  and  -1 , are  used,  the  need  for  multiplication  in  execu- 
tion of  the  FWT  is  obviated  by  taking  advantage  of  product 
identity  [Dz84 ] . 

However,  opinions  on  the  efficacy  of  using  the  FWT  in 
place  of  the  FFT  for  power  spectrum  estimation  are  split, 
advocates  versus  critics.  The  debate  centers  on  the  question 
of  whether  original  signal  information  can  be  obtained  as 
much  as  with  that  of  Fourier  analysis.  It  has  been  claimed 
that  Fourier  analysis  is  relevant  in  sinusoidally  based 
waveforms  and  the  Walsh  transform  is  more  appropriate  in 
cases  where  the  signal  contains  sharp  discontinuities  and  a 
limited  number  of  amplitude  levels  [Ja8l],  In  the  analysis 
of  the  truncation  and  roundoff  errors  associated  with  the 
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use  of  Fourier  and  Walsh  series,  it  is  claimed  that  for  long 
samples  of  smooth  signals,  far  more  terms  are  required  in 
the  Walsh  series  representation  and  greater  accuracy  is  re- 
quired of  their  coefficients  for  a given  root  mean  square 
total  error.  Even  for  discontinuous  signals  the  Walsh  series 
may  require  substantially  more  terms,  thus  counterbalancing 
the  computational  advantage  of  the  fast  Walsh  transform 
[B17M] . 

On  the  other  hand  it  is  reported  that  depending  on  the 
application,  useful  features  can  be  obtained  using  Walsh 
transform.  Though  the  discrimination  between  sleep  stages 
by  Walsh  power  spectra  is  not  as  good  as  the  Fourier  spectra 
for  sinusoidal  like  EEG  data  due  to  diffusion  throughout  the 
frequency  range,  the  activity  present  in  the  frequency  spec- 
tra was  evident  in  the  Walsh  spectra  plot  [Ye72],  In  an  EEG 
epoch  data  classification  problem  [La80]  with  predefined 
spectral  features,  features  selected  from  the  spectrum  of 
sleep  EEG  data  were  compared  to  corresponding  Fourier 
features.  Each  feature  set  was  used  to  classify  the  data  us- 
ing a minimum-distance  clustering  algorithm.  It  was  shown 
that  the  Walsh  spectral  features  classify  the  data  in  much 
the  same  way  as  the  Fourier  spectral  features.  It  is  claimed 
that  Walsh  spectral  features  in  place  of  Fourier  spectral 
features  enable  one  to  take  advantage  of  the  vast  computa- 
tional superiority  of  the  fast  Walsh  transform  over  the  fast 
Fourier  transform.  In  the  application  of  EEG  power  spectral 
estimation  to  on-line  monitoring,  in  which  it  is  desired  to 
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track  or  categorize  a subject's  state  versus  time,  the  FFT 
and  the  FWT  were  examined.  The  results  shows  that  when  the 
EEG  is  modeled  as  a Gaussian,  wide-sense  stationary  random 
process,  The  FWT  can  provide  useful  insight  into  performance 
with  actual  EEG  data,  with  care  in  interpretation  [Sm8l], 
In  the  statistical  evaluation  of  the  quantitative  and  dynam- 
ic differences  between  the  results  of  Walsh  and  Fourier 
analysis  of  the  EEG  done  for  the  purpose  of  characterizing 
the  effects  of  anesthesia,  the  efficacy  of  the  FWT  has  been 
reconfirmed  [Dz84].  However,  for  the  most  effective  applica- 
tions one  should  select  a proper  method  of  spectral  estima- 
tion and  the  interpretation. 

Dg_tg  M^els 

If  the  signal  is  locally  stationary,  it  may  be  more  ap- 
propriate to  derive  a model  of  a data  structure  which  has 
the  same  statistical  properties  as  the  signal  [Ch8l],  The 
signal  is  modeled  as  a linear  combination  of  its  past  output 
values  and  present  and  past  values  of  a hypothetical  input 
to  a system  whose  output  is  the  given  signal,  i.e.,  a linear 
system  output  with  a white  noise  input.  Because  the  signal 
is  predictable  from  the  past  and  present  values  of  inputs 
and  outputs,  increased  spectrum  resolution  can  be  obtained 
by  extrapolating  the  autocorrelation  function  beyond  the 
data  limit. 

There  are  three  basic  data  models,  autoregressive  (AR), 
moving  average  (MA),  and  autoregressive  moving  average 
(ARMA).  The  MA  model  and  AR  model  are  the  special  types  of 
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an  ARMA  model.  The  AR  model  output  is  given  as  a linear 
combination  of  past  outputs  and  present  input.  The  MA  model 
output  is  given  only  as  a linear  combination  of  past  inputs. 
In  the  data  model  there  are  two  problems  in  system  function 
approximation,  the  determination  of  the  order  and  coeffi- 
cients of  the  system  function.  In  general,  assuming  either 
the  deterministic  signal  or  random  signal  input,  the  method 
of  least  square  error  is  used  for  the  parameter  evaluation 
of  the  model.  Linear  prediction  [Ma76],  maximum  entropy 
[U175]  and  maximum  likelihood  [Ca67l  [Ca69l  methods  have 
been  used  widely  for  the  parameter  evaluation.  In  the  max- 
imum entropy  and  linear  prediction  methods,  the  autocorrela- 
tion function  is  extrapolated  in  such  a way  that  the  entropy 
of  the  corresponding  probability  density  function  is  maxim- 
ized in  each  step  of  the  extrapolation  of  the  autocorrela- 
tion function  [Bo71],  and  the  expected  mean  square  error 
between  the  next  sample  and  the  predicted  value  is  minimized 
[Ma76],  In  the  maximum  likelihood  method  the  input  is  as- 
sumed as  white  Gaussian  noise  and  output  is  the  maximum 
likelihood  estimate  of  data  sequence.  The  relationship 
between  maximum  entropy  spectra  and  maximum  likelihood  spec- 
tra was  shown  [Bu72]  as  the  reciprocal  of  the  maximum  likel- 
ihood spectrum  which  is  equal  to  the  average  of  the  recipro- 
cal of  the  maximum  entropy  spectra. 

A couple  of  criteria  have  been  used  to  determine  the 
order  of  the  data  model.  It  is  known  that  if  the  input  sig- 
nal is  a white  Guassian  noise  then  the  error  of  the  spectral 
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estimation  decreases  as  the  order  of  the  filter  increases. 
For  practical  purposes  the  minimum  order  of  the  AR  filter  is 
obtained  using  threshold  test 
1-Vp+1/Vp  < e 

where  Vp  is  the  normalized  error  of  the  pth  order  filter 
[Ma76].  Another  criteria  for  determining  the  order  of  the 
filter  is  Final  Prediction  error  [Ak74],  It  is  based  on  the 
fact  that  in  the  model-fitting  multiple  decision  criteria 
has  to  be  used  instead  of  the  hypothesis-testing  procedures. 
It  has  been  utilized  to  determine  the  order  of  the  filter 
for  the  case  of  autoregressive  signal  plus  noise  [To75] 
[To76]  . 

In  a method  of  fitting  an  AR  model  to  the  EEG  [Fe71]  a 
second  order  AR  model  is  initially  set  up.  After  calculating 
the  autocorrelation  function  and  the  power  spectral  density 
of  the  original  record  and  the  AR  model,  the  order  of  the 
filter  is  adjusted  to  satisfy  the  minimal  error  condition. 
The  results  show  that  the  optimal  epoch  period  for  the  EEG 
generator  is  20  to  30  seconds.  For  the  feature  extraction 
from  the  EEG  by  adaptive  segmentation  [Bo77]  [Pr77],  the 
linear  predictive  filter  was  used  to  split  EEG  record  into 
elementary  patterns  called  segments  and  transients.  Ap- 
propriate features,  representing  power  spectra  and  the  time 
structure  of  the  signal,  were  then  extracted  and  finally 
combined  into  a feature  set  representing  the  EEG  as  a whole. 
With  the  aid  of  second  order  differentiator  the  LP  filters 
could  give  a useful  degree  of  cluster  separation  between  the 
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epileptic  transients  and  background  activity  [Bi78b],  In 
the  hierarchical  model  of  the  EEG  signal  [Sa80] , a piecewise 
stationary  autoregressive  model  was  used  for  segment  clas- 
sification using  the  parameters  of  the  AR  model.  Though  it 
is  generally  understood  that  an  autoregressive  model  pro- 
vides a good  way  of  spectral  estimation  when  high  resolution 
is  needed  with  short  data  records  [Ma76],  one  must  be  care- 
ful in  choosing  the  order  of  an  autoregressive  model  and  the 
sampling  rate  for  the  statistical  characteristic  representa- 
tion of  the  signal. 

S££j£ AD^l^sis 

In  addition  to  an  appropriate  spectral  estimation,  an 
appropriate  spectral  interpretation  such  as  dominant  feature 
selection  is  the  key  to  spectral  analysis.  The  duration  and 
location  of  the  data  window  can  cause  improper  results  of 
the  analysis.  Adaptive  segmentation  has  been  utilized  to 
locate  the  spectral  window  properly.  Transient  activity 
could  be  detected  using  this  method.  Spectral  parameter 
analysis,  compressed  spectral  array  and  spatiotemporal 
analysis  have  been  utilized  for  proper  interpretation  of  the 
EEG  spectrum.  The  change  of  the  estimated  spectrum  over 
time  can  be  understood  in  terms  of  wave  parameters. 

S^ectxaJ,  Parameter  Arialyji,j 

Spectral  parameter  analysis  (SPA)  is  based  on  the  as- 
sumption that  the  spectral  parameters  describing  the  signal 
do  not  vary  during  the  analysis  epoch  which  is  usually  10  or 
20  sec.  In  other  words,  the  signal  is  regarded  as  a locally 
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stationary  random  process  during  such  a period  [Is75a] 
[Is75b].  SPA  provides  an  objective  description  of  the  dis- 
tribution of  spectral  power  in  the  EEG  signal  in  the  form  of 
a rationale  spectrum  with  no  more  than  8 parameters.  The 
spectrum  is  divided  into  1-3  components  that  are  described 
by  frequencies  and  power  parameters  such  as  bandwidth,  peak 
frequency  and  power.  These  spectral  parameters  are  deter- 
mined from  the  estimated  power  spectrum. 

SPA  was  used  for  describing  the  normal  EEG.  The  varia- 
tions with  time  in  the  center  frequency  of  the  alpha  activi- 
ty and  the  relation  between  the  bandwidths  of  delta  and  al- 
pha activity  were  studied [We7 1 1 . SPA  forms  the  basis  of  an 
EEG  signal  generator  [Ze69].  The  method  involves  determina- 
tion of  the  parameters  for  a data  model  of  the  EEG  signal. 
From  the  experimental  data  the  parameters  are  estimated  by 
first  computing  the  first  100  correlation  coefficients  and 
then  applying  an  AR  or  ARMA  model  to  generate  a stationary 
discrete  sequence.  The  parameter  estimation  procedure  was 
tried  on  a few  records  and  it  was  found  that  typically  a 
fifth  order  system  suffices  to  give  an  accuracy  of  about  5 
percent  for  the  correlation  coefficients  which  describe  the 
stationary  parts  of  EEG  signals. 

Adaptive  Segmentation 

Adaptive  segmentation  procedure  is  very  close  to  quan- 
tifying the  procedure  of  visual  analysis  of  the  EEG.  This 
method  was  used  for  transient  wave  detection  [Bo77]  [Mi79l 
[Pr77J.  For  adaptive  segmentation  it  is  assumed  that  the 
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EEG  is  a piecewise  stationary  time  series.  For  the  feature 
extraction  the  adaptive  segmentation  procedure  with  the  aid 
of  AR  model  was  used  [B077].  This  procedure  was  developed 
to  detect  the  segment  boundaries  and  locate  transients,  and 
to  represent  the  information  in  the  segments  in  terms  of  a 
set  of  parameters  specifying  their  power  spectra.  Here, 
"piecewise  stationary"  means  that  within  a segment,  the 
short-time  spectral  estimate  does  not  change  appreciably 
with  time.  Transients  are  short  time  nonstationarities . A 
segment  was  specified  by  its  length  and  power  spectra  esti- 
mate, a transient  by  a set  of  grapho-eleraents  and  time  of 
occurrence.  Though  it  showed  good  results  for  the  limited 
data,  the  selection  of  an  optimal  threshold  level,  which  is 
very  sensitive  and  dependent  on  the  subject,  is  not  general- 
ized . 

For  an  adaptive  segmentation  of  the  EEG  records  [Pr77J 
from  a two  second  data  record  an  AR  model  was  derived.  By 
passing  an  EEG  signal  through  this  filter  the  spectral  error 
is  measured.  When  the  spectral  error  exceeds  the  threshold  a 
new  AR  model  is  derived  and  the  transient  is  detected.  It  is 
reported  that  the  threshold  setting  for  the  detection  of  the 
fast  transients  is  more  difficult  and  needs  further  valida- 
tion on  more  material  if  it  is  to  be  used  for  clinical  pur- 
poses. 

Comores^e^  3££c tjjj].  Array 

The  compressed  spectral  array  (CSA)  is  a graphical 
display  method  that  can  compress  and  simplify  complex  data. 
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This  emerges  from  the  multiple  scalp  electrodes  over  consid- 
erable intervals  of  time.  How  and  in  what  location  the  EEG 
changes  in  frequency  and  amplitude  is  observed  easily  with 
the  CPA,  using  a hidden  line  suppression  algorithm  [Fo82] 
and  fast  Fourier  transform.  It  is  displayed  in  terms  of 
specttral  power  in  the  frequency  range  of  0 to  16Hz.  CSA  has 
been  practically  useful  in  monitoring  the  anesthesia  level 
of  the  patients  in  the  operating  room  [Bi73].  By  displaying 
the  spectra  using  a general  computer,  hemisphere  asymmetries 
can  also  be  detected.  However,  in  multichannel  application, 
because  of  excess  data  in  a short  period  there  are  difficul- 
ties in  the  interpretation  of  consecutive  spectral  arrays. 

It  is  important  to  recognize  a spatiotempor al  pattern 
since  the  occurrence  of  EEG  pattern  in  time  and  space  cannot 
be  separated.  In  a spatiotempor al  analysis  it  is  assumed 
that  there  are  a few  widely  spaced  EEG  signal  generators  in 
the  cortex  that  have  the  same  frequency  with  different 
phases  [J069].  That  is,  EEG  of  the  scalp  is  regarded  as  a 
traveling  wave  . As  a technique  for  displaying  the  EEG  spa- 
tiotemporal  characteristics,  a spatiotemporal  map  shows  how 
the  potential  distribution  along  the  lines  of  electrodes 
changes  with  time  [Re68].  A contour  map  presents  the  dis- 
tribution of  potential  over  an  area  of  the  head  at  one  in- 
stant of  time.  It  is  derived  from  a two  dimensional  elec- 
trodes array.  Using  interpolation,  points  of  the  same  poten- 
tial are  joined  together.  Spatiotemporal  analysis  is  used 
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to  locate  the  origin  of  abnormal  EEG  waves  [Br75l.  The  mul- 
tidimensional filtering  which  utilized  the  multidimensional 
Fourier  transform,  is  used  to  decompose  composite  wavefronts 
[Ch77l.  The  method  makes  it  possible  to  track  selected  wave 
propagations  across  the  array  of  scalp  electrodes. 

W3xe.fp.crg  PxCPIBPiPX  iUXlXSiP 

The  information  of  the  picture-like  EEG  data  can  be  un- 
derstood and  transformed  by  visual  recognition  better  than 
with  any  other  method.  Individual  waveform  recognition  or 
detection  is  the  key  to  waveform  analysis.  Once  a certain 
waveform  is  recognized,  it  is  interpreted  by  considering  the 
context  of  the  same  channel  or  the  other  channels.  The  gen- 
eral characterization  of  an  amplitude-time  pattern  of  each 
epoch  and  the  detection  of  the  individual  wave  patterns 
which  can  be  defined  in  advance  are  necessary  for  a qualita- 
tive description  of  the  sleep  EEG  data.  Several  kinds  of 
methods  have  been  utilized  to  describe  the  EEG  activity  in 
time  domain.  Amplitude  and  period  of  the  waveform  are  the 
basic  parameters  which  are  used  to  describe  waveform  activi- 
ties in  ways  similar  to  the  visual  evaluation.  The  first 
and  second  order  derivatives  of  the  waveform  have  also  been 
utilized  to  check  the  sharpness  and  the  local  peaks  of  the 
wave  as  the  description  parameters. 

Each  epoch's  activity  is  represented  by  the  total 
amount  of  the  absolute  amplitude  of  the  EEG.  In  a sleep 
staging  application  [Ag67l  [Ag73l,  it  was  reported  that  an 
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amplitude  analysis  of  the  EEG  alone  was  not  sufficient  for 
the  sleep  stage  determination  because  : the  contingency 
between  the  stage  of  sleep  and  the  electrogenesis  level  was 
not  a consistent  phenomenon,  and  only  sleep  stage  4 has  a 
level  which  can  be  clearly  distinguished  from  other  states. 
That  is,  the  electrogenesis  level  is  subject  dependent  and 
the  detection  of  the  specified  waveforms,  which  are  not 
identified  in  the  integrated  EEG,  is  essential  for  the  sleep 
staging. 

In  order  to  describe  an  EEG  data  quantitatively  three 
parameters  of  an  EEG  signal  known  as  slope  descriptors  were 
defined  in  the  time  domain.  They  are  equivalent  to  the 
zeroth,  second  and  fourth  moments  of  the  normalized  power 
spectral  density.  These  parameters  are  defined  as  follows: 

a)  activity,  the  variance  of  the  amplitude  or  mean 
power  in  the  epoch; 

b)  mobility,  the  average  power  of  the  normalized 
derivative  in  the  epoch; 

c)  complexity,  the  average  power  of  the  normalized 
second  derivative  in  the  epoch. 

Using  these  parameters  the  basic  properties  of  the  spectrum 
can  be  derived  without  any  frequency  analysis,  thus  requir- 
ing less  computation  time  [ H j7 0 ] . It  was  reported  [Bi78a] 
that  spectral  analysis  proved  more  reliable  than  Hjorth’s 
analysis  for  both  discriminating  and  predicting  the  site  of 
cerebral  pathology  between  patients  with  localized  pathology 
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and  normal  control.  This  method  has  not  come  into  general 
use . 

Zgri-siJ 

The  base  line  crossing  (zero  crossing)  of  the  EEG  and 
the  first  or  second  derivatives  of  the  EEG  are  determined  in 
order  to  select  a specified  frequency  band  activity  [It69l 
[Ro70]  [Sm74].  A high  frequency  component  superimposed  on  a 
slow  wave  can  be  detected  using  this  method.  The  relation- 
ship between  frequency  and  amplitude  gives  the  various 
parameters  which  describe  the  characteristics  of  the  EEG 
[ 1 16 9 1 [Fr69l.  It  was  shown  that  counting  the  number  of  zero 
crossings  over  an  epoch  is  effectively  equivalent  to  the 
more  complex  procedures  of  computing  zero,  second  and  fourth 
moments  of  the  power  spectral  density  [Sa71].  A comparison 
of  spectral  analysis  with  period-amplitude  analysis  (which 
utilizes  a period  and  absolute  amplitude  histogram),  applied 
to  the  narrow  band  EEG  activity  [Kt8l],  showed  that  while 
the  power  spectrum  efficiently  quantifies  the  overall  power 
trend  in  the  EEG  data,  period  amplitude  analysis  offers  more 
resolution  than  the  power  spectrum  in  detecting  electro- 
graphic details  in  amplitude  and  incidence  within  relatively 
narrow  frequency  bands, 
to eXSJCIS  ££iec_tj.or) 

This  method  [Go74]  [Sm75]  is  based  on  the  human  visual 
detection  of  the  waveform.  Specific  waveforms  which  are 
used  for  sleep  analysis  are  detected  by  checking  the  zero 
crossing  and  peak  amplitude  of  the  EEG  waveform.  It  consists 
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of  analog  filters,  amplitude  detectors,  period  analyzers  and 
wave  pattern  recognizers.  A filter  is  provided  for  the 
proper  operation  of  the  zero  crossing  detector.  Fast  activi- 
ty superimposed  on  the  big  slow  wave  has  to  be  filtered 
first  or  the  zero  crossing  detector  can  give  the  slow  wave 
period  instead  of  fast  wave  period.  The  pattern  recognizer 
selects  the  specific  wave  pattern  using  minimum  amplitude 
criteria  and  wave  criteria  [Sm78].  Over  92 % agreement  was 
obtained  between  human  and  automatic  detection  of  sleep 
spindle  waveforms  [Jo76]  C S i7 61. 

Pattern  £^Sfi£J3iiion 

If  the  actual  waveform  of  an  EEG  is  not  required  to  ob- 
tain the  information  required  from  a particular  measurement 
or  experiment,  mathematical  pattern  recognition  procedures 
can  be  used  to  extract  information.  There  are  two  different 
cases  in  pattern  recognition  problem  problems.  One  is  clus- 
tering and  the  other  is  a classification  problem.  When  pri- 
or knowledge  is  not  available,  the  problem  is  to  extract 
certain  features  which  are  successfully  able  to  classify  the 
given  data  set  into  subsets  for  proper  classification.  If  a 
prior  knowledge  is  not  available  then  the  problem  is  to  form 
clusters  which  share  certain  common  characteristics  [De77]. 

Usually  in  clinical  EEG  data  analysis,  a measurement 
data  set  is  classified  into  subsets  which  are  mutually  ex- 
clusive. The  pattern  recognition  rules  or  features  are 
derived  from  the  heuristic  rules  which  are  obtained  from  the 
observation  of  data  set  and  clinical  diagnoses,  or  from 
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previously  selected  training  data  set.  Many  features  from 
the  frequency  domain  and  the  time  domain  are  utilized  in  EEG 
analysis.  The  power  spectrum  obtained  from  FFT,  FWT  or  data 
models  is  used  as  basic  information  for  feature  extraction 
in  frequency  domain.  Wave  parameters,  such  as  Hjorth  param- 
eters [Hj70],  wave  period,  slope  or  time  of  occurrence  are 
used  for  time  domain  analysis.  Most  of  the  pattern  recog- 
nition techniques  have  been  applied  [Ge80]  [Mc8l],  including 
various  decision  functions,  linear  classifier,  quadratic 
classifier  and  syntactic  methods. 

A correlation  technique  was  used  for  recognition  of 
spontaneous  EEG  patterns  in  time  domain.  Correlations  were 
obtained  using  two  different  moving  templates  [We72],  The 
correlation  classifier  and  maximum  likelihood  quadratic 
classifier  were  compared  to  classify  single  evoked  poten- 
tials from  four  different  visual  stimuli.  Forward  sequential 
feature  selection  was  employed  to  reduce  data  dimensionali- 
ty. It  was  noted  that  if  the  covariance  matrices  of  the  nor- 
mal distributions  are  not  equal,  i.e.,  noise  is  correlated 
sample  to  sample,  the  quadratic  discriminant  function  gives 
better  performance  than  the  linear  one  [Se79].  Iteratively 
determined  discriminant  functions  were  used  in  a decision 
tree  for  all  night  sleep  staging  [Ma72].  Separate  systems 
were  designed  for  high  and  low  alpha  subjects  and  heuristic 
rules  were  also  included  in  the  decision  tree.  Linear  step- 
wise discriminant  analysis  were  also  applied  to  sleep  stag- 
ing [La70].  Forward  sequential  feature  selection  is  also 
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often  used  to  select  features  in  a sequential  manner.  To 
date,  no  satisfactory  techniques  have  been  found  for  an  op- 
timum feature  selection  for  the  general  classification  prob- 
lem. A multivariate  parametric  classifier  is  compared  with 
a K-nearest  neighborhood  classifier  [Yu80].  It  was  reported 
that  the  optimized  K-nearest  neighborhood  classifier  gave 
slightly  better  results  than  the  parametric  one,  the  over 
all  level  of  success  was  below  70  percent. 

The  syntactic  pattern  recognition  analysis  was  also 
used  [Gi79].  Training  of  a classifier  was  conducted  with 
well  defined  data  epochs,  selected  by  electroencephalogra- 
phers.  A standard  token  set  was  obtained  from  these  epochs. 
Preparse  contextual  analysis  was  used  to  avoid  ambiguities 
caused  by  the  training  set  of  broad  spectral  peaks.  Post- 
parse contextual  analysis  was  also  used  to  evaluate  cross 
channel  differences.  Although  this  method  gave  over  80  per- 
cent correctness,  the  selection  of  primitives,  and  produc- 
tion rules  which  govern  the  merging  of  primitives  were  not 
generalized. 

Each  method  has  advantages  and  limitations  in  applica- 
tion to  sleep  EEG  analysis.  In  spectral  analysis,  the  time 
domain  information  such  as  the  duration  of  specific  wave 
bursts  or  the  number  of  the  occurrences,  which  is  essential 
for  the  clinical  EEG  analysis,  is  not  well  provided.  Though 
the  waveform  analysis  can  provide  effective  time  domain  in- 
formation, the  overall  power  trend  in  each  epoch  or  contex- 
tual information  of  the  specific  waveforms  are  not  well 
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provided.  In  pattern  recognition  analysis,  because  of  the 
variety  in  feature  vectors,  an  individual  specific  event 
such  as  transient  peak  or  a specified  waveforms  cannot  be 
identified  properly.  It  is  clear  that  none  of  the  methods 
described  to  this  point  is  good  enough  to  perform  the  com- 
plex procedure  of  the  sleep  EEG  analysis. 

In  conclusion,  in  computerized  sleep  analysis,  indivi- 
dual waveforms  such  as  alpha,  beta,  delta,  sigma,  etc.  can 
be  detected,  but  certain  waveforms  and  sleep  states,  such  as 
REM,  slow  wave  or  K-complex  can  be  recognized  and  interpret- 
ed only  by  considering  multiple  channel  activities  and/or 
the  previous  and  following  epoch  activities.  That  is,  after 
detecting  an  individual  wave,  a contextual  analysis  of  mul- 
tichannel EEG/EOG  signals  is  required. 


CHAPTER  3 

EXPERT  SYSTEM  APPROACH 


The  previous  chapter  discussed  the  advantages  and  the 
difficulties  or  weaknesses  of  the  conventional  methods  for 
sleep  EEG  data  analysis.  We  notice  that  the  procedure  for 
the  sleep  EEG  data  analysis  needs  two  processes-namely  the 
visual  perception  of  the  waveforms  and  the  contextual  under- 
standing of  each  waveform. 

For  visual  perception  flexible  waveform  criteria  which 
can  be  adjusted  depending  on  the  subject  or  required 
analysis  are  needed.  For  instance,  the  amplitude  threshold 
value  for  a certain  waveform  detection  can  be  different 
depending  on  the  subject.  For  the  contextual  understanding 
of  each  waveform  or  contextual  analysis  of  each  channel  ac- 
tivity, a set  of  rules  is  required.  These  rules  form  the 
domain  specific  knowledge  base  for  EEG  signal  analysis  using 
all  the  channel  information  and  subject  dependent  data.  The 
channel  activity  for  an  epoch  is  described  using  these 
rules.  The  flexibility  of  the  waveform  detection  and  expert 
rules  for  contextual  understanding  of  the  waveforms  are  the 
key  elements  to  be  encoded  in  the  present  EEG  analysis  sys- 
tem. 

Recently,  in  the  field  of  Artificial  Intelligence, 
various  kinds  of  expert  systems  were  developed  to  solve  the 
problems  which  require  the  knowledge  base  of  a human  expert 
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and  a deduction  mechanism  for  decision  making  [Ba82a] 
[Ba82b]  [ Co8 1 ] [Mi8l]  [Ge83]  [St82]  [We84]  [Mc82].  In  the 
following  section  Expert  Systems  are  briefly  discussed  to 
show  their  present  status  and  how  and  to  what  extent  such 
techniques  can  be  utilized  for  the  sleep  EEG  data  analysis. 

BjctoP-yUii 

Most  of  the  early  work  in  AI  research  centered  around 
an  attempt  to  find  a general  (domain  independent)  problem- 
solving strategy  (GPS)  in  a belief  that  a few  reasoning 
mechanisms  with  powerful  computers  can  give  the  solutions  to 
domain-independent  problems.  Though  this  GPS  research  demon- 
strated the  important  role  of  the  fundamental  techniques 
such  as  various  search  methods,  knowledge  representation 
methods,  and  problem  decomposition  into  subgoals,  ultimately 
it  led  to  the  view  that  they  could  not  offer  sufficient 
power  to  solve  most  complex  problems  [Ne69].  Many  of  the 
researchers  have  instead  focused  on  simple  or  single  domain 
problems  which  require  incorporation  of  large  amounts  of 
domain  specific  knowledge.  These  are  the  Expert  Systems. 

Expert  systems  deal  with  various  problems  which  fall 
into  one  of  the  following  categories:  data  interpretation, 
prediction,  diagnosis,  monitoring,  planning,  design,  debug- 
ging, repair,  instruction,  control,  etc.  [Ha83b],  Expert 
systems  solve  domain  specific  problems  at  the  level  of  human 
expertise  utilizing  accumulated  expert  knowledge,  comprehen- 
sive knowledge  representation  and  inference  schemes.  Expert 
knowledge  provides  the  key  to  expert  performance,  while 
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comprehensive  knowledge  representation  and  inference  schemes 
provide  the  mechanism  for  its  use  [Fe77].  Most  of  the  work 
in  building  an  expert  system  involves  capturing  the  expert 
knowledge  and  then  forming  a knowledge  structure  for  the 
deduction  mechanism  or  inference  schemes  to  use  that  expert 
knowledge  efficiently. 

Various  knowledge  representation  techniques  and  infer- 
ence schemes  [Wi79)  [Ri83J  [We84]  [Ni80]  have  been  developed 
for  different  domain  problems.  The  structure  of  the  infer- 
ence scheme  is  usually  closely  related  to  the  method  of 
knowledge  representation.  Once  the  rule  language  has  been 
defined,  the  inference  engine  can  be  built  in  terms  of  the 
general  framework.  Rules  for  the  domain-specific  knowledge 
are  written  in  the  previously  defined  language  format  or 
grammar.  In  an  ideal  expert  system,  the  inference  engine 
and  knowledge  base  are  separated  for  system  development.  For 
instance,  if  a certain  domain  knowledge  is  represented  using 
production  rules,  the  knowledge  base  can  be  changed  without 
affecting  the  inference  engine.  R1  (an  expert  system  for 
configuring  computers)  [Mc82]  is  a good  example.  Though  the 
inference  engine  and  the  knowledge  base  are  written  as 
rules,  the  knowledge  base  is  changed  without  affecting  the 
control  strategy  of  the  inference  engine. 

Recently,  some  of  the  expert  systems  and  general- 
purpose  representation  languages  have  been  developed  as  ex- 
pert system  design  tools.  These  expert  systems  evolved  from 
their  previous  expert  systems,  and  the  specific  languages 


have  been  used  or  developed  for  their  own  expert  systems 
[Ha82b].  Although  each  of  these  expert  systems  has  been 
generalized  to  adapt  different  domain  problems,  their  per- 
formances are  limited  by  the  previously  defined  rule 
language  and  the  control  structure  embodied  in  the  inference 
engine.  Although  the  general  purpose  representation 
language  has  fewer  limitations,  the  appropriate  rule 
language  and  control  structure  have  to  be  defined  in 
correspondence  with  appropriate  knowledge  structure  to 
represent  their  domain  specific  knowledge. 

-Same  . 

Several  expert  systems  are  reviewed  in  this  section  to 
provide  an  insight  into  the  knowledge  representation  tech- 
niques and  inference  schemes  that  can  be  utilized  in  build- 
ing expert  systems.  These  systems  (CASNET,  MYCIN,  PIP,  IN- 
TERNIST, HEARSAY-II)  which  were  developed  in  the  early  age 
of  AI  for  medical  consultations  are  basically  based  on  the 
production  system  and  are  classification  models.  Though  the 
following  expert  systems  could  not  find  their  way  into  every 
day  use,  they  demonstrate  that,  depending  on  the  problem 
domain,  an  appropriate  knowledge  representation  technique 
and  inference  scheme  should  be  utilized. 

CjSNET/gjaucprgj 

A causal  associational  network  [Ku80]  was  developed  to 
represent  knowledge  of  reasoning  about  diagnosis,  prognosis 
and  treatment  selection  in  opthalmology . Production  rules 
were  used  to  represent  the  relationships  between  observation 
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and  diagnosis  in  a particular  type  of  semantic  network.  The 
representation  of  the  domain  knowledge  is  composed  of  three 
layers:  observation,  pathophysiological  state,  and  disease 
category.  Each  observation  is  mapped  into  the  corresponding 
pathophysiological  state  through  the  implication  links.  The 
classification  links  connect  each  pathophysiological  state 
to  a disease  with  a one  to  many  relation.  In  pathophysio- 
logical states,  the  causal  links  are  used  to  relate  cause 
and  effect.  The  causal  relations,  with  associated  degrees 
of  strength,  express  the  mechanism  of  a disease  and  their 
modifications  under  various  regimens  of  treatment  [Ku80] . 
With  causal  ordering  information,  the  most  likely  cause  for 
a particular  set  of  patient  data  can  be  easily  deduced  by 
traversing  the  causal  network. 

The  reasoning  control  strategy  of  CASNET  starts  in 
event-driven  fashion.  Utilizing  the  incoming  clinical  data 
and  a thresholding  evaluation  mechanism,  patient-specific 
interpretation  models  or  hypotheses  are  formed  for  confirmed 
and  undetermined  states.  At  these  states  the  causal  model 
is  used  to  constrain  the  search  space  for  possible  hy- 
pothesis. The  next  question  choice  is  "hypothesis-driven". 
When  all  the  available  data  are  accumulated,  the  system 
makes  a final  decision  about  diagnostics,  prognostics,  and 
treatment  categories  in  a purely  deterministic  fashion. 
However,  the  CASNET/Glaucoma  model  never  became  widely  used. 
It  was  pointed  out  that  the  problem  domain  was  related  in 
scope  to  many  additional  opthalmological  problems  which  were 
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not  covered  fully  in  the  original  model.  Also,  many  consul- 
tation problems  did  not  easily  fit  into  a causal  framework. 

EXPERT  as  an  expert  system  design  tool  evolved  from 
CASNET.  This  production  rule  language  maps  findings  onto 
hypothesis  and  solves  problems  using  a simple  classification 
model.  In  this  system  the  production  rules  are  classified 
in  terms  of  the  three  types  of  logical  relationships  between 
findings  and  hypothesis:  FF  (finding  to  finding),  FH  (find- 
ing to  hypothesis),  and  HH  (hypothesis  to  hypothesis)  rules. 
FF  rules  are  used  to  establish  local  control  over  the 
sequencing  of  the  questions  i.e.  to  specify  deterministic 
logical  relations  between  findings  whose  truth  value  is  al- 
ready found.  FH  rules  are  used  to  assign  a certainty  factor 
to  a hypothesis  from  the  findings.  FH  and  HH  rules  are  used 
to  infer  hypothesis  from  findings  and  hypothesis  respective- 
ly. Using  these  simple  logical  relations  the  reasoning  rules 
and  expert  knowledge  are  described  in  production  rules. 
While  these  simple  rule  classifications  make  EXPERT  easy  to 
use,  it  is  limited  in  the  application  to  the  different 
domain  problems  which  have  different  knowledge  structures. 
However,  EXPERT  has  been  used  in  numerous  applications  in- 
cluding medicine,  oil  well-log  analysis,  laboratory  instru- 
mentation, and  computer  fault  diagnosis  [We84]. 

MYCIN  [Sh73J  [Sh75l  C Sh76 ] [Sh79]  was  developed  to 
determine  the  likely  identities  of  pathogens  in  patients 
with  infections,  and  to  assist  physicians  in  the  selection 
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of  a therapeutic  regimen  appropriate  for  opposing  the  organ- 
ism under  consideration.  Production  rules  with  certainty 
factors  were  used  to  represent  the  control  strategy  and 
domain  specific  knowledge.  The  certainty  factors,  which  are 
an  approximation  of  conditional  probabilities,  are  used  to 
model  the  inexact  reasoning  process  of  medical  experts.  The 
certainty  factor  is  defined  as  MB  - MD  (measure  of  belief 
and  disbelief).  MB  and  MD  are  defined  in  terms  of  condi- 
tional and  a priori  probabilities  [Sh75l . The  certainty 
factor  relates  possible  conditions  to  associated  interpreta- 
tions with  weighted  relations  ranging  from  -1  to  +1 . In  the 
predicate  structure  (IF-THEN)  representation  of  the 
knowledge  the  antecedents  and  actions  are  expressed  with  the 
attribute-triple  pair  with  certainty  factor  or  attached  pro- 
cedures. The  ■ associative  triple  pair  consists  of  "attri- 
bute", "object”,  and  "value".  In  MYCIN  the  objects  in  the 
triple  pairs  are  the  domain  entities  which  are  organized 
into  a hierarchical  context  tree.  Each  finding  or  observa- 
tion is  related  to  the  value  of  the  triple  pair.  The  triple 
pairs  represent  medical  facts  and  hypotheses  or  related 
diseases.  This  uniformity  of  knowledge  representation  made 
it  system  possible  to  separate  domain  specific  knowledge  and 
control  strategy  in  this  system. 

As  with  most  production  systems,  a goal-directed  back- 
ward chaining  strategy  was  used.  To  find  the  highest  level 
goal,  such  as  treatments  for  all  the  infections  of  the  pa- 
tient, MYCIN  starts  by  testing  the  condition  of  the  highest 
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level  rule  against  available  data,  or  requests  data  from  the 
physician.  Every  possible  condition  is  checked  for  validity 
and  other  rules  can  be  inferred  recursively  at  this  stage, 
generating  subgoals  to  be  confirmed  with  findings.  MYCIN 
also  utilized  certainty  factors  in  the  deduction  process. 
It  is  a hierarchical  AND/OR  tree  with  certainty  factors, 
combining  uncertain  assertions  and  heuristic  cumulative 
function  for  decision  making  [Sh75] , using  a fuzzy  logic 
function.  By  looking  into  a stack  of  production  rules  used, 
MYCIN  can  answer  questions  about  decision  making  such  as  how 
and  why.  The  premises  and  the  conclusions  of  the  rules  are 
utilized  to  generate  a natural  language  like  answer. 

The  MYCIN  system  has  been  generalized  to  EMYCIN  (Essen- 
tial MYCIN)  [Me8l]  for  domain-independent  use.  EMYCIN  has  a 
rule  compiler  which  takes  advantage  of  knowledge  available 
from  static  analysis  of  the  rule  set  to  produce  a more  effi- 
cient representation.  It  constructs  a LISP  program  which 
eliminates  redundant  computation  of  a same  premise  clause  or 
conditions  in  different  rules.  Thus  the  rule  compiler  makes 
the  system  use  an  efficient  deductive  mechanism.  This  rule 
compiler  works  like  a rule  manager  to  make  rule  addition, 
debugging,  and  explanation  easier.  EMYCIN  has  been  used  to 
build  various  consultation  programs.  In  an  application 
[Br80]  of  EMYCIN,  some  problems  about  backward  chaining  and 
parallel  searching  with  a negative  certainty  factor  were  re- 
ported. With  the  backward  chaining  control  structure,  the 
grouping  of  questions  is  determined  by  the  order  of  the 
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appearance  of  the  parameters  or  variables  in  rule.  That  is, 
data  has  to  be  used  in  a serial  rather  than  parallel 
fashion.  And  in  parallel  search,  the  certainty  factor  has  to 
be  adjusted  to  be  dependent  between  subtrees  with  the  same 
node . 

PIP/P£es£ni  XE&SMSl 

PIP  was  developed  to  take  the  history  of  the  present 
illness  of  a patient  with  edema  [Pa76],  The  present  illness 
forms  a basis  upon  which  a physician  builds  his  diagnosis 
and  bases  his  treatment  schedule.  PIP  consists  of  four 
building  blocks:  patient  specific  data  (input  data),  super- 
visory program  (inference  engine),  short  term  memory  (work- 
ing memory,  hypothesis  under  testing),  and  associative 
memory  (knowledge  base).  The  supervisory  program  is  a con- 
trol program  which  guides  the  system  in  taking  the  present 
illness  and  controls  the  operation  of  various  subprocesses, 
such  as  selecting  questions  or  relevant  rule  selection  to 
build  a coherent  model  of  the  case.  The  short  term  memory 
is  the  site  in  which  the  patient’s  data  interact  with  the 
knowledge  base  in  associative  memory.  The  associative 
memory  as  a system  knowledge  base  is  organized  into  many 
frames.  Each  frame  is  a structure  with  a name  and  a number 
of  slots,  which  can  be  filled  by  various  properties,  logical 
and  semantic  relations,  and  associated  inference  rules. 
Frames  in  the  associative  memory  represent  the  diseases, 
clinical  states  or  physiological  states.  The  frames  are 
linked  into  a complex  network.  The  links  of  the  network 
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represent  the  kinds  of  association  between  the  prototypical 
findings,  such  as  "caused-by ” , "must  exclude",  "component 
of",  etc.. 

The  reasoning  starts  with  the  presentation  of  the 
patient’s  complaint,  which  is  characterized  by  the  supervi- 
sor utilizing  the  procedure  kept  in  the  knowledge  base  of 
long  term  memory.  The  characterized  patient’s  complaint  ac- 
tivates the  corresponding  rules.  Then  a hypothesis  is  gen- 
erated, moving  frames  from  the  associative  memory  to  the 
short  term  memory,  where  each  frame  recursively  induces  re- 
lated frames  in  the  network.  To  evaluate  each  frame  or  hy- 
pothesis a scoring  process  is  performed.  In  the  scoring 
process,  certainty  factors  contained  in  the  frame  are  used 
to  reflect  the  likelihood  that  various  clinical  findings  ex- 
ist in  the  given  disorder.  This  score  has  two  components:  a 
measure  of  fit,  i.e.,  matched  findings  to  expected  findings, 
and  the  ratio  of  the  number  of  findings  matched  to  the  total 
number  of  findings.  In  this  evaluation  process,  if  neces- 
sary, more  information  is  collected  to  fill  each  slot  of  the 
frame.  It  has  been  reported  [Sz79]  that  the  generalized 
scoring  function  and  the  threshold  value  for  termination 
procedure  caused  difficulties  in  generating  lines  of  the 
reasoning  process. 

JXXE&mSJ-J  / c jUmaidAant  in  ini&CJiai 

INTERNIST-I  [Mi82]  has  been  developed  to  aid  the  physi- 
cian with  a patient’s  workup  in  order  to  make  multiple  and 
complex  diagnosis  in  internal  medicine.  The  INTERNIST-I 
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attempts  to  form  an  appropriate  differential  diagnosis  in  a 
selected  group  of  observed  findings,  the  differential  diag- 
nosis of  which  forms  what  is  assumed  to  be  a mutually  closed 
set  of  diagnoses.  The  search  space  of  INTERNIST-I  is  very 
large  and  highly  structured  as  in  general  internal  medi- 
cine. The  individual  disease  is  a building  block  of  the 
knowledge  base.  For  each  diagnosis  entered  into  the  system, 
a disease  file  is  generated  to  include  findings  that  have 
been  reported  to  occur  in  association  with  disease. 

INTERNIST-I  has  knowledge  representation  of  weighted 
links  between  findings  and  diseases.  The  diseases  have 
hierarchical  structure  from  general  to  specific.  The 
"evoking-strength"  links  finding  to  diseases  and  "frequency" 
links  diseases  to  findings.  The  evoking-strength  links 
describe  how  strongly  the  findings  explain  the  diagnosis. 
The  frequency  is  an  estimation  of  how  often  patients  with 
the  disease  have  the  findings.  The  other  links  are  causal 
links  which  connect  diseases  to  diseases.  Each  finding  is 
assigned  a disease  independent  "import"  weight.  The  import 
is  the  global  importance  of  the  findings  or  manifestation, 
that  is,  the  extent  to  which  one  is  compelled  to  its  pres- 
ence in  any  patient.  Most  of  the  rules  are  expressed  as 
link  relations. 

The  reasoning  strategy  of  INTERNIST-I  begins  in  data 
driven  fashion.  After  partitioning  general  findings  into 
problem  areas,  a set  of  selected  diseases  are  tested  against 
the  partitioned  findings  to  obtain  highly  weighted  diseases. 
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Each  disease  is  scored  positively  and  negatively  based  on 
the  findings  found  and  findings  not  found  respectively. 
Heuristic  principles,  such  as  diagnosis  by  exclusion,  are 
employed  to  resolve  each  of  the  differential  diagnostics. 
To  control  the  number  of  induced  hypotheses,  different  stra- 
tegies are  used  depending  on  the  number  of  competing  hy- 
potheses. However,  the  most  highly  ranked  hypothesis  is 
selected  using  explained  findings.  Then  the  program  repeats 
the  cycle  with  unexplained  findings  until  all  the  findings 
have  been  accounted  for.  Several  problems  were  reported, 
including  the  inability  to  construct  differential  diagnoses 
spanning  multiple  problem  areas,  the  inability  to  explain 
its  reasoning,  inability  to  reason  anatomically  or  tem- 
porarily, and  the  occasional  attribution  of  findings  to  im- 
proper causes.  The  successor  to  this  system  CADUCES  is 
still  being  developed  [Mi82],  (see  [We84]). 

& -]j£ S£S$£IP 

HEARSAY-II  was  developed  to  understand  connected  speech 
[Er80],  Speech  understanding  problem  needs  multilevel  in- 
formation transformations  from  electrical  voice  signal  to 
recognized  sentence.  The  levels  include  parameter,  segment, 
syllable,  word,  word-sequence,  phrase,  and  data  base  inter- 
face levels.  At  each  level  different  knowledge  sources  are 
distributed  and  utilized  to  induce  higher  level  knowledge 
sources  or  hypotheses.  The  knowledge  sources  perform  a 
variety  of  functions,  such  as  extracting  acoustic  parame- 
ters, classifying  acoustic  segments  into  phonetic  classes, 
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recognizing  words,  parsing  phrase,  and  generating  and 
evaluating  predictions  for  undetected  words  or  syllables. 
Production  rules  are  used  to  represent  each  level  of 
knowledge . 

HEARSAY-II  attempts  to  reconstruct  a speaker's  inten- 
tion from  hypothetical  interpretations  obtained  at  various 
levels  of  abstraction.  The  goal  is  to  construct  the  most 
credible  overall  interpretation.  The  basic  cycle  of  the 
operations  in  construction  are  hypothesis  generation,  hy- 
pothesis combination,  and  hypothesis  evaluation.  At  each 
step  in  the  construction,  each  knowledge  source  builds  a 
larger  interpretation,  adding  their  constraints  to  the  in- 
terpretation. These  constraints  are  used  to  resolve  uncer- 
tainty in  the  data  and  in  the  knowledge  sources  themselves. 
For  this  incremental  problem  solving  method,  the  independent 
knowledge  sources  need  to  communicate  with  each  other 
through  a ’’blackboard''  which  contains  currently  competing 
hypotheses  and  data  of  various  knowledge  sources.  From  the 
blackboard  each  knowledge  source  gets  information  about  the 
intermediate  state's  partial  interpretation  and  generates  a 
hypothesis  on  the  blackboard.  The  information  on  the  black- 
board is  also  divided  into  various  levels  depending  on  the 
level  of  the  knowledge  sources.  Whenever  a hypothesis  is 
recorded  on  the  blackboard,  a scheduler  with  data-directed 
control  regime  which  utilizes  opportunistic  problem  solving 
method  [Ni82]  makes  a schedule  for  the  next  knowledge  source 
activation.  To  avoid  possible  combinatorial  explosion  of 
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knowledge  source  invocation  the  computing  resources  are  al- 
located to  the  most  important  and  most  promising  action.  A 
heuristic  schedular  [Ha77]  calculates  a priority  for  each 
action  and  executes  the  waiting  action  with  the  highest 
priorities.  Calculation  of  priority  is  based  on  the  black- 
board information  and  the  effects  of  each  action  of  the  hy- 
pothesis. 

Based  on  the  HEARSAY-II  architecture,  HEARSAY-III 
[Ha83b]  was  developed  as  a domain-independent  framework.  The 
HEARSAY-III  structure  is  similar  to  that  of  HEARSAY-II,  con- 
sisting of  a relational  data  base  system  and  its  control 
mechanism.  It  also  utilizes  a blackboard  structure  for  sys- 
tem wide  knowledge  source  communication. 

We  have  seen  several  examples  of  expert  systems  and  ex- 
pert system  design  tools  (more  such  systems  are  described 
elsewhere  [Ba82a]  [Ba82b]).  Most  of  the  systems  discussed 
work  only  in  a narrow  domain,  and  their  performances  are  not 
so  successful  as  to  have  led  to  their  daily  use.  In  CASNET, 
for  example,  the  problem  domain  was  not  clearly  understood 
in  the  initial  design  stage.  The  problem  domain  of  glaucoma 
itself  is  overlapped  with  the  other  opthalmological  problems 
in  which  the  different  knowledge  representation  schemes  or 
the  different  control  mechanisms  are  required.  Also  the 
knowledge  base  about  the  other  opthalmological  problems  was 
not  fully  supported  in  the  original  causal  network.  In  MY- 
CIN  and  PIP  there  are  no  general  scoring  rules  which  assign 
a precise  certainty  factor  to  a specific  finding.  The 
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heuristic  rules  cannot  be  always  properly  represented  using 
certainty  factors  or  probabilities.  The  effects  of  a cer- 
tainty factor,  i.e.,  the  propagation  of  a certainty  factor 
in  a decision  tree,  were  also  not  covered  completely.  The 
linear  relationship  between  different  nodes  in  a decision 
tree  is  not  always  linear.  Therefore,  it  is  very  difficult 
to  change  the  knowledge  basis,  i.e.,  the  system  refinement 
which  is  essentially  required  in  the  development  process  is 
very  difficult.  Also,  in  HEARSAY-II  the  uncertainty  of  mul- 
tilevel knowledge  sources  causes  the  same  scoring  problem 
and  propagation  problem  of  a certainty  factor  of  heuristics. 
In  INTERNIST  and  HEARSAY-II  the  problem  domain  of  internal 
medicine  or  of  the  vocabulary  to  be  recognized  is  too  big  to 
handle  using  specific  methods  and  limited  computational 
power.  The  large  search  space  has  to  be  divided  into  sub- 
spaces of  the  high  levels  of  expertise.  There  is  too  much 
knowledge  to  be  encoded  in  the  knowledge  base  for  these  sys- 
tems. The  multi-symptoms  of  complex  diseases  could  not  be 
fully  covered  in  INTERNIST.  The  limitation  of  the  knowledge 
base  in  HEARSAY-II  permits  only  about  one  thousand  spoken 
words  to  be  recognized  with  a speaker's  trained  speech  data. 

However,  each  system  has  a different  architecture 
corresponding  to  their  domain  specific  problem.  That  is,  the 
required  expert  knowledge  and  the  way  the  human  experts  use 
their  knowledge  to  solve  problems  are  different  and  can  be 
structured  appropriately  depending  on  the  characteristics  of 
the  problem  domain.  Surely,  this  is  the  reason  why  expert 
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systems  have  different  architectures  and  the  key  point  that 
has  to  be  considered  in  any  expert  system  design. 

Sys^enig  Design 

In  this  section  the  general  problems  of  the  expert  sys- 
tem and  general  design  steps  are  discussed  to  characterize 
the  design  problem  of  an  EEG  expert  system. 

The  major  problems  of  expert  systems  are  [Da82] 

1)  To  separate  the  domain  dependent  knowledge  base  from  the 
inference  mechanism  used  for  reasoning  and  controlling  the 
system  (separation  between  knowledge  and  meta  knowledge  for 
knowledge  accumulation):  In  general,  the  separation  of  the 
inference  mechanism  from  a knowledge  base  is  not  always 
clear.  The  inference  mechanism  itself  as  a knowledge  base 
may  be  embedded  into  the  system  knowledge  base. 

2)  To  develop  logically  powerful  and  expressive  representa- 
tions for  describing  concepts,  facts,  and  procedures 
(knowledge  representation):  The  knowledge  representation 
techniques  are  chosen  based  on  the  usage  of  the  particular 
knowledge  in  the  system.  There  is  no  best  way  for  knowledge 
representation.  Depending  on  the  characteristics  of  the 
problem  domain  the  proper  knowledge  representation  tech- 
niques have  to  be  used  for  system  efficiency. 

3)  To  capture  the  expert  knowledge  about  specific  inferences 
or  decisions  in  the  form  of  modular  rules  without  conflic- 
tion  or  redundancies  among  the  rules  (knowledge  base  manage- 
ment): Modularity  of  the  knowledge  is  dependent  on  the 
structure  of  the  knowledge.  Modular  rules  are  not  always 
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possible. 

4)  To  develop  a general  inference  mechanism  to  control  rea- 
soning (deduction  mechanism  for  decision  making):  The  more 
the  framework  is  generalized,  the  more  difficult  it  is  to 
implement  the  domain  specific  expertise. 

5)  To  develop  methods  of  facilitating  user  interaction  with 
the  system  by  a way  of  normal  conversation  (natural  language 
interface):  Natural  language  processing  is  another  big  area 
in  AI.  A natural  language  like  user  interface  makes  the  sys- 
tem more  intelligent. 

As  previously  discussed,  the  problems  of  expert  systems 
result  in  a compromise  between  system  performance  and  the 
ease  of  the  computer  implementation,  or  between  the  system 
expertise  and  the  size  of  the  problem  domain.  That  is,  the 
more  the  system  is  generalized,  the  less  is  the  expertise. 

The  process  of  building  expert  systems  can  be  described 
by  the  following  steps: 

1)  Problem  identification:  Problem  characteristics  are 
determined,  such  as  area,  scope,  possible  problem  solving 
methodologies,  required  resources,  constraints,  and  specifi- 
cation of  goals. 

2)  Conceptualization:  A detailed  description  of  the  informa- 
tion flow  is  determined,  such  as  information  transformation 
in  each  intermediate  stage,  data  and  hypotheses  control 
mechanisms,  and  possible  implementation  methodologies. 

3)  Formalization:  The  knowledge  representation  techniques, 
system  rule  language,  and  inference  mechanism  are  selected 
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appropriately . 

4)  Implementation:  The  rules  that  embody  the  knowledge  are 
implemented  on  a computer. 

5)  System  refinement:  The  system  is  evaluated  using  field 
data  and  modified  according  to  the  rule  modifications. 

Ej&  &£st.eiS  _D.gs_i.gn  .Cgncej}^ 

In  this  section  the  basic  concepts  of  the  EEG  expert 
system  design  are  described  based  on  the  design  steps 
described  in  the  previous  section.  A detailed  description 
of  the  system  design  is  given  in  the  following  chapters. 

1)  Problem  Identification:  The  system  will  interpret  mul- 
tichannel EEG  (three  EEG  and  one  EOG  channels)  for  sleep 
analysis.  Epilepsy  waveform  detection  and  sleep  disorder 
diagnosis  will  be  included  in  the  future.  The  knowledge 
base  of  the  system  consists  of  specified  waveform  detection 
techniques,  EEG  signal  understanding  rules,  and  clinical 
knowledge  about  the  corresponding  diagnosis. 

2)  Conceptualization:  The  system  has  three  different 

knowledge  levels:  the  waveform  detector  outputs  of  the  EEG 

signals,  the  channel  activity  description,  and  the  synthetic 
report  level.  The  specific  waveforms  are  detected  from  the 
EEG  signals  using  digital  waveform  detectors.  The  channel 
activities  are  described  based  on  the  waveform  detector  out- 
puts. A synthetic  report  is  generated  using  the  channel 
description  list  and  corresponding  clinical  knowledge  bases. 

3)  Formalization:  At  the  waveform  level  each  waveform  is  ex- 
pressed as  a hexadecimal  character.  All  channel  activities 
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at  an  instant  form  a word.  Using  these  words  the  channel 
activities  are  described  in  the  channel  description  level. 
The  channel  description  list  is  utilized  for  synthetic  re- 
port generation.  To  realize  these  processes  in  the  main 
processor,  the  language  understanding  mechanism  is  utilized. 
The  EEG  signal  language  is  defined  and  the  EEG  signal  com- 
piler is  designed  using  existing  software  of  the  UNIX 
operating  system.  To  display  the  line  of  reasoning,  the 
channel  description  list  is  saved  in  a system  utility  file. 
On  the  synthetic  report  level,  interpretive  knowledge  such 
as  EEG  rule  understanding  is  incorporated  with  the  channel 
description  list  to  generate  a synthetic  report.  The  pro- 
duction system  and  frame  scheme  are  used  to  represent  the 
EEG  signal  understanding  rules  and  the  clinical  knowledge. 

4)  Implementation:  One  signal  processor  is  used  as  a fron- 
tend  processor  to  generate  "tokens".  The  frontend  processor 
is  implemented  using  a TI-990/101MA  board  level  microcomput- 
er and  an  A/D  board  (RTI-1241,  Analog  Device).  A PDP-11/23+ 
microcomputer  will  be  used  as  a main  frame.  The  token  gen- 
erator is  written  in  assembly  language  after  considering  the 
data  processing  speed.  For  the  token  processor  the  channel 
description  module  is  written  in  "C"  language.  This  module 
is  automatically  generated  using  a compiler  generator 
mechanism.  The  synthetic  reporter  is  written  in  LISP  for 
processing  a channel  description  list. 

5)  System  Refinements:  The  system  performance  is  tested  us- 
ing previously  recorded  EEG  data  tapes.  Once  the  prototype 
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is  built  the  system  can  be  tested  against  human  experts. 
System  refinements  such  as  rule  modifications  can  be  contin- 
ued up  to  the  level  of  a human  expert.  Various  user  friend- 
ly modules  can  be  added  at  this  stage. 

Although  these  are  the  expert  system  design  processes 
for  the  target  system,  in  the  present  system  not  all  the 
features  are  implemented.  In  the  next  two  chapters  the 
design  details  of  the  EEG  expert  system  are  discussed  and 
the  system  status  is  described  precisely. 


CHAPTER  4 

DESIGN  AND  IMPLEMENTATION  - I 
This  chapter  describes  the  human  EEGer's  sleep  data 
recognition  procedures  in  detail  to  provide  an  appropriate 
EEG  expert  system  model.  After  discussing  the  system  infor- 
mation flow  and  structure,  each  system  block  design  and  im- 
plementation of  the  frontend  processor  are  described  in  de- 
tail. The  main  processor,  the  system  status  and  future  sys- 
tem expansion  are  discussed  in  the  next  chapter. 

Design 

The  human  EEGer  analyzes  the  sleep  EEG  data  by  looking 
at  an  epoch  (one  minute  or  thirty  seconds  data)  of  the  data. 
In  epoch,  the  specified  waveforms  are  recognized  by  consid- 
ering all  the  channel  activities.  Sleep  staging  is  per- 
formed by  looking  at  all  the  activities  in  an  epoch.  If  am- 
biguity occurs  within  an  epoch  of  data  a couple  of  the  pre- 
vious or  following  epochs  are  checked  to  see  the  change  of 
waveform  activities,  such  as  the  amplitude  variation  or  wave 
rhythm  change  between  consecutive  epochs.  In  the  epoch 
analysis  specified  waveforms,  such  as  alpha,  beta,  delta, 
spindle,  etc.,  are  recognized  as  a group  of  waveforms  by  the 
number  of  local  peaks  and  their  duration.  For  instance,  to 
detect  an  alpha  waveform  (8  to  12  Hz),  the  number  of  the 
waves  is  counted  in  a certain  interval,  which  may  be  a one 
or  two  second  interval.  If  at  least  six  consecutive  waves, 
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whose  duration  lie  in  the  range  of  the  alpha  wave  are  found, 
then  they  are  regarded  as  alpha  waveforms.  Of  course,  these 
numbers  do  not  precisely  represent  the  visual  analysis  cri- 
teria. Also,  the  selected  alpha  waveform  itself  cannot  be 
clearly  defined.  They  can  be  consecutive  waveforms  or  a 
group  of  split  waveforms  [Sm78].  In  some  aspects,  it  is  sub- 
jective whether  a certain  waveform  is  regarded  as  an  alpha 
waveform  or  not.  However,  by  adding  up  the  alpha  waveforms 
in  an  epoch  the  total  alpha  time  of  an  epoch  is  obtained. 
In  the  adding  procedure,  because  of  the  ambiguous  waveforms 
the  human  EEGer  cannot  make  a clear  distinction  of  the 
waveforms  and  the  exact  amount  of  the  total  alpha  activity 
time  cannot  be  obtained.  Usually  an  approximation  of  each 
second  is  made  and  added  up  to  get  the  total  alpha  time  in 
an  epoch.  The  sleep  staging  rules  are  applied  on  this  total 
alpha  time. 

On  the  other  hand,  in  the  case  of  the  delta  wave  (0.5 
to  2.0  Hz),  usually  the  number  of  peaks  in  a certain  inter- 
val are  counted  or  the  period  between  the  local  peak  points 
is  compared  to  the  interval  of  0.4  seconds,  denoted  by  lines 
on  the  data  sheet.  In  this  case  the  local  peaks  of  the  high 
frequency  components  or  the  fast  activities  are  just  ignored 
or  smoothed  to  count  only  the  slow  wave  components.  The  to- 
tal delta  activity  time  is  obtained  in  much  the  same  way  as 
the  alpha  time. 

While  the  alpha  wave  and  delta  wave  recognition  pro- 
cedures require  the  checking  of  only  one  channel,  the  REM 
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waveform  recognition  procedure  needs  to  check  all  the  chan- 
nel activities.  When  a slow  wave  appears  in  the  EOG  chan- 
nel, the  REM  waveform  recognition  is  totally  dependent  on 
the  other  channel  activities.  The  slow  wave  activities  of 
the  frontal  channel  usually  affect  the  EOG  channel  waveform. 
When  a K-complex  waveform  occurs  it  appears  on  more  than  two 
channels  simultaneously.  The  EOG  channel  can  be  one  of 
them.  When  delta  activities  appear  on  the  occipital  channel, 
the  EOG  channel  also  usually  shows  slow  wave  activities. 
Thus  the  EOG  channel  slow  wave  has  to  be  interpreted  by  con- 
sidering all  the  channel  activities.  On  the  other  hand  in 
the  REM  period  detection,  the  epoch  which  is  between  REM 
periods  is  regarded  as  a REM  period  without  considering  the 
slow  wave  existence  on  EOG  channel.  Therefore,  for  sleep 
data  interpretation  all  the  channel  activities  have  to  be 
checked  and  in  addition  the  preceding  and  the  following 
epochs  have  to  be  considered  simultaneously. 

Then  how  can  these  processes  be  realized  on  the  comput- 
er? As  previously  shown,  the  procedures  or  steps  of  the  hu- 
man EEGer's  sleep  EEG  data  analysis  require  individual 
waveform  recognition  and  contextual  analysis  of  these 
waveforms.  Especially  for  the  contextual  analysis  of  the 
waveform  the  descriptive  knowledge  or  procedural  knowledge 
has  to  be  utilized  with  the  channel  information.  To  ap- 
proach these  kinds  of  waveform  recognition  and  interpreta- 
tion problems,  an  expert  system  is  designed.  Each  step  of 
the  sleep  EEG  analysis  is  simulated  using  a hybrid  system. 
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A signal  processor  which  consists  of  an  A/D  board  and  a 
microcomputer  detects  sleep  waveforms.  A main  processor 
which  is  a software  program  of  a minicomputer  performs  con- 
textual analysis  of  the  detected  waveforms.  Though  this 
system  does  not  exactly  simulate  the  EEG  data  processing 
procedures  of  the  human  EEGer,  in  every  step  of  the  recogni- 
tion procedure  an  attempt  is  made  to  incorporate  the 
corresponding  techniques  of  waveform  detection  or  knowledge 
base  of  waveform  interpretation  to  the  system  design.  Each 
specified  waveform  is  detected  using  a waveform  detector 
which  simulates  the  detection  process  of  the  human  EEGer  and 
then  the  waveforms  are  interpreted  utilizing  a contextual 
analysis  which  utilizes  all  the  channel  information  and  the 
subject  dependent  data.  Although  there  are  minor  changes  in 
the  waveform  detection  criteria,  the  waveform  detection 
techniques  are  mainly  based  on  the  techniques  of  the  SAHC 
[Sm78]  [Pr79l  which  performed  well  in  the  sleep  EEG  waveform 
analysis. 

In  the  present  system  all  the  analog  filters  and 
hardware  detectors  of  the  SAHC  are  converted  to  software 
programs  (assembly  language  on  a TI-990/101MA  board  level 
microcomputer).  Slow  wave  detectors  are  added  to  each  chan- 
nel for  the  REM  period  detection.  For  the  contextual 
analysis  of  the  sleep  EEG  signal  the  signal  language  is  de- 
fined and  semantic  analysis  is  performed.  The  computer  com- 
piler mechanism  is  utilized  for  signal  language  understand- 
ing and  the  standard  software  (LEX,  YACC)  for  the  compiler 
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generation  is  adapted  for  signal  compiler  generation.  The 
detected  waveform  patterns  and  the  relationships  between 
patterns  are  matched  to  the  tokens  and  the  syntactic  struc- 
ture of  the  compiler,  respectively.  The  analysis  of  these 
syntactic  patterns  in  the  detected  waveforms  corresponds  to 
the  semantic  action  of  the  compiler.  Furthermore,  the  au- 
tomated procedure  of  compiler  generation  gives  great  flexi- 
bility in  rule  management.  The  waveform  interpretation 
rules  in  signal  grammar  and  semantic  action  file  can  be 
changed  independently  of  the  other  part  of  the  system.  Au- 
tomatic system  generation  is  attempted  to  simplify  the  fre- 
quent rule  changes  which  essentially  occur  in  the  expert 
system  development  processes. 

Although  the  target  system  is  to  be  implemented  on  a 
PDP-11/23+  microcomputer  with  UNIX  operating  system,  the 
current  system  is  developed  on  a GOULD  Concept  32/8780  mini- 
computer using  the  UNIX  operating  system.  The  portability 
of  the  UNIX  operating  system  and  the  use  of  the  standard 
software  and  language  "C"  make  system  conversion  a simple 
process.  For  the  data  pipeline  in  the  system  the  system 
shell  command  is  developed  using  the  UNIX  system  shell  com- 
mand. The  language  "LISP"  is  selected  for  the  processing  of 
the  channel  description  list  and  the  system  knowledge  base 
for  the  sleep  disorder  diagnosis.  LISP  will  enhance  the  ca- 
pability of  the  interactive  user  interface  and  system  rea- 
soning. 
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Fig.  4.1  shows  the  system  information  flow.  The  key 
concept  of  designing  the  system  is  to  simulate  the  interpre- 
tation procedures  of  the  electroencephalographer  in  analyz- 
ing sleep  data.  From  the  multichannel  EEG/EOG  inputs,  a 
synthetic  report  of  the  sleep  data  analysis  is  generated. 
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Fig.  4.1  Block  Diagram  of  System  Information  Flow. 


The  system  is  characterized  as  a distributed  multilevel 
(or  cooperative)  knowledge  base  expert  system.  In  each  step 
of  the  data  understanding  or  the  interpretation  procedure 
the  various  knowledge  sources  interact  with  each  action 
module,  such  as  the  filter  and  detector  bank,  channel 
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descriptor  and  synthetic  reporter.  Each  action  module  ac- 
cepts the  data  from  the  previous  module  and  processes  the 
data  utilizing  its  knowledge  source.  In  this  process  the 
data  flow  occurs  only  in  one  direction  between  the  modules, 
from  low  level  to  high  level. 

The  knowledge  base  is  distributed  throughout  the  system 
at  each  information  processing  level  such  as  parameter  lev- 
el, channel  level  and  synthetic  reporter  level.  In  the 
parameter  level  the  information  about  the  data  sampling  such 
as  channel  assignment,  sampling  rate  and  channel  multiplex- 
ing is  incorporated  in  a data  handling  module  (interrupt 
handler  and  main  program  of  the  token  processor  of  the 
present  system,  refer  to  Appendix  2).  In  the  filter  bank 
(Fig.  4.2)  the  required  cascaded  filter  information  such  as 
amplitude  characteristics,  phase  response,  and  transient 
response  are  incorporated  in  each  filter  module.  All  the 
detector  criteria  for  the  specified  waveforms  such  as  the 
required  window  size  of  the  specified  waveform  detection, 
wave  periods,  number  of  consecutive  waves,  amplitude  thres- 
hold value,  hysteresis,  minimum  display  time,  wave  period 
averaging  process,  etc.  are  incorporated  in  each  detector 
module.  In  the  channel  level  the  channel  description  rules 
for  the  specified  waveforms,  such  as  waveform  significance 
resolution,  waveform  relation  resolution  in  time  domain  and 
waveform  smoothing  for  total  activity  time  calculation  are 
incorporated  in  the  channel  descriptor.  The  channel 
description  rules  are  in  the  rule  files  (lex  source  files, 
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grammar  files  and  semantic  action  files  of  the  present  sys- 
tem) of  the  channel  descriptor  module.  In  the  synthetic  lev- 
el an  epoch  description  rules  such  as  REM  period  recognition 
and  sleep  staging  rules  are  incorporated  in  the  synthetic 
reporter  module. 

The  system  consists  of  two  functional  blocks:  the  fron- 
tend  processor  and  the  main  processor.  The  frontend  proces- 
sor (which  is  also  called  a token  generator),  as  a signal 
processor,  includes  an  A/D  converter,  a filter  and  a detec- 
tor bank.  The  main  processor  (which  is  also  called  a token 
processor),  as  an  interpreter,  includes  an  interpreter, 
knowledge  sources,  system  utilities,  and  controller.  The 
frontend  processor  detects  the  specific  waveforms  and  gen- 
erates the  token  stream.  In  the  detector  bank  all  the 
waveform  detector  criteria  are  embedded  into  the  different 
detector  modules.  The  main  processor  interprets  the  token 
stream  and  generates  a written  report.  The  knowledge  base  of 
the  main  processor  consists  of  channel  description  rules  and 
REM  detection  rules  for  sleep  staging.  The  knowledge  base 
will  be  expanded  to  include  the  sleep  disorder  diagnosing 
rules. 

The  signal  flow  of  the  present  system  is  as  follows. 
The  analog  input  signals  from  the  EEG/EOG  amplifiers  are  di- 
gitized by  an  A/D  converter  and  passed  through  the  digital 
filter  section.  The  cascaded  filters  (Fig.  4.2)  with  dif- 
ferent sampling  rates  [Pr79]  approximate  the  first  step  of 
waveform  recognition  of  the  human  EEGer  by  selecting  the 
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Fig.  4.2  The  Frontend  Structure  of  the  System 
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specified  frequency  component  from  the  EEG/EOG  signals.  The 
corresponding  detectors  then  detect  the  specified  waveforms 
appearing  in  the  corresponding  channel  utilizing  waveform 
criteria  [Sm78],  All  the  detector  outputs  are  sampled 
simultaneously  and  coded  into  a two  byte  word.  This  two  byte 
word  forms  the  input  token  to  the  channel  descriptor  (Fig. 
4.3).  This  process  in  which  all  the  channel  information  of 
the  same  instant  are  multiplexed,  provides  the  way  of  infor- 
mation checking  across  the  channels  at  an  instant.  These  to- 
kens are  then  interpreted  by  the  interpreter  of  the  main 
processor  using  heuristic  rules  derived  to  mimic  the  human 
expert. 

The  interpreter  consists  of  a channel  descriptor  and  a 
synthetic  reporter.  Each  of  these  modules  utilizes  the 
corresponding  knowledge  sources  to  process  the  token  stream. 
In  the  channel  descriptor  each  channel  activity  is 
described,  based  on  the  total  activity  time  of  the  specified 
waveform  in  each  channel.  The  token  stream  which  includes 
all  the  detected  waveform  information  is  regarded  as  a sig- 
nal sentence  and  interpreted  to  generate  a synthetic  report 
in  each  epoch.  The  descriptive  words  for  the  channel  ac- 
tivities are  previously  selected  depending  on  the  waveform 
characteristics  and  the  required  analysis. 

Fig.  4.3  shows  the  interpreter  structure  of  the  main 
processor  in  the  system.  In  the  channel  decoder  the  lexical 
analysis  of  the  proper  tokens  is  performed  based  on  the 
waveform  or  channel  priority,  namely  selecting  a correct 
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detector  output.  The  channel  priority  and  waveform  priority 
are  resolved  in  the  decoder.  The  syntactic  and  semantic 
analysis  of  output  relationships  between  different  tokens  of 
the  same  channel  is  performed  on  the  output  of  the  decoder. 
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Fig.  4.3  Block  Diagram  of  Main  Processor 
( Token  Processor  ) 

Then  the  translator  generates  the  description  list  of  all 
the  waveform  activities  of  the  corresponding  channel. 
Heuristic  rules  for  interpreting  these  relations  are  applied 
to  this  stage  of  analysis.  The  channel  description  list  is 
used  for  the  synthetic  reporter. 

The  synthetic  reporter  makes  the  final  synthetic  report 
after  considering  all  the  previous  channel  descriptions. 
Many  of  the  technical  rules  for  interpreting  the  EEG  data 
are  implemented  in  the  synthetic  reporter  based  on  the  hy- 
brid production  system  and  frame  scheme  [ Wi7 9 3 . The  rules 
for  the  synthetic  channel  description  are  implemented  using 
a production  system  [Ni80],  The  structure  of  the  channel 
description  rules  in  the  production  system  simplifies  the 
connection  between  the  EEG  data  and  the  other  knowledge 
base,  such  as  sleep  disorder  diagnosing  rules.  Both  the 
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channel  description  rules  and  sleep  disorder  diagnosing 
rules  will  be  put  into  a rule  file.  This  will  unify  the  sys- 
tem rule  management.  The  frame  scheme  [ Mi7 6 3 is  considered 
for  the  sleep  disorder  diagnosing  rules  of  EEG  diagnosis 
(not  used  in  the  present  system,  where  only  the  basic  infer- 
ence mechanism  is  implemented).  A clear  distinction  between 
the  knowledge  base  for  the  channel  descriptor  and  the 
knowledge  base  for  the  synthetic  reporter  cannot  be  made 
though  they  have  different  mechanical  structures  for  infor- 
mation processing.  Some  of  the  rules  can  be  implemented  in 
either  place  without  affecting  the  final  decision.  The  de- 
tailed discussion  of  each  block  design  and  implementation 
are  given  in  the  following  section. 

The  frontend  processor  generates  tokens  every  quarter 
second  utilizing  an  A/D  converter,  a filter  bank,  a detector 
bank,  and  a multiplexor.  The  frontend  processor  except  for 
the  A/D  converter  is  implemented  on  TI-990/101MA  board  level 
microcomputer.  The  architecture  of  the  independent  work 
space  pointer  of  the  TI-990/101MA  microcomputer  causes  in- 
dependent operation  of  each  filter  and  detector  and  speeds 
up  the  linking  operation  of  each  module.  The  expansion 
memory  board  of  32K  bytes  memory  is  used  for  the  program 
space  (4K  bytes)  and  the  generated  token  storage  space  ( 28 K 
bytes).  The  on-board  memory  of  4k  bytes  is  also  used  for 
the  program  space.  The  functional  flowchart  of  the  frontend 
processor  is  shown  in  Fig.  4.4.  The  program  consists  of  a 
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main  program,  12  filter  modules,  10  detectors,  and  an  inter- 
rupt handler. 
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Fig.  4.4  Functional  Flowchart  of  the  Frontend. 


In  the  beginning  of  the  main  program  parameter  initializa- 
tion is  performed.  After  setting  up  the  parameters,  each 
filter,  detector,  and  interrupt  handler  are  called  in  the 
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main  program  according  to  the  previously  selected  order. 
Once  the  program  starts,  the  third  level  interrupt  is  gen- 
erated every  1/480  second.  In  that  interval  the  related 
filters  and  detectors  detect  the  specified  waveforms.  Then 
the  multiplexor  generates  tokens  every  quarter  second,  using 
detected  waveforms.  Each  functional  block  is  discussed  in 
detail  in  the  following  sections. 

In  the  initial  setup  procedure  various  addresses  are 
defined,  such  as  that  of  each  filter  and  detector  module, 
I/O  port  of  each  module,  interrupt  handler,  A/D  converter 
base  address,  starting  address  of  the  processed  data 
storage,  and  various  clocks  and  counters.  All  the  variables 
are  also  initialized,  such  as  sampling  rate,  control  parame- 
ters of  the  A/D  converter,  and  various  clocks  and  counters 
(see  main  program  in  Appendix  1 ) . To  control  the  A/D  con- 
verter, first  the  CRU  (communication  register  unit)  base  ad- 
dress is  set  to  A/D  converter  base  address  and  various  con- 
trol registers  of  the  A/D  board  are  set.  The  CRU  of  the  CPU 
board  is  used  to  control  the  mapped  I/O  of  the  system.  Us- 
ing a register  in  the  main  program  the  auto  scanning  mode  of 
the  multiplexor  is  disabled  to  check  only  four  input  chan- 
nels. The  sampling  port  address  is  handled  in  the  interrupt 
handler  by  one  of  the  address  pointers  in  the  interrupt 
handler.  The  timer  is  initialized  to  generate  interrupt  at 
480Hz.  The  gain  of  the  sampler  is  set  to  one  in  this  step. 
For  the  proper  operation  of  the  interrupt  handler  the  jump 
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address  of  the  third  level  interrupt  is  set  to  a prefixed 
address  of  the  system.  The  starting  address  of  the  pro- 
cessed data  is  loaded  by  setting  one  of  the  interrupt 
handler  registers.  That  register  is  used  as  an  address 
pointer  in  autoincrement  address  mode.  Most  of  the  registers 
in  the  main  program  are  used  for  the  subroutine  pointers, 
making  the  jump  to  the  subroutine  in  the  main  program  a 
quick  one.  After  setting  up  all  the  variables  the  interrupt 
masks  of  the  PSI  (programmable  system  interface)  and  the 
frontend  CPU  are  enabled.  The  sampling  operation  starts 
after  getting  a third  level  interrupt  signal  from  the  timer 
on  the  A/D  board.  Then  A/D  conversion  is  initiated  by  writ- 
ing a dummy  number  on  the  conversion  start  register  of  the 
A/D  converter  board.  Once  the  sampling  procedure  starts 
then  it  repeats  sampling  operation  at  480  Hz. 
i/D  Conveptep 

The  signals  from  three  EEG  channels  and  one  EOG  channel 
are  converted  into  digital  signals  using  an  A/D  converter 
(RTI-1241  Analog  Device).  A 12  bit  A/D  converter  was 
selected  to  provide  an  adequate  signal  to  noise  ratio 
[PRI79].  The  base  address  of  the  A/D  board  is  selected  us- 
ing jumpers  present  on  the  board.  The  A/D  board  operation 
is  controlled  by  the  initial  setup  procedure  and  by  the  in- 
terrupt handler.  In  the  initial  setup  procedure,  all  the 
parameters  are  initialized  for  the  A/D  board  to  operate  in 
the  multiplexed  normal  sampling  mode  after  setting  up  the 
CRU  of  the  CPU  board.  For  the  four  channel  multiplex 
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sampling  the  auto  scanning  of  the  input  channel  is  disabled. 
The  input  channel  addresses  are  written  in  the  multiplexor 
address  register  every  four  channel  sampling  cycle  of  the 
interrupt  handler.  A register  of  the  interrupt  handler  is 
used  for  the  data  input  channel  control.  This  register  is 
used  as  an  address  pointer  in  the  autoincrement  address 
mode.  The  polled  status  mode  of  the  normal  operation  modes 
which  checks  the  end  of  conversion  bit  to  determine  the  com- 
pletion of  the  A/D  conversion,  is  used  to  maximize  the  flex- 
ibility and  the  frontend  CPU  efficiency.  In  the  normal  mode 
operation,  once  the  multiplexor  address  and  the  gain  are  es- 
tablished the  conversion  is  proceeded  by  writing  a conver- 
sion command.  Then  the  frontend  CPU  is  freed  to  perform  oth- 
er tasks. 

The  A/D  board  starts  the  operation  with  the  interrupt 
generated  by  the  timer.  The  third  level  interrupt  is  gen- 
erated every  sampling  interval  by  the  external  timer  of  the 
TMS-9901  programmable  system  interface  of  the  A/D  board. 
Then  the  interrupt  handler  resets  the  external  timer  and 
makes  the  A/D  converter  collect  samples  at  every  channel, 
and  then  transfers  the  control  over  the  filter  module. 

The  interrupt  handler  as  an  interrupt  service  routine, 
is  called  every  sampling  interval.  First,  the  interrupt 
handler  controls  the  data  sampling  and  stores  the  generated 
tokens  in  the  memory.  Four  input  data  channels  are  sampled 
in  a prefixed  order  by  keeping  the  input  channel  addresses 
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and  are  sent  to  the  corresponding  filter  input.  The  first 
filter  bank  consisting  of  four  120  Hz  low  pass  filters  in 
each  channel,  is  implemented  in  the  interrupt  handler  to 
reduce  the  subroutine  switching  time  of  the  cascaded  low 
pass  filters.  These  low  pass  filter  outputs  are  supplied  to 
each  waveform  filter  with  different  sampling  frequency.  In 
the  interrupt  handler  several  timers,  clocks,  and  address 
pointers  are  kept  for  token  generation  and  acquisition.  The 
detector  outputs  are  multiplexed  into  a two  byte  word  token 
every  quarter  second  and  stored  in  the  extended  memory 
boards.  The  control  signal  for  the  data  display  on  paper  is 
generated  every  two  seconds  with  a half  second  turn-on  time. 
The  statistics  of  the  processed  data,  such  as  total  alpha 
time,  beta  time,  number  of  spindles,  etc.  are  also  stored  in 
memory  every  minute,  After  servicing  all  these  jobs,  the  in- 
terrupt handler  resets  the  interrupt  mask  which  enables  the 
third  level  interrupt  of  the  CPU  and  the  programmable  system 
interface . 

Filter;  Banjs 

The  filter  bank  consist  of  12  filters  (except  four  low 
pass  filters  in  the  interrupt  handler).  Alpha  and  delta 
filters  consist  of  three  different  low  pass  and  band  pass 
filters  in  cascade  (see  Fig.  4.2).  The  others  consists  of 
two  filters  in  cascade.  The  filtering  process  is  necessary 
in  order  for  the  zero  crossing  detector  to  identify  a high 
frequency  component  superimposed  on  slow  waves,  or  to  elim- 
inate the  effects  on  the  zero  crossings  of  high  frequency 
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components  superimposed  on  a slow  wave.  Actually,  in  visual 
analysis  of  the  EEG  data  the  human  EEGer  also  detects  a 
specified  waveform  having  these  kind  of  filtering  operations 
in  mind.  For  instance,  the  number  of  local  peaks  on  a slow 
wave  is  counted  to  select  a high  frequency  component  or  the 
local  peaks  of  the  high  frequency  component  on  slow  wave  are 
just  ignored  to  select  the  slow  wave. 

The  filter  corner  frequencies  in  each  channel  are 
selected  depending  on  the  periods  of  the  waves  to  be  select- 
ed. Linear  phase  finite  impulse  response  (FIR)  filters  were 
designed  using  algorithms  which  eliminate  multiplications  in 
the  filter  implementation  by  selecting  filter  coefficients 
that  are  the  sum  of  2**(-n ) [Ly77l  and  using  only  zeros  in 
the  filter  transfer  function  to  obtain  the  specified  fre- 
quency characteristic  of  the  digital  filter  C Pr 8 5 1 . The  out 
of  band  frequency  gain  is  suppressed  by  putting  zeros  close- 
ly on  the  unit  circle  of  the  Z plane.  This  algorithm  made 
the  real  time  operation  of  the  filter  bank  possible.  The 
digital  filters  are  realized  using  the  work  space  registers 
(up  to  12  registers).  The  registers  are  used  as  delay  units 
with  fixed  periods.  The  delay  time  corresponds  to  each 
filter’s  sampling  rate.  The  coefficient  of  the  filter  is 
realized  by  ADD  and  SHIFT  instruction.  For  instance,  1 .25  * 
A is  obtained  by  the  operation  of  (A  + (SHIFTRIGHT 
(SHIFTRIGHT  A))).  Filter  phase  linearity,  which  can  minim- 
ize the  distortion  of  zero  crossing,  is  obtained  by  posi- 
tioning the  zeros  at  the  specified  positions  on  the  unit 
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circle  of  the  Z-plane.  In  case  of  narrow  band  low  frequency 
filters,  the  undersampl ing  effect  is  minimized  by  utilizing 
a maximum  of  three  cascaded  filters,  which  have  different 
sampling  rates.  The  first  filter  of  the  cascaded  filters  has 
a higher  sampling  rate  and  a wider  bandwidth  than  the  oth- 
ers. It  makes  the  first  folding  frequency  range  wide.  The 
second  filter  has  a bandwidth  which  is  just  larger  than  the 
half  of  the  first  one  (if  there  is  a second  filter  in  the 
cascaded  filter).  Then  the  folding  frequency  bandwidth  can 
be  larger  than  the  case  of  one  filter,  even  with  a low  sam- 
pling frequency.  Also,  using  this  structure  the  low  fre- 
quency filters  with  narrow  bandwidth  are  designed  with  a 
smaller  number  of  transmission  zeros.  This  means  simpler 
structures  and  faster  filter  operation  for  the  real  time 
data  processing.  Usually,  the  narrow  pass  band  filter  im- 
plementation with  high  sampling  frequency  requires  many 
zeros  on  the  unit  circle  to  suppress  the  gain  of  out  band 
frequency . 

Though  the  filters  are  satisfactory  for  practical  ap- 
plications, the  complete  elimination  of  frequency  folding 
from  high  frequency  to  low  frequency  due  to  aliasing  cannot 
be  achieved  (see  the  alpha  filter  frequency  response  in  Ap- 
pendix 1).  This  is  mainly  due  to  the  maximum  possible  sam- 
pling frequency,  which  is  determined  by  the  speed  of  fron- 
tend  processor.  Although,  using  this  method  the  linear 
phase  shift  can  be  obtained,  the  linear  phase  intercept  dis- 
tortion [Pe77]  is  not  always  achieved.  The  DINAP  [Ba78] 
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program  was  utilized  to  check  the  frequency  response  and 
phase  characteristic  of  each  filter.  (See  Appendix  2 for  the 
filter  characteristics  and  implementation). 

P£tec.fcO£  S^nk 

The  detector  algorithm  (Fig.  4.5)  used  for  this  system 
is  based  on  an  analog  EEG  waveform  detector  which  performs 
well  in  sleep  analysis  [Sm753.  There  was  a minor  change  in 
step  (5)  (see  Table  4.1).  Step  (6)  is  added  for  the  token 
generation  which  samples  detector  outputs  every  quarter 
second.  Once  a certain  waveform  is  detected  then  the  detec- 
tor output  should  be  displayed  to  be  sampled  correctly.  Six 
detectors  for  alpha,  beta,  delta,  sigma,  theta  and  artifact 
are  assigned  (as  shown  in  Fig.  4.2). 

(1)  Check  zero  crossing 

(2)  If  necessary,  check  minimum  amplitude  of  wave 

(3)  Check  the  period  of  a half  wave  for  slow  waves  and 
a full  wave  for  the  others 

(4)  If  necessary,  use  the  wave  averaging  method  for 
the  higher  resolution 

(5)  Apply  the  wave  criteria  to  each  wave 

(6)  Display  the  detector  output  for  each  minimum  period 

Fig.  4.5  Waveform  Detector  Algorithm 

A slow  wave  detector  in  each  EEG  channel  plus  a slow  wave 
detector  with  different  bandwidth  in  the  EOG  channel  are 
used  for  REM  detection.  The  waveform  detector  first  checks 
the  zero  crossing  of  each  waveform.  The  sign  of  the  sampled 
data  value  is  checked  for  the  zero  crossing  detector. 
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Table  4.1 
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Wavef  orm : 
Sr  : 
A : 
Period  : 
Wmax  : 

Win  : 
Wout  : 
Dm  in  : 

Average  : 

Pr  : 
* 


Sleep  waveforms  and  slow  waveforms 
Sampling  rate  (Hz) 

Detector  amplitude  threshold  value  (uV) 
Wave  period  (number  of  sampling  intervals) 
Maximum  number  of  waves  kept  in  the 
detector  window 

Number  of  waves  for  detector  on 
Number  of  waves  for  detector  off 
Minimum  detector  display  time  (number  of 
sampling  interval) 

Average  wave  period  (number  of  sampling 
interval ) 

Previous  wave  period 
No  criteria 


Then  the  maximum  amplitude  of  each  waveform  is  tested.  If 
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the  amplitude  is  too  small,  the  wave  is  ignored.  After 
checking  the  amplitude,  the  period  criteria  are  applied  for 
zero  crossing  detection.  For  low  frequency  waveform  detec- 
tion, i.e.,  delta  (0. 5-2.0  HZ)  or  slow  wave(0.5  - 6.0  HZ), 
half  cycle  detection  is  used  instead  of  full  cycle  detec- 
tion. The  specified  waveforms  are  detected  based  on  the 
number  of  wave  criteria  [Sm783.  For  zero  crossing  detec- 
tion, the  minimum  sampling  rate  fs  for  a frequency  resolu- 
tion fr  [Sm79l  [Sm80]  is  given  by 

fr=fi*fi/fs,  where  fi  is  the  input  frequency. 
Although  this  condition  cannot  be  satisfied  because  of  the 
limit  on  the  maximum  sampling  rates,  an  averaging  method  is 
utilized  to  increase  the  frequency  resolution  of  the  high 
frequency  waveforms.  In  the  averaging  method  the  specified 
number  of  the  consecutive  zero  crossing  intervals  are  summed 
together.  If  the  consecutive  waveforms  satisfy  the  wave  cri- 
teria then  the  total  zero  crossing  intervals  are  averaged 
and  compared  to  the  standard  wave  period  of  the  correspond- 
ing waveform. 

For  REM  detection,  three  slow  wave  detectors  in  three 
different  EEG  channels  check  the  existence  of  slow  waves  in 
each  channel.  If  a slow  wave  exists  in  any  EEG  channel,  the 
REM  like  slow  wave  of  the  EOG  channel  is  reexamined  in  a 
subsequent  context  analysis.  If  a delta  wave  or  a slow  wave 
exists  in  the  occipital  or  central  channel  then  the  slow 
wave  in  EOG  channel  is  not  regarded  as  a REM.  Also,  if  a 
delta  wave  exists  then  a slow  wave  in  the  previous  and  the 
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following  2 seconds  is  not  regarded  as  a REM.  Once  the 
specified  waveform  is  detected,  the  detector  output  is 
turned  on  for  a minimum  period  of  each  waveform,  providing 
for  the  token  generation.  Because  the  token  is  generated 
every  quarter  second,  the  waveform  detector  output  should  be 
longer  than  this  period  to  be  sampled  correctly.  This 
minimum  display  time  differs,  depending  on  the  waveforms. 
For  instance,  the  delta  detector  displays  the  previous  delta 
wave  period  and  the  spindle  detector  displays  7/16  second 
(see  each  detector  program  in  Appendix  1). 

.Bs-tecj^-E  SmeXi-Rg 

The  waveform  detector  outputs  are  sampled  every  quarter 
second  and  stored  in  a two  byte  word.  Two  bytes  are  used 
for  each  sample.  Ten  bits  are  used  for  the  detector  output, 
one  bit  is  used  for  time  marking,  and  five  bits  are  avail- 
able for  future  expansion.  The  relative  bit  positions  of 
each  waveform  detector  output  in  a two  byte  word  are 
represented  as  follows 

[ , , ,time  mark]  [alpha,  beta,  delta,  sigma] 

[theta,  artifact,  , ] [rslow,  oslow,  pslow,  qslow] . 

Each  bracket  denotes  a hexadecimal  number  of  a four  bit 
binary  notation.  For  instance,  the  alpha  waveform  detector 
output  is  denoted  in  the  first  bit  of  the  second  nibble. 
The  simultaneous  sampling  of  the  detector  outputs  makes  it 
possible  to  analyze  waveform  relationships  between  the  dif- 
ferent channels  or  of  the  same  channel.  The  waveform  detec- 
tor outputs  form  the  input  token  to  the  main  processor.  With 
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this  sampling  scheme  the  waveform  priority  is  given  easily 
in  a subsequent  main  processor  context  analysis.  This  sam- 
pled output  of  a two  byte  word  is  converted  to  four  hexade- 
cimal ASCII  characters  for  data  processing  of  the  main  pro- 
cessor . 

In  the  next  chapter  the  main  processor  design  is  dis- 
cussed and  the  system  status  is  described  with  an  example  of 
the  data  format  of  each  process. 


CHAPTER  5 

DESIGN  AND  IMPLEMENTATION  - II 
The  main  processor  design  and  its  implementation  are 
described  in  this  chapter.  Not  all  of  the  main  processor 
features  described  in  this  chapter  are  implemented  in  the 
present  system.  The  basic  frame  which  demonstrates  how  the 
expert  system  approaches  the  sleep  EEG  analysis  has  been  im- 
plemented. REM  period  detection  is  shown  as  an  example. 
The  system  design  concepts  describe  how  the  system  is  struc- 
tured and  how  the  system  works.  The  sleep  data  analysis  ex- 
ample and  the  current  system  status  are  discussed  in  detail 
at  the  end  of  this  chapter. 

Main  Pnoce^or  p^igr]  imcJ.mexi.U-U.as 
The  main  processor  reads  the  token  stream  and  generates 
interpreted  results  every  minute.  In  the  main  processor  the 
token  stream  is  regarded  as  a signal  sentence  which 
expresses  the  contained  information  in  a signal  language. 
The  signal  sentence  is  interpreted  using  a language  compiler 
mechanism  such  as  lexical  analysis,  syntactic  analysis,  and 
semantic  analysis.  The  interpreted  results  will  be  incor- 
porated with  the  subject  dependent  data  or  knowledge  source 
about  sleep  diagnosis  to  generate  a synthetic  report. 

In  the  present  system  only  the  interpreted  result  of 
the  sleep  EEG  is  utilized  in  the  synthetic  report  genera- 
tion. Though  different  mechanisms  can  be  utilized  for  the 
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syntactic  pattern  recognition,  a computer  compiler  mechanism 
is  utilized  for  system  implementation  since  the  automated 
procedure  of  compiler  generation  provides  several  useful 
features. 

For  instance,  the  standard  software  for  the  compiler 
generation  is  used  as  an  expert  system  development  tool  of 
the  present  sleep  EEG  expert  system.  The  separated  file 
structure  of  the  compiler  generation  makes  rule  management 
independent  in  the  sleep  EEG  expert  system.  This  gives 
great  flexibility  to  rule  implementation  and  makes  the 
inference  engine  independent  of  the  sleep  EEG  signal  under- 
standing rules.  In  addition,  the  automated  procedure  of 
compiler  generation  makes  system  development  rapid.  Whenev- 
er rule  adjustments  occur  in  the  system  refinement  stage  of 
expert  system  development,  the  whole  system  can  be  automati- 
cally regenerated. 

For  the  present  system  development,  a signal  compiler 
is  designed  for  the  channel  descriptor  and  a shell  program 
is  written  for  the  data  pipeline  of  the  channel  descriptor. 
The  lexical  part  of  a compiler  is  used  to  select  a proper 
waveform  token.  The  syntactic  and  semantic  part  of  a com- 
piler are  used  to  implement  the  sleep  EEG  signal  understand- 
ing rules.  Although  the  computer  compiler  mechanism  is  not 
sophisticated  enough  for  natural  language  processing,  a lim- 
ited number  of  specified  EEG  waveforms  make  a high  level 
language  compiler  mechanism  useful  for  EEG  signal  language 
understanding. 
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The  information  flow  of  the  main  processor 
Fig.  5.1.  The  main  processor  takes  inputs  from 
processor  and  the  operator. 
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Fig.  5.1  Information  Flow  of  Main  Processor 

The  necessary  information  for  setting  the  parameters  of  the 
interval  and  epoch  level  knowledge  sources  can  be  obtained 
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from  the  operator  at  the  beginning  of  the  session.  The 
parameters  set  previously  cannot  be  changed  interactively  in 
the  present  system.  Some  of  the  rule  parameters  of  the 
channel  descriptor  can  be  adjusted  by  the  knowledge  source 
procedures  utilizing  this  information.  After  the  channel 
description  rules  are  adjusted,  the  channel  description  list 
is  generated  from  the  output  token  stream  of  the  frontend 
processor.  Then  the  channel  description  list  is  used  for 
synthetic  report  generation.  In  synthetic  report  generation, 
if  the  patient  data  or  diagnosing  information  is  required, 
the  corresponding  knowledge  sources  are  activated  and  demand 
more  information  from  the  operator.  The  information  is  used 
when  the  decision  has  to  be  made  due  to  a lack  of  sufficient 
data.  Only  the  channel  description  list  is  utilized  for  syn- 
thetic report  generation  in  the  present  system. 

The  block  diagram  of  the  main  processor  is  shown  in 
Fig.  5.2.  The  main  processor  consists  of  a channel  descrip- 
tor, a synthetic  reporter,  and  a knowledge  source  about  di- 
agnosis (actual  diagnostic  rules  are  not  implemented  in  the 
present  system).  The  channel  descriptor  describes  channel 
activities  of  each  channel.  The  synthetic  reporter  gen- 
erates a synthetic  report  utilizing  the  channel  description 
list  and  the  knowledge  source  about  diagnosis.  The  design 
and  implementation  details  of  each  block  are  explained  in 
the  following  block  description  and  implementation. 
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Fig.  5.2  The  Block  Diagram  of  the  Main  Processor 

SJOSWSl  Des^iiiptipx 

The  channel  descriptor  describes  each  channel  activity 
in  terms  of  the  specified  EEG  waveforms.  For  instance,  if 
the  total  activity  time  of  the  alpha  wave  in  an  epoch  is 
larger  than  30  seconds  a channel  description  list  "alpha  ac- 
tivity in  channel  2 is  high"  is  generated.  Fig.  5.3  shows 
the  functional  structure  of  the  channel  descriptor,  consist- 
ing of  a decoder  and  several  translators.  Input  tokens  are 
decoded  into  a set  of  character  streams  in  the  decoder 
depending  on  the  corresponding  detector  outputs  of  each 
channel.  Every  waveform  detector  output  is  reexamined  in 
each  decoder  to  see  if  there  is  a false  detection  due  to 
other  channel  activities  and  if  so,  it  is  reconstructed  to 
denote  correct  waveform  information.  Each  channel  activity 
is  described  by  the  translator  using  the  decoded  output  of 
each  waveform  such  as  alpha,  beta,  delta,  etc.. 


Two  levels 
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of  rules  are  applied  in  this  process.  The  first  level  rules 
are  the  first  priority  rules  which  are  applied  independently 
to  the  token  itself  in  the  decoder.  The  second  level  rules 
describe  the  relationships  between  different  tokens. 
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Fig.  5.3  Functional  Structure  of  Channel  Descriptor. 

In  the  first  pass  of  the  channel  descriptor  first  level 
rules  are  applied  to  all  tokens  simultaneously.  All  the 
channels  are  tested  together  by  simultaneously  looking  at 
each  detector  output  horizontally.  That  is,  the  relation- 
ships between  different  channels  or  different  detector 
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outputs  of  the  same  channel  are  examined.  For  example,  if 
artifact  activity  exists  in  the  EOG  channel,  all  the  data  of 
the  same  interval  are  not  reliable.  In  that  case,  depending 
on  the  previous  data,  the  channel  information  has  to  be 
reexamined.  If  an  artifact  is  detected  then  the  delta  wave 
detector  output  of  the  same  interval  is  ignored  in  the  first 
pass  analysis.  Also,  if  there  is  a slow  wave  in  the  occipi- 
tal or  parietal  channel  during  the  same  interval,  the  slow 
wave  output  of  the  EOG  channel  is  ignored.  If  a spindle 
waveform  is  detected  with  the  artifact,  it  is  also  ignored 
and  so  on.  The  first  pass  analysis  outputs  are  the  consecu- 
tive token  stream  outputs  on  which  the  first  level  rules  al- 
ready have  been  applied  to  resolve  the  priority  between  the 
channels  or  the  detectors. 

In  the  second  pass,  the  second  level  rules  are  applied 
to  resolve  the  vertical  relationships  of  each  channel,  i.e., 
the  output  relationships  between  different  tokens.  The  re- 
lationships between  the  ouputs  of  the  same  detector  or  dif- 
ferent detectors  are  examined.  To  check  these  relations 
certain  rules  get  first  priority.  For  instance,  to  decide 
whether  the  detected  slow  waveform  output  of  the  EOG  channel 
is  truly  a REM  of  the  sleep  stage  1 or  not,  the  delta  wave 
detector  output  has  to  be  examined  first.  If  there  is  a del- 
ta wave  in  the  previous  or  the  following  two  seconds  then 
the  detector  output  of  the  EOG  channel  is  ignored.  The  to- 
tal activity  time  or  the  number  of  each  specified  waveform 
is  also  obtained  in  this  process.  For  the  total  activity 
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time  calculation  a smoothing  process  is  performed  on  the 
detector  output.  For  example,  in  visual  analysis  if  a 
discontinuous  interval  is  detected  in  the  alpha  waveform 
detector  output,  depending  on  the  duration  of  discontinuity 
it  may  be  regarded  as  a continuous  waveform  in  visual  total 
activity  time  calculation.  Three  quarter  second  smoothing 
is  used  for  slow  wave  (rslow,  oslow,  qslow,  pslow)  and  one 
second  smoothing  is  used  for  other  waveforms  in  the  present 
system. 

The  selected  word  representing  a specified  waveform  ac- 
tivity in  the  translator  depends  on  the  total  activity  time 
of  the  waveform  in  an  epoch.  Each  waveform  activity  time  or 
the  number  of  the  occurrence  for  the  corresponding  descrip- 
tion word  is  given  in  Fig.  5.4. 


high  medium  low  none 


alpha, 

30 

12 

2.5 

beta , 

30 

12 

2.5 

delta , 

30 

12 

2.5 

sigma , * 

2 

** 

1 

theta , 

30 

12 

2.5 

artifact, 

20 

1 0 

2.5 

rslow 

4 

2 

.25 

oslow , 

30 

1 2 

2.5 

pslow , 

30 

1 2 

2.5 

qslow, 

30 

1 2 

2.5 

* : number  of  occurrence,  others  are  seconds 
**:  ’’medium"  is  not  defined  for  sigma 


Fig.  5.4  Activity  Time  (number)  vs  Descriptive  Word 


To  describe  the  alpha  wave  activity,  the  total  alpha  time 
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per  one  minute  interval  is  calculated.  The  words  "high" 
refer  to  thoses  case  longer  than  30  seconds,  "medium”  to 
cases  between  12  seconds  and  30  seconds  and  "low"  for  cases 
less  than  12  seconds.  For  spindle  activity,  the  number  of 
occurrences  is  counted  rather  than  the  total  activity  time. 
The  activity  time  or  the  number  of  occurrences  of  a certain 
waveform,  which  are  used  for  decision  criteria,  may  be  ad- 
justed after  considering  the  applied  rules  for  the  particu- 
lar detector. 

CpanrjpJ.  Descri.pj.P£  InLpiprpppipiiop 

As  shown  in  the  previous  section  lexical  analysis,  syn- 
tactic and  semantic  analysis  are  required  for  the  EEG  signal 
language  understanding,  that  is  for  the  channel  descriptor. 
In  this  section  how  the  channel  descriptor  is  automatically 
generated  and  how  the  EEG  signal  understanding  rules  are  re- 
lated with  the  various  program  modules  for  channel  descrip- 
tor generation  are  discussed. 

Fig.  5.5  shows  a block  diagram  of  each  program  module, 
and  the  input/output  relationship  for  automatic  program  gen- 
eration for  the  channel  descriptor.  "Lex"  and  "Yacc"  are 
standard  software  library  packages  on  the  UNIX  operating 
system  which  supports  the  automatic  design  of  language  com- 
pilers. "Main"  is  a hand  written  program  in  "C"  language  for 
calling  "yyparse"  which  is  a parsing  function.  "Lex"  [LES75] 
is  the  program  which  generates  a lexical  analyzer  whose  con- 
trol is  directed  by  instances  of  regular  expressions  in  the 
input  stream.  To  generate  a lexical  analyzer  using  the  Lex 
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program,  lexical  rules  are  provided  as  a lex  source.  The  lex 
source  (lexical  rules)  consists  of  two  parts:  a set  of  regu- 
lar expressions  and  the  corresponding  program  fragments 
which  derive  lex  action.  A set  of  lexical  rules  is 
translated  to  a program  which  reads  an  input  stream,  coping 
it  to  an  output  stream  and  partioning  the  input  stream  which 
matches  the  given  expression.  As  each  such  string  is  recog- 
nized the  corresponding  program  fragment,  such  as  lex  ac- 
tions, are  executed.  The  recognition  of  the  expression  is 
performed  by  a deterministic  finite  automation  generated  by 
’Lex*.  In  the  present  system  'Lex*  generates  a decoder  pro- 
gram which  reads  the  input  token  stream  and  converts  it  into 
an  output  token  stream  (yylex  in  Fig.  5.5)  for  which 
waveform  priority  rules  are  already  resolved. 
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Fig  5.5  Block  Diagram  for  Automatic  Program  Generation 
for  Channel  Descriptor. 
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In  Lex  source  all  the  specific  waveform  tokens  are  de- 
fined in  the  regular  expression  and  the  waveform  priority 
resolution  rules  are  implemented  in  lex  action  program  seg- 
ments. More  than  60  nonterminal  tokens  are  defined  in  a lex 
source;  namely  10  waveform  tokens,  40  corresponding  wave 
pattern  tokens  and  others.  Four  rules  (seventeen  cases)  are 
implemented  in  the  Lex  action  program  fragments  of  the  first 
pass.  More  than  40  rules  are  implemented  in  the  second  pass 
(see  Appendix  3 for  all  source  input  files,  grammar  files, 
parser  and  semantic  action  files). 

',Yacc,,  [Jo75]  is  a program  which  takes  the  language 
grammar  and  corresponding  semantic  actions  as  input  source 
and  generates  a parsing  tree.  Yacc  reads  yylex  as  input 
from  Lex  and  takes  semantic  action  when  specific  rules  are 
recognized.  The  rules  are  specified  in  the  grammar  in  which 
the  input  structure  is  described  by  the  user.  The  actions 
are  described  in  the  semantic  action  file. 

The  input  source  programs  in  which  the  channel  activity 
is  described  are  the  key  elements  in  the  translator  of  the 
channel  descriptor.  In  the  present  system  context  free 
grammar  [Ah79l  [Ho79l  is  used  to  read  various  EEG  signal 
waveform  inputs.  Different  grammars  can  be  used  depending 
on  the  signal  characteristics  of  the  waveforms  to  be  select- 
ed. The  channel  description  rules  are  implemented  in  a se- 
mantic action  file.  To  describe  each  waveform  activity 
several  counters  are  kept  in  each  semantic  action  file.  One 
of  the  counters  keeps  track  of  the  continuous  off  time  of  a 
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waveform  detector  and  is  then  used  for  a waveform  smoothing. 
Another  counter  is  used  to  determine  the  total  activity 
time  of  an  epoch  that  describes  waveform  activity.  Whenever 
the  specified  input  token  is  recognized  the  counter  is  in- 
creased by  one.  A certain  waveform  counter  is  used  to  check 
different  waveform  activity  which  is  related  to  its  own 
waveform.  For  instance,  in  the  rslow  semantic  action  file 
three  counters  are  used  to  check  the  other  channel’s  slow 
activities.  They  check  the  delta  activity  in  the  previous 
and  the  following  two  second  intervals  or  slow  wave  activity 
in  the  other  channels  to  see  whether  a slow  wave  on  the  EOG 
channel  is  truly  a REM  wave  or  not. 

In  the  translator  implementation  of  the  present  system 
the  ”Slrgen"  [Lo84]  and  ’’Parser”  programs  are  used  instead 
of  Yacc.  Yacc  consists  of  a parsing  table  generator  and  a 
parser  which  handles  parsing  actions  such  as  shift,  reduce, 
accept,  and  error  action.  Slrgen  is  a simple  left  right 
parsing  table  generator  which  does  not  accept  ambiguous 
grammar.  Parser  is  a hand  written  program  to  handle  the 
parsing  table  in  "C"  language.  The  separated  program  Slrgen 
and  Parser  give  more  flexibility  to  handle  the  various  rule 
implementations.  However,  to  generalize  the  system  genera- 
tion in  UNIX  operating  system  Yacc  has  to  be  used. 

The  actual  pass  of  the  each  waveform  descriptor  is 
shown  in  Fig  5.6.  Each  channel  descriptor  has  a common  pass 


of  Lex  and  its  individual  pass  of  Lex,  Slrgen,  and  Parser. 
Although  Lex,  Slrgen,  and  Parser  in  the  second  pass  could  be 
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combined  into  a single  Lex,  Slrgen  and  Parser  pass,  this 
structure  is  selected  to  give  great  flexibilities  in  in- 
dependent rule  management  for  each  waveform  description. 
The  corresponding  channel  description  rules  are  evoked 
depending  on  the  waveform  in  each  channel.  For  automatic 
program  generation,  a shell  command  program  of  UNIX  [B078] 
[B083]  operating  system  is  written  in  a Makefile  (see  Appen- 
dix 3).  All  the  file  management  procedures  are  shown  in  this 
file. 
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Fig.  5.6  Program  Modules  used  for  Channel  Descriptor. 

Some  of  the  rules  are  coded  into  a grammar  of  the  sig- 
nal specification  language  and  some  of  the  rules  are  coded 


- 94 


into  the  semantic  action  part  of  the  language  compiler.  For 
instance,  alpha  period  discrimination  is  encoded  in  the 
grammar  file  and  the  total  alpha  time  smoothing,  alpha  time 
calculation  and  the  description,  are  encoded  in  the  semantic 
action  file.  In  the  present  system  most  of  the  channel 
description  rules  are  implemented  in  semantic  action  file. 

Although  it  is  not  clear  where  the  rules  have  to  be 
coded  into  and  although  the  rule  implementation  procedures 
are  complex,  due  to  the  frequent  rule  changes  in  expert  sys- 
tem development,  standard  software  provides  an  excellent 
tool  for  designing  the  present  expert  system.  Specifically, 
when  the  waveform  in  the  signal  data  is  well  defined  it  can 
be  considered  a problem  of  understanding  the  simple  signal 
language.  The  prominent  signal  shape  to  the  human  EEGer  can 
be  regarded  as  a word  in  a signal  language.  The  interrela- 
tionships of  each  waveform  on  various  channels  can  provide  a 
proper  grammar  and  semantic  action  for  interpreting  the  EEG 
signal.  All  the  files  for  generating  the  channel  descriptor 
are  given  in  Appendix  3. 

Syntijgtic  R^oxi^r 

The  synthetic  reporter  generates  a synthetic  report 
based  on  the  EEG  signal  analysis  results  (which  is  now  a 
channel  description  list),  subject  dependent  data,  and  a 
knowledge  source  about  diagnostic  rules  (which  are  not  im- 
plemented yet).  The  synthetic  reporter  will  be  the  heart  of 
the  system.  The  channel  description  list  is  one  of  the 
knowledge  sources  in  the  synthetic  reporter.  The  present 
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system  only  contains  the  REM  period  detection  knowledge. 

The  structure  of  the  synthetic  reporter  has  to  be 
dependent  on  the  knowledge  source  to  be  used  because  the 
knowledge  structure  will  affect  the  structure  of  inference 
engine  and  the  reasoning  mechanism  of  the  system.  Although 
the  structure  of  the  sleep  disorder  diagnosing  rules  are  not 
thoroughly  studied  and  cannot  be  clearly  defined,  the  frame 
design  of  the  synthetic  reporter  is  described  to  see  how  the 
channel  description  list  can  be  used  in  the  synthetic  re- 
porter . 

Fig.  5.7  shows  the  decision  tree  of  a synthetic  report- 
er. 
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Fig.  5.7  Decision  Tree  of  Synthetic  Reporter 
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User  supplied  information  forms  the  input  to  the  screening 
modules.  A screening  module  is  first  selected  depending  on 
the  analysis  required.  The  REM  screening  module  is  the 
first  one  in  the  present  system.  After  a screening  module  is 
selected,  a set  of  production  rules  in  the  screening  module 
gives  a best  selection  of  a frame. 

For  example,  depending  on  the  required  analysis  screen- 
ing module  1 (REM)  can  select  a proper  set  of  rules  to  in- 
terpret the  channel  description  list  (i.e  rules  for  the  al- 
coholics or  normal  subjects).  In  this  level  selection,  the 
branching  ambiguity  can  be  minimized  by  proper  grouping  of 
the  goal  states  without  overlapping.  Once  a frame  is  select- 
ed, all  the  information  is  utilized  to  select  a final  goal 
state.  A synthetic  report  is  then  generated  based  on  the 
information  of  the  final  goal  state. 

For  the  reasoning  process  both  of  the  directions  of  the 
reasoning  need  to  be  used  depending  on  the  level  of  the 
tree.  In  the  lower  level  of  the  decision  tree,  a branching 
decision  is  made  using  some  of  the  key  information  of  a 
channel  description  list  and  user  supplied  information.  In 
the  upper  level  of  the  tree,  i.e.,  the  subtree  which  has  a 
group  of  final  goal  states,  a model  driven  or  reasoning 
backward  strategy  is  used  to  score  every  member  of  the  group 
such  as  a clinical  name  of  a sleep  disorder.  A confidence 
factor  also  can  be  given  to  every  property  in  the  property 
list,  i.e.,  a certain  symptom  of  sleep  disorder.  If  there 
is  more  than  one  member  which  has  a score  over  threshold 
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value,  then  all  the  selected  members  will  be  displayed  in 
decreasing  order.  The  selection  of  the  confidence  factor  of 
the  property  and  the  threshold  value  for  decision  making  are 
based  on  both  the  human  expert  knowledge  and  the  other  sys- 
tem parameters.  This  structure  was  selected  to  aid  future 
system  expansion  which  will  include  the  analysis  of  sleep 
disorders.  The  synthetic  module  is  written  in  ’LISP' 
language  to  handle  symbolic  data  and  natural  language  like 
user  interface. 

KnowJ,e4_g£ 

Generally,  problems  in  different  domains  need  different 
knowledge  base  structures  and  knowledge  base  management 
schemes.  Though  there  exists  no  general  theory  of  knowledge 
representation,  depending  on  the  application  and  the  kind  of 
knowledge,  a good  combination  of  declarative  structure  and 
interpretive  procedure  [ R i8 3 3 can  be  utilized  for  proper  ac- 
quisition and  retrieval  of  a knowledge  base. 

In  the  present  system,  the  procedural  representation 
method  is  used  to  handle  the  channel  descriptor.  A hybrid 
system  of  a frame  scheme  and  a production  system  is  selected 
to  process  the  channel  description  list  and  the  diagnosing 
rules  about  sleep  disorder.  A set  of  production  rules  is 
used  in  screening  modules  as  branching  rules.  For  instance, 
for  a frame  selection,  "If  REM  analysis  is  required  for  the 
normal  subject  then  select  frame  1"  is  used.  General  signal 
characteristics  of  each  channel  plus  some  of  the  user-  sup- 
plied information  constitute  the  condition  part  of  the 
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production  rules.  A screening  module  contains  a certain 
number  of  frames.  A primitive  semantic  network  in  the  frame 
structure  is  used  to  represent  each  group  of  analysis  or 
diseases.  A frame  can  have  different  level  structure 
depending  on  the  classification  level.  A frame  contains  a 
group  of  final  states  which  share  some  common  properties. 
Each  final  state  has  a property  list  which  describes  a goal 
state.  Basically,  a generalized  property  list  is  implement- 
ed in  a frame  structure.  Each  symptom  of  a sleep  disorder 
can  be  described  in  the  property  list.  If  there  is  a slot 
of  unknown  property,  information  can  be  obtained  from  the 
users  in  hypothesis  assurance  session. 

The  present  system  has  no  distinct  characteristics  of 
the  knowledge  representation  methods.  Depending  on  the 
knowledge  base  which  will  be  included  in  the  future  and  the 
reasoning  mechanism,  the  knowledge  representation  techniques 
need  to  be  refined. 
tolajia_y,_er)  pi  tjae  Rs&s&aijig 

The  explanation  of  the  reasoning  is  one  of  the  prom- 
inent features  in  the  expert  system.  Although  a good  reason- 
ing mechanism  usually  slows  down  the  system,  it  is  necessary 
for  communication  with  the  user.  There  has  to  be  a tradeoff 
between  good  explanation  capability  and  the  system  speed. 
The  explanation  mechanism  is  heavily  dependent  upon  the  sys- 
tem knowledge  structure  and  the  language  used  for  building 
the  system. 
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In  the  present  system,  only  a primitive  explanation  ca- 
pability is  provided  in  the  synthetic  reporter.  Whenever  a 
branching  occurs  in  the  synthetic  reporter,  the  branching 
node  is  stacked  in  the  explanation  stack.  If  the  question 
is  given,  then  the  stack  is  dumped  from  the  top,  depending 
on  the  level  of  the  question.  The  level  of  the  answer  is 
the  same  as  the  level  of  the  stack.  The  current  explanation 
mechanism  needs  to  be  enhanced  to  include  the  reasoning  base 
of  the  sleep  disorder  in  the  future  system  expansion. 

o£  Pj-o.sps.sijxg 

Some  examples  of  the  sleep  data  processing  done  by  the 
system  are  given  in  this  section.  Fig.  5.8  shows  the  data 
format  in  each  step  of  the  frontend  processor.  (a)  shows 
one  of  the  four  channel  inputs  (occipital  channel)  of  10 
seconds.  (b)  is  the  same  channel  analog  output  after  pass- 
ing through  a linear  phase  FIR  alpha  filter.  Eight  filters 
are  assigned  to  four  input  channels  . (c)  shows  an  alpha 
wave  counter  output  which  counts  the  number  of  alpha  waves 
using  wave  criteria.  (d)  shows  an  alpha  wave  detector  out- 
put which  displays  a certain  amount  of  time  depending  on 
each  waveform  activity.  (e)  shows  10  seconds  of  decoded 
sampled  outputs  of  each  waveform  detector.  A four  hexade- 
cimal number  contains  a quarter  second  channel  information. 
The  first  10  bits  of  four  hexadecimal  number  denotes  each 
waveform  existence  and  the  11th  bit  is  used  for  the  minute 


indicator . 
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Fig.  5.9  shows  the  data  format  in  each  step  of  the  main 
processor.  (a)  shows  a preprocessed  data  format  of  10 
seconds  of  the  decoder  output  in  the  channel  descriptor.  In 
the  first  pass  analysis,  the  priority  of  the  detector  is 
resolved.  The  encircled  number  shows  an  artifact  priority 
over  a delta  wave,  which  ignores  the  delta  wave  detector 
output  if  the  artifact  exists  in  the  eye  channel.  (b)  is  a 
20  minutes  second  pass  analysis  output  of  the  translator. 
Each  row  denotes  a different  waveform  activity  and  each 
column  denotes  one  minute  waveform  activities.  The  charac- 
ters "H",  nM" , "L",  and  "N"  are  used  to  denote  "high", 
"medium",  "low",  and  "none",  respectively.  (c)  is  the  final 
output  of  the  synthetic  reporter.  Sleep  staging  rules  and 
heuristic  REM  detection  rules  are  implemented  in  the  REM 
module  of  the  synthetic  reporter.  Using  a four  minute  mov- 
ing window  REM  sleep  stage  is  detected.  A minute  state 
description  shown  in  example  is  available  depending  on  user 
request. 
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Fig.  5.8  Data  Format  in  Token  Generator 
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Fig.  5.9  Data  Format  in  Token  Processor 
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System  Status 

The  present  EEG  expert  system  is  in  a prototype  stage. 
The  whole  system  evaluation  is  not  possible  without  rule  re- 
finement. The  REM  period  detection  module  is  currently  be- 
ing tested.  The  intermediate  results  of  the  REM  period 
analysis  are  discussed  in  the  next  chapter.  It  demonstrated 
a good  result  for  limited  data.  More  data  have  to  be 
analyzed  to  refine  REM  period  detection  rules.  The  frontend 
processor,  which  consists  of  filters,  detectors,  and  a mul- 
tiplexor, is  in  refinement  stage.  The  main  processor,  which 
consists  of  a channel  descriptor  and  a synthetic  reporter, 
is  not  complete.  Although  the  channel  descriptor  is  in  a 
refinement  stage,  the  synthetic  reporter  is  far  from  com- 
plete. To  complete  the  synthetic  reporter  design,  the 
knowledge  structure  of  the  sleep  disorder  diagnosing  rules 
first  have  to  be  defined  and  the  knowledge  representation 
techniques  need  to  be  used  properly. 

For  the  frontend  processor  refinement,  each  module, 
such  as  filter,  detector,  and  multiplexor,  needs  to  be  test- 
ed. All  the  responses  affect  the  zero  crossing.  In  the 
filter  design,  the  effects  of  linear  phase  intercept  distor- 
tion, phase  linearity,  and  skirt  characteristics  to  zero 
crossing  have  to  be  studied  more  quantitatively  and  have  to 
be  considered  in  the  filter  design.  Generally,  the  zero  in- 
tercept distortion  condition  gives  constraints  to  the  zero 
positioning  in  the  linear  phase  filter  design.  Also,  the 
zero  positioning  of  the  linear  phase  filter  affects  the 
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amplitude  characteristics  of  the  filter.  In  the  detector  re- 
finement, the  waveform  criteria  and  the  amplitude  criteria 
need  to  be  adjusted  to  closely  simulate  visual  analysis. 
This  can  only  be  achieved  by  analyzing  a lot  of  data  with 
different  criteria.  For  the  multiplexor,  the  sampling  rate 
for  token  generation  needs  to  be  considered  to  decrease  the 
use  of  on-board  memory  and  for  the  compact  data  acquisition. 

In  the  channel  descriptor  of  the  main  processor,  the 
signal  grammar  needs  to  be  clearly  defined.  Although  the 
current  grammar  takes  any  kind  of  waveform  token,  more  syn- 
tactic patterns  can  be  simply  recognized  by  a proper  gram- 
mar. This  will  decrease  the  semantic  actions  in  the  present 
system  and  can  provide  the  capability  of  erroneous  data  han- 
dling. In  the  synthetic  reporter,  the  knowledge  base  is 
currently  composed  of  only  the  channel  list.  The  connection 
link  between  the  channel  description  list  and  sleep  disorder 
is  not  provided  yet.  In  the  synthetic  reporter  design, 
first  the  task  definition  has  to  be  clearly  defined.  Once 
the  task  is  defined  the  synthetic  reporter  can  be  designed 
using  a pseudo  channel  description  list. 

System  expansion  demands  an  extensive  examination  of 
the  knowledge  structure.  As  the  knowledge  base  grows  with 
system  expansion,  the  knowledge  representation  method, 
inference  engine,  and  the  explanation  of  the  reasoning  base 
have  to  be  restructured.  For  the  data  base  management  the 
binary  data  transfer  and  store  program  need  to  be  developed 
in  the  present  system.  A binary  file  handling  capability 
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will  provide  a global  data  pool  for  the  expert  system. 
For  the  user  friendly  system  the  controller,  which  provides 
the  interactive  module,  control  needs  to  be  added.  Currently 
the  interactive  rule  management  or  system  parameter  adjust- 
ments capability  is  not  provided.  The  system  should  be 
enhanced  with  an  interactive  rule  management  module.  This 
module  will  enable  the  system  to  explain  the  reasoning  base 
at  any  state  and  allow  the  user  to  easily  access  the 
knowledge  base  in  order  to  add  and  modify  the  rules. 


CHAPTER  6 

SYSTEM  EVALUATION  AND  CONCLUSIONS 

Today  the  system  is  in  the  prototype  stage.  The  present 
system  evaluation  emphasis  is  qualitative  rather  than  quan- 
titative. As  an  initial  refinement,  the  basic  system  rules 
were  evaluated  and  modified  corresponding  to  the  analysis 
results.  The  rule  modifications  are  performed  on  the  rule 
files  and  the  the  system  is  automatically  regenerated  and 
evaluated.  The  simplicity  of  the  rule  modification  is 
essential  for  an  expert  system  development.  The  system  re- 
finements need  flexible  rule  management  and  supporting 
tools. 

Although  it  is  difficult  to  check  the  validity  and  re- 
liability of  the  present  system  as  there  is  no  standard  data 
base  with  which  to  compare  the  results,  the  REM  period 
detection  of  the  present  system  is  evaluated  to  show  how  the 
system  performs  and  can  be  refined,  by  comparing  the  results 
obtained  from  the  present  system  with  those  obtained  from 
the  visual  analysis  minute  by  minute.  The  results  demon- 
strate that  the  analyzing  method  of  the  present  system  can 
interpret  the  sleep  data  as  much  in  the  same  way  as  the 
visual  analysis  method  and  that  the  rule  modifications  are 
properly  handled  in  the  present  system. 
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One  subject  was  selected  in  each  of  the  following  age 

groups:  Groupl : 25  to  34  years,  Group2:  43  to  53  years, 

Group3:  67  to  79  years.  Six  hours  sleep  data  of  each  sub- 
ject were  processed  off-line.  A SANGAMO  3500  recoder  with 
sixteen  channels  and  a GRASS  Model  78  polygraph  system  were 
utilized  to  replay  and  display  the  sleep  data  prerecorded  on 
the  magnetic  tapes.  The  recoder  output  was  connected  to  SAHC 
[Sm78],  The  gain  of  the  SAHC  external  output  of  each  chan- 
nel was  adjusted  to  1 Volt  peak  to  peak  for  50  uV  EEG/EOG 

input.  The  external  outputs  of  the  SAHC  were  connected  to 
EEG  polygraph  system  and  the  A/D  converter  input  of  the  sys- 
tem. 

The  rules  used  for  the  channel  description  are: 

1)  "N",  "L",  ”M",  and  "H"  denote  the  total  corresponding  ac- 
tivity time  range  (see  Fig.  5.4)  in  an  epoch. 

2)  The  alpha  activity  in  presence  of  a delta  wave  is  not 
counted  for  the  total  alpha  activity  time  calculation. 

3)  The  slow  wave  on  the  EOG  channel  with  a delta  activity  is 
not  regarded  as  a REM. 

4)  The  slow  wave  on  the  EOG  channel  with  a slow  wave  on  the 
occipital  or  central  channel  is  not  regarded  as  a REM. 

5)  A slow  wave  which  has  a delta  wave  in  the  previous  or  the 
following  2 seconds  is  not  counted  for  the  total  REM  activi- 
ty time  calculation. 

6)  The  tokens  with  artifact  activity  are  ignored. 

The  REM  period  detection  rules  applied  to  the  channel 


description  list  are: 
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1)  If  alpha  activity  is  H (larger  than  30  seconds)  in  an 
epoch,  it  is  not  regarded  as  a REM  period. 

2)  An  epoch  surrounded  by  H alpha  minutes  is  not  regarded  as 
a REM  period. 

3)  If  delta  activity  in  an  epoch  is  equal  to  or  greater  than 
L (larger  than  2.5  seconds),  the  epoch  is  not  regarded  as  a 
REM  period. 

4)  An  epoch  surrounded  by  M delta  minutes,  or  H delta 
minutes,  or  M and  H delta  minutes  is  not  regarded  as  a REM 
period. 

5)  If  the  slow  wave  activity  (0.3  HZ-  6.0Hz)  in  an  epoch  is 
L,  M or  H(larger  than  2.5  seconds)  in  occipital  or  parietal 
channel,  it  is  not  regarded  as  a REM  period. 

6)  An  epoch  which  has  L,  M,  or  H REM  is  regarded  as  a REM 
period  with  the  exception  of  the  REM  period  detection  rules 
from  1 to  5 . 

7)  A period  between  REM  periods  is  also  regarded  as  a REM 
period. 

Fig.  6.1  shows  the  intermediate  results  of  the  system. 
The  upper  row  and  lower  row  denote  non-REM  period  data  and 
REM  period  data,  respectively.  The  first  element  in  each 
four  tuple  denotes  the  total  non-REM  or  REM  minutes  in  one 
hour  data.  The  second  element  denotes  the  total  correct 
detection  minutes.  The  third  element  denotes  the  total 
missed  minutes  during  detection.  The  fourth  element  denotes 
the  total  false  detection  minutes. 
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Table  6.1 


Data  Number: 
Hour 
1 st 

10071 

(60,58,2,0) 
(0,0, 0,2) 

11735 

(60,60,0,0) 
(0,0, 0,0) 

11747 

(60,59,0,0) 
(0,0,0, 1 ) 

2nd 

(58,58,0,1) 
(2,1 ,1 ,0) 

(60,56,4,0) 
(0,0, 0,4) 

(43,42,1,3) 
(17,14,3,1 ) 

3rd 

(55,55,0,3) 
(5, 2, 3,0) 

(52,52,0,3) 

(13,10,3,0) 

(45,45,0,2) 

(15,13,2,0) 

4 th 

(53,53,0,3) 
(7, 4, 3,0) 

(45,42,3,7) 

(15,8,7,3) 

(55,55,0,1 ) 
(5, 5, 0,1) 

5 th 

(29,18,1  ,4) 
(31,27,4,1) 

(51  ,51  ,0,6) 
(9, 3, 6,0) 

(34,32,2,5) 
(26,21 ,5,2) 

6 th 

(44,44,0,3) 

(16,13,3,0) 

(45,45,0,15) 

(15,0,15,0) 

(53,53,0,2) 
(7, 5, 2,0) 

Total  for 
REM  period 

(61 ,47,14,3) 

(52,21 ,31,7) 

(70,58,12,5) 

In  the  first  run  of  the  data,  a 77%  (10071),  40% 
(11735)  and  83%  (11747)  true  REM  period  were  obtained  in  the 
360  minutes  record  for  each  subject.  In  the  first  data 
analysis,  the  following  points  were  observed. 

1)  The  number  of  REM  clusterings  is  mostly  detected  correct- 
ly in  all  the  age  groups. 

2)  Most  REM  periods  not  detected  bu  the  system  are  due  to 
the  channel  descriptor  and  the  synthetic  reporter.  The  rule 
which  ignores  the  previous  and  the  following  2 seconds  of  a 
delta  wave  and  the  rule  which  does  not  regard  a slow  wave  on 
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EOG  channel  with  a delta  wave  as  a REM  (channel  description 
rule  5 and  3)  incorrectly  removes  REM  waves.  The  rule  which 
does  not  regard  a REM  minute  when  the  total  delta  activity 
time  is  larger  than  2.5  seconds  incorrectly  removes  REM 
minutes.  (9  minutes,  2156  in  data  1 1735  of  Group3)  . When  a 
clear  REM  period  cluster  is  detected,  two  minute  REM  period 
smoothing  can  decrease  the  number  of  the  missed  REM  periods 
almost  without  increasing  false  alarms.  One  subject  (11735) 
had  more  than  30  seconds  alpha  time  with  clear  REM  waves.  In 
this  subject,  the  high  alpha  activity  (channel  description 
rule  1)  removes  REM  period  16  minutes  (31 %)• 

3)  The  optimum  threshold  value  for  each  detector  is  not  easy 
to  determine.  Once  the  threshold  value  is  changed,  all  the 
tokens  need  to  be  regenerated,  i.e.,  the  data  has  to  be  run 
again  on-line.  For  setting  up  the  threshold  values,  four 
different  records  in  different  age  group  were  utilized. 

4)  The  results  are  sensitive  to  delta  wave  activities.  The 
dropout  of  the  tape  and  bad  electrode  contacts  of  (11735) 
generates  false  detections  of  the  delta  waves,  which  in  turn 
removes  the  REM  waves  in  the  REM  period  detection. 

The  problems  described  previously  indicate  more  rules 
need  to  be  added  and  refined.  For  instance,  the  rules  for 
the  delta  wave  and  artifact  need  to  be  changed  and  evaluated 
again  utilizing  the  generated  tokens.  The  channel  descrip- 
tion list  can  also  be  smoothed  to  get  rid  of  false  detec- 
tions, i.e.,  a one  minute  false  REM  period  detection. 
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The  same  data  tokens  were  analyzed  using  different 
rules.  All  the  rule  changes  have  been  made  with  the  lex- 
specs  or  semantic  action  file  of  the  corresponding  channel 
descriptor.  After  changing  the  rules,  the  system  is  easily 
regenerated  by  executing  a command  file.  The  simplicity  of 
the  rule  changes  or  the  system  regeneration  gives  great 
flexibility  to  the  present  system  rule  management. 

Table  6.2  shows  the  results  after  removing  the  rules 
which  ignore  the  tokens  with  artifact  and  which  do  not  re- 
gard a slow  wave  on  the  EOG  channel  as  a REM  when  there  ex- 
ists a slow  wave  on  the  occipital  or  parietal  channel. 
Table  6.4  shows  the  results  after  changing  the  effective  REM 
activity  calculation  (rule  5 of  the  channel  description) 
range  from  2 seconds  to  1 second. 
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Table  6.2 


Data  Number: 
Hour 
1 St 

10071 

(60,58,2,0) 
(0,0, 0,2) 

11735 

(60,60,0,0) 
(0,0, 0,0) 

11747 

(60,59,1,0) 
(0,0, 0,1 ) 

2nd 

(58,58,0,1 ) 
(2,1,1 ,0) 

(60,56,4,0) 
(0,0, 0,4) 

(43,42,1,3) 

(17,14,3,1) 

3rd 

(55,55,0,2) 
(5, 3, 2,0) 

(47,46,1,3) 

(13,10,3,1) 

(45,45,0,1 ) 
(15,14,1,0) 

4 th 

(53,53,0,3) 
(7, 4, 3,0) 

(45,42,3,7) 

(15,8,7,3) 

(55,55,0,0) 
(5, 5, 0,0) 

5 th 

(29,28,1 ,4) 
(31,27,4,1) 

(51 ,51 ,0,5) 
(9, 4, 5,0) 

(34,32,2,4) 

(26,22,4,2) 

6 th 

(44,44,0,2) 

(16,14,2,0) 

(45,45,0,15) 

(15,0,15,0) 

(53,53,0,0) 
(7, 7, 0,0) 

Total  for 
REM  period 

(61 ,49,12,3) 

(52,22,30,8) 

(70,62,8,4) 
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Table  6.3 


Data  Number: 
Hour 
1 st 

10071 

(60,58,2,0) 
(0,0, 0,2) 

11735 

(60,60,0,0) 
(0,0, 0,0) 

11747 

(60,59,1 ,0,0) 
(0,0, 0,1 ) 

2nd 

(58,58,0,0) 
(2, 2, 0,0) 

(60,56,4,0) 
(0,0, 0,4) 

(43,42,1,3) 

(17,14,3,1) 

3rd 

(55,55,0,2) 
(5, 3, 2,0) 

(47,46,1  ,2) 
(13,11,2,1) 

(45,45,0,1 ) 
(15,14,1,0) 

4 th 

(53,53,0,2) 
(7, 5, 2,0) 

(45,42,3,5) 

(15,10,5,3) 

(55,54,1 ,0) 
(5, 5, 0,1) 

5 th 

(29,28,1 ,2) 
(31,29,2,1) 

(51 ,51  ,0,3) 
(9, 6, 3,0) 

(34,32,2,4) 

(26,22,4,2) 

6 th 

(44,44,0,2) 

(16,14,2,0) 

(45,45,0,14) 

(15,1,14,0) 

(53,53,0,0) 
(7, 7, 0,0) 

Total  for 
REM  period 

(61 ,53,8,3) 

(52,28,24,8) 

(70,62,8,5) 
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The  Table  6.4  shows  the  results  of  the  rule  change  in 
the  system. 


Table  6.4 


1st  run 
2nd  run 
3rd  run 


d 1 007  1 

(61 .47.14.3) 

(61 .49.12.3) 

(61.53.8.3) 


d 1 1735 
(52,21 ,31  ,7) 
(52,22,30,8) 
(52,28,24,8) 


dl 1747 

(70.58.12.5) 

(70.62.8.5) 

(70,62,8 ,5) 


The  results  obtained  from  the  Table  6.5  shows  87?,  54%,  and 
89%  HEM  period  detection.  Although  the  results  are  not  sa- 
tisfactory, it  was  clearly  demonstrated  that  the  system  per- 
formance level  can  be  increased  as  the  system  rules  grow  and 
are  refined. 

The  basic  system  frame  design  and  REM  module  are  com- 
pleted. Although  it  is  far  from  a final  form,  it  is  shown 
that  the  data  driven  system  based  on  the  time  domain 
waveform  detection  method  (Sra78)  and  heuristic  rules  for 
symbolic  signal  processing  are  feasible,  especially  for 
simulating  the  gestalt  perception  of  low  frequency  biomedi- 
cal signals.  The  experience  of  building  an  expert  system 
has  shown  that  software  for  design  aids  are  essential  for 
the  expert  system  development,  especially  for  the  frequently 
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changing  rule  implementation.  Further  system  development, 
particularly  automatic  waveform  detector  generation  and 
proper  knowledge  base  management  in  synthetic  reporter  will 
enhance  the  system  to  be  used  for  biomedical  signal  analysis 
applications. 


APPENDIX  1 

TOKEN  GENERATOR  PROGRAM 


The  main  part  of  the  token  generator  program  is  given 
to  show  the  structure  of  the  frontend  processor.  The  fol- 
lowing is  the  TI-9900  assembly  program  of  interrupt  handler 
(HANDLER),  alpha  filter  (ALPHA),  alpha  detector  (ALPHADET), 
and  main  (FRONT).  The  other  parts  of  the  program  (the  other 
filters  and  detectors)  are  available  in  the  EEG  laboratory 
of  Electrical  Engineering,  University  of  Florida. 
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PDOS  ASM  R2.4 

14:26  09/14/84  FILE:  HANDLER / 1 , SLEEP 

A 

* MAIN  CONTROL  OF  DATA  SAMPLING  AND  AQUISITION 

* LOW  PASS  FILTERING  EVERY  CHANNEL  WITH  ZERO 

* AT  -1  FOR  480HZ  SAMPLING  RATE  TO  ELIMINATE 

* ALIASING  PROBLEM 

A 


IDT 

'HANDLER' 

CFF8 

MUX 

EQU 

>CFF8 

CHANNEL  MULTIPLEXING 

CFFC 

EOC 

EQU 

>CFFC 

END  OF  CONVERSION 

CFFA 

CON 

EQU 

>CFFA 

A/D  CONDERSION  STARTER 

CFFE 

ADC 

EQU 

>CFFE 

OUTPUT  OF  A/D  CONVERTER 

F088 

AEX 

EQU 

>F088 

ALPHA  EXIST 

F288 

BEX 

EQU 

>F288 

BETA  EXIST 

F488 

DEX 

EQU 

>F488 

DELTA  EXIST 

F688 

SEX 

EQU 

>F688 

SIGMA  EXIST 

F888 

TEX 

EQU 

>F888 

THETA  EXIST 

FA88 

FEX 

EQU 

>FA88 

ARTIFACT  EXIST 

15B8 

OEX 

EQU 

> 15B8 

OREM  EXIST 

17B8 

PEX 

EQU 

> 17B8 

PREM  EXIST 

19B8 

QEX 

EQU 

> 19B8 

QREM  EXIST 

1BB8 

REX 

EQU 

> 1BB8 

REM  EXIST 

11F0 

CHO 

EQU 

> 11F0 

CHO  OUTPUT  PORT 

11F2 

CHI 

EQU 

> 11F2 

CHI  OUTPUT  PORT 

11F4 

CH2 

EQU 

> 11F4 

CH2  OUTPUT  PORT 

11F6 

CH3 

EQU 

> 11F6 

CH3  OUTPUT  PORT 

11CF 

CHOO 

EQU 

> 11CF 

DELAY  UNIT  FOR  CHO  LPF 

11CD 

CH11 

EQU 

> 11CD 

DELAY  UNIT  FOR  CHI  LPF 

11CB 

CH22 

EQU 

> 11CB 

DELAY  UNIT  FOR  CH2  LPF 

11C9 

CH33 

EQU 

> 11C9 

DELAY  UNIT  FOR  CH3  LPF 

11F8 

CLOCK 

EQU 

> 11F8 

CLOCK 

UFA 

NDEL 

EQU 

> UFA 

# OF  DELTA  IN  A MIN. 

11FC 

NSPI 

EQU 

> 11FC 

# OF  SPINDLE  IN  A MIN. 

11FE 

PSPI 

EQU 

> 11FE 

PREVIOUS  SPINDLE  CHECKER 

11E0 

AT 

EQU 

> 11E0 

TOTAL  ALPHA  TIME 

11E2 

BT 

EQU 

> 11E2 

TOTAL  BETA  TIME 

11E4 

DT 

EQU 

> 11E4 

TOTAL  DELTA  TIME 

11E6 

DN 

EQU 

> 11E6 

TOTAL  # OF  DELTA  OCCURANCE 

11E8 

SN 

EQU 

> 11E8 

TOTAL  # OF  SPINDLE 

11EA 

TT 

EQU 

> 11EA 

TOTAL  THETA  TIME 

11D2 

ATM 

EQU 

> 11D2 

ALPHA  TIMER 

11D4 

BTM 

EQU 

> 11D4 

BETA  TIMER 

11D6 

DTM 

EQU 

> 11D6 

DELTA  TIMER 

11D8 

DTN 

EQU 

> 11D8 

# OF  DELTA  WAVE 

11DA 

STN 

EQU 

> 11DA 

# OF  SPINDLE 

11DC 

TTM 

EQU 

> 11DC 

THETA  TIMER 

11DE 

FTM 

EQU 

> 11DE 

ARTIFACT  TIMER 

DDDD 

DACD 

A 

EQU 

>DDDD 

DUMMY 

0000 'FFAA 

AORG 

>FFAA  ; 

INT  LEVEL  3 VECTOR  ADDRESS 

0460  1000 

B 

@>1000  ; 

TO  HANDLER 

FFAE  1000 

AORG 

>1000 

A 


* CHANNEL  MULTIPLEXING  AND  LOW  PASS  FILTERING 

* Rl:  CHANNEL  SELECTOR 

A 

04C1  CLR  Rl  ; CLEAR  CHANNEL  SELECTOR 
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PAGE:  2 

14:26 

09/14/84 

FILE:  HANDLER /l. SLEEP 

1 

1002:  04E0 

CFF8 

CLR 

0MUX 

GET  CHO 

2 

1006:  C801 

CFFA 

MOV 

R1 ,0CON 

INITIATE  A/D  CONVERSION 

3 

100A:  0581 

INC 

R1 

GET  CHI 

4 

100C:  C820 

11CF 

11F0 

MOV 

0CHOO , 0CHO 

; SHIFT  FOR  CHO  LPF 

5 

1012:  C020 

CFFC 

WAITO 

MOV 

0EOC.RO 

;GET  EOC  TEST  BIT 

6 

1016:  COOO 

MOV 

RO  ,R0 

;END  OF  CONVERSION? 

7 

1018:  15FC 

JGT 

WAITO 

8 

101A:  C820 

CFFE 

11CF 

MOV 

0ADC.0CHOO 

;PUT  VALUE  AT  CHOO 

9 

1020:  C801 

CFF8 

MOV 

R1,0MUX 

SELECT  CHI 

10 

1024:  C801 

CFFA 

MOV 

R1 ,0CON 

INITIATE  A/D  CONVERSION 

11 

1028:  0581 

INC 

R1 

GET  CH2 

12 

102A:  A820 

11CF 

11F0 

A 

0CHOO ,0CHO 

;LPF  OUT  OF  CHO 

13 

1030:  C820 

11CD 

11F2 

MOV 

0CH11 , 0CH1 

.•SHIFT  FOR  CHI  LPF 

14 

1036:  C020 

CFFC 

WAIT1 

MOV 

0EOC ,R0 

;GET  EOC  TEST  BIT 

15 

103A:  COOO 

MOV 

RO  ,R0 

;END  OF  CONVERSION? 

16 

103C : 15FC 

11CD 

JGT 

WAIT1 

;PUT  VALUE  AT  CH11 

17 

103E:  C820 

Ltt  £ 

MOV 

0ADC .0CH11 

18 

1044:  C801 

CFF8 

MOV 

R1 ,0MUX 

SELECT  CH2 

19 

1048:  C801 

CFFA 

MOV 

R1 ,0CON 

INITIATE  A/D  CONVERSION 

20 

104C : 0581 

INC 

R1 

GET  CH3 

21 

104E:  A820 

11CD 

11F2 

A 

0CH11 ,0CH1 

; OUTPUT  OF  CHI 

22 

1054:  C820 

11CB 

11F4 

MOV 

0CH22 ,0CH2 

; SHIFT  FOR  CH2  LPF 

23 

105A:  C020 

CFFC 

WAIT2 

MOV 

0EOC.RO 

;GET  EOC  TEST  BIT 

24 

105E:  COOO 

MOV 

RO.RO 

;END  OF  CONVERSION? 

25 

1060:  15FC 

JGT 

WAIT2 

26 

1062:  C820 

CFFE 

11CB 

MOV 

0ADC.0CH22 

.-PUT  VALUE  ATCH22 

27 

1068:  C801 

CFF8 

MOV 

R1.0MUX 

SELECT  CH3 

28 

106C : C801 

CFFA 

MOV 

R1.0CON 

INITIATE  A/D  CONVERSION 

29 

1070:  AQ20 

11CB 

11F4 

A 

0CH22.0CH2 

.•OUTPUT  OF  CH2 

30 

1076:  C820 

11C9 

11F6 

MOV 

0CH33.0CH3 

.•SHIFT  FOR  CH3 

31 

107C : C020 

CFFC 

WAIT3 

MOV 

0EOC.RO 

GET  EOC  TEST  BIT 

32 

1080:  COOO 

MOV 

RO.RO 

END  OF  CONVERSION 

33 

1082:  15FC 

JGT 

WAIT3 

34 

1084:  C820 

CFFE 

11C9 

MOV 

0ADC.0CH33 

.-PUT  VALUE  AT  CH33 

35 

108A:  A820 

11C9 

11F6 

A 

0CH33.0CH3 

.•OUTPUT  OF  CH3  LPF 

36 

* 

37 

* MINUTE  COUNT  AND  DATA  AQUISITION 

38 

* 

R2 : MINUTE  COUNTER 

39 

A 

R3 : STORAGE  POINTER 

40 

A 

R4 : ALPHA 

TIMER 

41 

A 

R5 : BETA  TIMER 

42 

A 

R6 : DELTA 

TIMER 

43 

A 

R7 : PREVIOUS  DELTA  CHECKER 

44 

A 

R8 : QUARTER  SEC  DATA 

STARTING  ADDR. 

45 

A 

R9 : DATA 

MULTIPLEXING  REG. 

46 

A 

RIO:  QUARTET  SEC.  TIMER 

47 

A 

48 

1090:  0582 

INC 

R2 

EVERY  SAMPLING  INTERVAL 

49 

1092:  058A 

INC 

RIO 

QUARTER  SEC  TIMER 

50 

1094:  028A 

0078 

Cl 

RIO. 120 

QUARTER  SEC  ELAPSED? 

51 

1098:  1143 

JLT 

MDATA 

KEEP  CHECKING 

52 

109A:  04C9 

CLR 

R9 

DATA  MUX  REG  CLEAR 

53 

109C : 8282 

C 

R2.R10 

FIRST  QUARTER  SEC 

54 

109E:  1502 

JGT 

NOMARK 

IN  MINUTE  INTERVAL 

55 

10A0 : 0269 

1000 

ORI 

R9 , > 1000 

TIME  MARK  SET 

56 

10A4 : 04CA 

NOMARh 

CLR 

RIO 

RESET  QUARTER  SEC  TIMER 

PDOS  ASM  R2.4 
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14:27 

09/14/84 

FILE:  HANDLER /l, SLEEP 

1 

2 

10A6 : C820 

F088 

F088 

* 

MOV 

0AEX,0AEX 

;TEST  ALPHA 

3 

10AC : 1302 

JEQ 

QBX 

; IF  NOT,  TEST  BETA 

4 

10AE : 0269 

0800 

ORI 

R9 , > 0800 

;SET  12TH  BIT 

5 

10B2 : C820 

F288 

F288 

QBX 

MOV 

0BEX,0BEX 

;TEST  BETA 

6 

10B8 : 1302 

JEQ 

QDX 

; IF  NOT,  TEST  DELTA 

7 

10BA:  0269 

0400 

ORI 

R9 , > 0400 

;SET  11TH  BIT 

8 

10BE:  C820 

F488 

F488 

QDX 

MOV 

0DEX , 0DEX 

;TEST  DELTA 

9 

10C4 : 1302 

JEQ 

QSX 

; IF  NOT,  TEST  SIGMA 

10 

10C6 : 0269 

0200 

ORI 

R9 , >0200 

;SET  10TH  BIT 

11 

10CA:  C820 

F688 

F688 

QSX 

MOV 

0SEX,0SEX 

;TEST  SIGMA 

12 

10D0 : 1302 

JEQ 

QTX 

; IF  NOT,  TEST  THETA 

13 

1002 : 0269 

0100 

ORI 

R9 , >0100 

;SET  9TH  BIT 

14 

10D6 : C820 

F888 

F888 

QTX 

MOV 

0TEX,0TEX 

;TEST  THETA 

15 

10DC : 1302 

JEQ 

QFX 

; IF  NOT,  TEST  ARTIFACT 

16 

10DE:  0269 

0080 

ORI 

R9 , >0080 

;SET  8TH  BIT 

17 

10E2:  C820 

FA88 

FA88 

QFX 

MOV 

0FEX,0FEX 

;TEST  ARTIFACT 

18 

10E8:  1302 

JEQ 

QOX 

; IF  NOT,  TEST  OREM 

19 

10EA:  0269 

0040 

ORI 

R9 , >0040 

;SET  6TH  BIT 

20 

10EE:  C820 

15B8 

15B8 

QOX 

MOV 

0OEX,0OEX 

;TEST  OREM 

21 

10F4 : 1302 

JEQ 

QPX 

; IF  NOT,  TEST  PREM 

22 

10F6:  0269 

0008 

ORI 

R9,>0008 

;SET  4TH  BIT 

23 

10FA:  C820 

17B8 

17B8 

QPX 

MOV 

0PEX,0PEX 

;TEST  PREM 

24 

1100:  1302 

JEQ 

QOX 

; IF  NOT,  TEST  QREM 

25 

1102:  0269 

0004 

ORI 

R9 , >0004 

;SET  3RD  BIT 

26 

1106:  C820 

19B8 

19B8 

QQX 

MOV 

0QEX,0QEX 

;TEST  QREM 

27 

HOC:  1302 

JEQ 

QRX 

; IF  NOT,  TEST  REM 

28 

110E:  0269 

0002 

ORI 

R9 , >0002 

;SET  2ND  BIT 

29 

1112:  C820 

1BB8 

1BB8 

QRX 

MOV 

0REX,0REX 

;TEST  REM 

30 

1118:  1302 

JEQ 

DMUX 

; IF  NOT,  PROCESS  MIN  DATA 

31 

111A:  0269 

0001 

ORI 

R9 , > 0001 

;SET  1ST  BIT 

32 

111E:  CE09 

DMUX 

MOV 

R9,*R8+ 

STORE  MULTIPLEXED  DATA 

33 

1120:  020C 

0120 

MDATA 

LI 

R12 , > 120 

;GET  TMS  9901  PIO  CRU  ADDR 

34 

1124:  0282 

03C0 

Cl 

R2 , 960 

2 SEC  INTERVAL 

35 

1128:  1502 

JGT 

PULSE 

36 

112A:  1D06 

SBO 

6 

OUTPUT  1 AT  PIO  P6 ( P4-12 ) 

37 

112C : 1001 

JMP 

PULSE1 

38 

112E:  1E06 

PULSE 

SBZ 

6 

OUTPUT  0 AT  PIO  P6 ( P4-12 ) 

39 

1130:  0282 

7080 

PULSE1 

Cl 

R2 ,480*60 

1 MINUTE  ELAPSED? 

40 

1134:  111F 

JLT 

COUNT 

KEEP  COUNTING 

41 

1136:  05A0 

11F8 

INC 

0CLOCK 

CLOCK 

42 

113A:  CCEO 

11F8 

MOV 

0CLOCK , *R3+  ; STORE  CLOCK  VALUE 

43 

113E:  CCC4 

MOV 

R4,*R3+ 

STORE  ALPHA  TIME 

44 

1140:  CCC5 

MOV 

R5 , *R3+ 

STORE  BETA  TIME 

45 

1142:  CCC6 

MOV 

R6 , *R3+ 

STORE  DELTA  TIME 

46 

1144:  CCEO 

11FA 

MOV 

0NDEL , *R3+ 

; STORE  # OF  DELTA  OCCURANCE 

47 

1148:  CCEO 

1 1FC 

MOV 

0NSPI,*R3+ 

; STORE  # OF  SIGMA  SPINDLE 

48 

114C : A804 

11E0 

A 

R4 , 0AT 

TOTAL  ALPHA  TIME 

49 

1150:  A805 

11E2 

A 

R5 ,0BT 

TOTAL  BETA  TIME 

50 

1154:  A806 

11E4 

A 

R6,0DT 

TOTAL  DELTA  TIME 

51 

1158:  A820 

11FA 

11E6 

A 

0NDEL , 0DN 

TOTAL  # OF  DELTA 

52 

115E:  A820 

11FC 

11E8 

A 

0NSPI.0SN 

TOTAL  # OF  SPINDLE 

53 

1164:  04C2 

CLR 

R2 

RESET  MINUTE  CONUTER 

54 

1166:  04C4 

CLR 

R4 

RESET  ALPHA  TIMER 

55 

1168:  04C5 

CLR 

R5 

RESET  BETA  TIMER 

56 

1 16A:  04C6 

CLR 

R6 

RESET  DELTA  TIMER 

PAGE:  4 

1 116C : 

2 1170: 

3 1174: 

4 1178: 

5 117A: 

6 117C : 

7 117E: 

8 1182: 

9 1184: 

10  1186: 
11  1188: 

12  118C : 

13  118E: 

14  1190: 

15  1192: 

16  1194: 

17  1196: 

18  119A: 

19  119C : 

20  11A0 : 

21  11A2: 

22  11A4 : 

23  11AA: 

24  11AC : 

25  11B0 : 

26 

27  11B4 : 

28  11B8 : 

29  11BA: 

30  11BC: 

31  11C0 : 

32  11C2: 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43  F020: 

44  F024 : 

45  F026 : 

46  F028 : 

47  F02A: 

48  F02C : 

49  F02E: 

50  F030 : 

51  F032 : 

52  F034 : 

53  F038 : 

54  F03E : 

55 

56  F042 : 
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14:27 

09/14/84 

FILE:  HANDLER /l, SLEEP 

04E0 

UFA 

CLR 

0NDEL 

RESET  DELTA  # COUNTER 

04E0 

11FC 

CLR 

0NSPI 

RESET  SIGMA  # COUNTER 

C060 

F088 

COUNT 

MOV 

0AEX,R1 

GET  ALPHA  EXIST 

C041 

MOV 

R1.R1 

IS  THERE  ALPHA? 

1301 

JEQ 

BX 

IF  NO,  TEST  BETA 

0584 

INC 

R4 

C060 

F288 

BX 

MOV 

0BEX,R1 

GET  BETA  EXIST 

C041 

MOV 

R1,R1 

IS  THERE  BETA? 

1301 

JEQ 

DX 

IF  NO,  TEST  DELTA 

0585 

INC 

R5 

BETA  TIME 

C060 

CO 

CO 

DX 

MOV 

0DEX,R1 

GET  DELTA  EXIST 

C041 

MOV 

R1  ,R1 

IS  THERE  DELTA? 

1305 

JEQ 

NODEX 

IF  NOT 

0586 

INC 

R6 

INCREASE  DELTA  TIME 

C1C7 

MOV 

R7  ,R7 

DELTA  EXISTED  PREVIOUSLY? 

1603 

JNE 

SX 

IF  SO,  TEST  SIGMA 

05A0 

UFA 

INC 

0NDEL 

OR  INCREASE  DELTA  COUNTER 

C1C1 

NODEX 

MOV 

R1  ,R7 

RECORD  EXISTENCE  OF  DELTA 

C060 

F688 

SX 

MOV 

0SEX,R1 

GET  SIGMA  EXIST 

C041 

MOV 

R1,R1 

IS  THERE  A SIGMA? 

1306 

JEQ 

NOSEX 

IF  NOT 

C820 

11FE 

11FE 

MOV 

0PSPI , 0PSPI  ; SIGMA  EXISTED  PREVIOUSLT 

1604 

JNE 

GO 

IF  PREVIOUSLY  EXISTED 

05A0 

11FC 

INC 

0NSPI 

RECORD  SIGMA  OCCURANCE 

C801 

11FE 

NOSEX 

MOV 

R1,QPSPI 

RECORD  SIGMA  EXISTENCE 

* RESET  INTERRUPT  MECHANISM 


020C 

0100 

GO 

LI 

R12 , > 100  ;PSI  CRU  BASE  ADDR.  SET 

1E00 

SBZ 

0 ;SET  9901  IN  INTERRUPT  MODE 

1D03 

SBO 

3 ; ENABLE  INT3  IN  9901 

0300 

0003 

LIMI 

3 ; ENABLE  INT3  IN  9900 

05CE 

INCT 

R14  ; SUPPORT  ' JMP  $'  INSTRUCTION 

0380 

RTWP 

LINK 

ALPHA/ 1 

IDT 

'ALPHA' 

* 

ALL  THE 

FILTERS  HAVE  DUMMY  OUTPUT  PORT  0DDDD 

A 

0CFFO , 

0CFF2  ARE  USED  AS  OUTPUT  PORTS 

*ADC 

EQU 

>CFFE  ; SIGNAL  INPUT  PORT  ADDR. 

F082 

ALPHA 1 

EQU 

>F082  jALPHADET  R1  PRESENT  INPUT 

F084 

ALPHA2 

EQU 

>F084  ;ALPHADET  R2  PREVIOUS  INPUT 

F1B4 

ALPOUT 

EQU 

>F1B4  ;ALPD  OUTPUT,  RIO  OF  ALPD 

11C4 

F020 

AORG 

>F020 

C060 

F1B4 

MOV 

0ALPOUT,R1 

0811 

SRA 

R1 , 1 ;MAX  GAIN  5,  ADJUST  BY  1/4 

6181 

S 

R1.R6 

A206 

A 

R6.R8 

C287 

MOV 

R7,R10 

081A 

SRA 

RIO , 1 

A20A 

A 

RIO ,R8 

A248 

A 

R8,R9 

0829 

SRA 

R9 , 2 

C809 

DDDD 

MOV 

R9.0DACD 

C820 

F082 

F084 

MOV 

0ALPHA1 , 0ALPHA2  ; STORE  PREVIOUS  VALUE 

C809 

F082 

MOV 

R9,0ALPHA1  ; STORE  PRESENT  VALUE 

C248 

MOV 

R8.R9 

121 


PAGE:  5 


16 

17 

ia 

19 

20 
21 
22 

23 

24 


35 

36 

37  F0A8 : 

38  FOAA: 

39  FOAC : 

40  FOAE: 

41  FOBO: 

42 

43 

44  F0B2 : 

45  F0B8 : 

46  FOBE: 

47  F0C4 : 

48  FOCA: 

49  FODO: 

50 

51 

52  F0D4 : 

53  F0D8 : 

54 

55 

56 


C041 

135D 

115C 

C082 

155A 
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FILE:  ALPHA /l, SLEEP 


1 

F044 : 

C 207 

MOV 

R7,R8 

2 

F046 : 

C1C6 

MOV 

R6,R7 

3 

F048 : 

C185 

MOV 

R5,R6 

4 

F04A: 

C144 

MOV 

R4,R5 

5 

F04C : 

C103 

MOV 

R3,R4 

6 

F04E: 

C0C2 

MOV 

R2,R3 

7 

F050 : 

C081 

MOV 

R1.R2 

8 

F052 : 

0380 

RTWP 

9 

LINK 

ALPHADET/ 1 

10 

IDT 

'ALPHADET' 

11 

F054  FOAO 

AORG 

>F0A0 

12 

A 

13 

*EACH  REGISTERS  DENOTE  FOR 

14 

* PRESENT 

INPUT:  R1 

15 

* PREVIOUS 

INPUT:  R2 

* 

* 

A 

A 

A 

A 

A 


EXIST:  R4 


HISTO: 
MESH  : 
TIMER: 
PEAK  : 
DTIMER 


R5 

R6 

R7 

R8 

R9 


: DELAY  TIMER  FOR  DISPLAY 


’“'PARALLEL  I/O  PORT  P0(P4-20)  FOR  DETECTOR  OUTPUT 
’“'AVERAGE  WAVE  TIME  6 CONSECUTIVE  WAVES  IS  USED 


25 

’“'FOR 

DETECTOR 

CRITERION 

26 

F070 

A1 

EQU 

>F070 

27 

F072 

A2 

EQU 

>F072 

28 

F074 

A3 

EQU 

>F074 

29 

F076 

A4 

EQU 

>F076 

30 

F078 

A5 

EQU 

>F078 

31 

F07A 

A6 

EQU 

>F07A 

32  FOAO: 

OOFF 

ONE 

DATA 

>FF 

33  F0A2 : 

020C  0120 

LI 

R12 , >120 

34  F0A6 : 

0609 

DEC 

R9 

ZERO  CROSSING  CHECK 
MOV  R1,R1 

JEQ  ANOCRO 

JLT  ANOCRO 

MOV  R2 ,R2 

JGT  ANOCRO 


MAGNIFIED  FOR  DAC  OUT 
GET  TMS  9901  PIO  CRU  ADDR. 
DTIMER  ADJUSTMENT 


WAVE (NT) >0? 

NO  CROSSING 

NO  CROSSING,  NEGATIVE  GOING 
WAVE(  (N-l)TXO? 

IF  NOT,  NO  CROSSING 


C820  F078  F07A 
C820  F076  F078 
C820  F074  F076 
C820  F072  F074 
C820  F070  F072 
C807  F070 


0288  0030 
110A 


AVERAGE  WAVE  DURAION  CHECK 


MOV 

MOV 

MOV 

MOV 

MOV 

MOV 


0A5 , @A6 
0A4 , @A5 
0A3 , 0A4 
0A2 ,0A3 
0A1 ,0A2 
R7,0A1 


FODA:  C145 


MINIMUM  AMPLITUDE  CHECK 

Cl  R8 , 48  ; LARGER  THEN  MIN  INPUT? 

JLT  ASENDO  ;IF  NOT,  NO  WAVE 

HISTORY  CHECK:  ZERO  CROSSING  AFTER  OVER  DURATION 
MOV  R5.R5  ; 1 MEANS  HAVING  HISTORY 
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1 

FODC : 1607 

JNE 

ASENDH 

; IF  IT  HAS  HISTORY,  SEND  0 

2 

* 

X 

;MESH  AND  SET  HISTORY  0 

4 

* WAVE  LENGTH (DURATION) 

CHECK 

5 

FODE:  0287  0009 

Cl 

R7,9 

LARGER  THEN  MIN  INTERVAL? 

6 

F0E2 : 1A05 

JL 

ASENDO 

IF  NO  WAVE,  SEND  0 TO  MESH 

7 

F0E4 : 0287  OOOF 

Cl 

R7 , 15 

SMALLER  THEN  MAX  INTERVAL? 

8 

F0E8 : 1B02 

JH 

ASENDO 

IF  NOT,  SEND  0 TO  MESH 

9 

FOEA:  1008 

JMP 

ASEND1 

IF  SO,  SEND  1 TO  MESH 

10 

A 

11 

* HISTORY  VALUE  RESET 

12 

FOEC : 04C5 

ASENDH 

CLR 

R5 

; SET  HISTO  TO  0 

13 

* 

14 

* SEND  0 TO 

MESH  THERE  IS  NO  WAVE 

15 

FOEE:  0286  0003 

ASENDO 

Cl 

R6 , 3 

MESH  > 3? 

16 

F0F2 : 1202 

JLE 

AMESHO 

MESH<4?,  THEN  SET  MESH  TO 

17 

F0F4 : 0606 

DEC 

R6 

MESH) 3?,  THEN  MESH:=MESH-1 

18 

F0F6 : 1006 

JMP 

AWAVE 

19 

F0F8 : 04C6 

AMESHO 

CLR 

R6 

SET  MESH  TO  0 

20 

FOFA:  100E 

JMP 

AOUTO 

21 

* 

22 

* SEND  1 TO 

MESH  IF  THERE  IS  A SATISFIED  WAVE 

23 

FOFC : 0286  0006 

ASEND1 

Cl 

R6 ,6 

REACHES  MAXIMUM  # OF  WAVE? 

24 

F100:  1401 

JHE 

AWAVE 

IF  SO,  KEEP  MAX  VALUE 

25 

F102:  0586 

INC 

R6 

IF  NOT,  MESH : =MESH+1 

26 

A 

27 

* IS 

THERE  A 

WAVE  WITH  PREVIOUS  OUTPUT? 

28 

F104 : 0284  OOFF 

AWAVE 

Cl 

R4 , >FF 

WAS  THERE  A WAVE:EXIST=FF? 

29 

F108 : 1A04 

JL 

ANONEX 

IF  NOT 

30 

A 

31 

* WHEN  THERE 

WAS  A WAVE 

32 

F10A:  0286  0004 

Cl 

R6,4 

CHECK  HYSTERESIS  CONDITION 

33 

F10E:  1409 

JHE 

A0UT1 

IF  SATTISFIED,SET  EXIST  TO 

34 

F110 : 1003 

JMP 

AOUTO 

IF  NOT  SET  EXIST  TO  0 

35 

A 

36 

* WHEN  THERE 

WAS  NO  WAVE 

PREVIOUSLEY 

37 

F112 : 0286  0006 

ANONEX 

Cl 

R6 ,6 

MESH) =WAVEIN? 

38 

F116 : 1405 

JHE 

A0UT1 

IF  SATISFIED,  SET  EXIST  TO 

39 

F118 : C249 

AOUTO 

MOV 

R9.R9 

DTIMER  ON? 

40 

F11A:  151B 

JGT 

AS ETC 0 

THEN  CONTINUE  TO  DISPLAY 

41 

FI 1C:  04C4 

CLR 

R4 

SET  EXIST  TO  0 

42 

FI IE:  1E00 

SBZ 

0 

OUTPUT  0 TO  PO ( P4-20 ) 

43 

F120 : 1018 

JMP 

ASETCO 

44 

F122 : 04C3 

A0UT1 

CLR 

R3 

45 

F124 : AOEO  F070 

A 

0A1 ,R3 

46 

F128 : AOEO  F072 

A 

0A2 ,R3 

47 

F12C : AOEO  F074 

A 

0A3 ,R3 

48 

F130 : AOEO  F076 

A 

0A4 ,R3 

49 

F134 : AOEO  F078 

A 

0A5 , R3 

50 

F138 : AOEO  F07A 

A 

0A6 ,R3 

51 

F13C : 0283  0038 

Cl 

R3 , 56 

52 

F140 : 1AEB 

JL 

AOUTO 

53 

F142 : 0283  0054 

Cl 

R3 , 84 

54 

F146 : 1BE8 

JH 

AOUTO 

55 

F148 : C120  FOAO 

MOV 

0ONE , R4  ; 

SET  EXIST  TO  1 

56 

F14C : 0209  0038 

LI 

R9 , 56  ; 

SET  MIN  ALPFA  TIME 
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FILE:  ALPHADET/1, SLEEP 

1 

F150 : 

1D00 

SBO 

0 

OUTPUT  1 TO  PO ( P4-20 ) 

2 

F152 : 

C804 

DDDD 

ASETCO 

MOV 

R4 , GDACD 

OUTPUT  EXIST  TO  DAC2 

3 

F156 : 

C0C6 

MOV 

R6  ,R3 

GET  MESH  VALUE  AT  R3 

4 

F158: 

0A73 

SLA 

R3 , 7 

MESH* 128 

5 

F15A: 

C803 

DDDD 

MOV 

R3 , GDACD 

OUTPUT  OF  MESH*128 

6 

F15E: 

04C7 

CLR 

R7 

INITIALIZE  TIMER 

7 

F160 : 

04C8 

CLR 

R8 

SET  PEAK  0 

8 

q 

F162 : 

0460 

F184 

B 

GAFINE 

10 

* WHEN  THERE 

IS  NO  ZERO  CROSSING 

11 

F166 : 

0287 

OOOF 

ANOCRO 

Cl 

R7 , 15 

CHECK  WAVE  LENGTH  (DURATION 

12 

F16A: 

1B08 

JH 

AOVERT 

IF  SO 

13 

F16C : 

0587 

INC 

R7 

TIMER : =TIMER+1 

14 

* 

15 

* PEAK  VALUE 

CHANGE  WHEN 

NOT  IN  OVER  TIME 

16 

F16E: 

C281 

MOV 

R1.R10 

GET  INPUT  AT  RIO 

17 

F170 : 

074A 

ABS 

RIO 

GET  ABSOLUTE  VALUE 

18 

F172 : 

8288 

C 

R8 ,R10 

PEAK  VALUE  > ABSOLUTE  INPUT 

19 

F174 : 

1507 

JGT 

AFINE 

THEN  NO  CHANGE 

20 

F176 : 

C20A 

MOV 

RIO ,R8 

GET  NEW  PEAK  VALUE 

21 

F178 : 

0460 

F184 

B 

GAFINE 

22 

* 

23 

* WHEN  WAVE  LENGTH  EXCEEDS  MAX  TIME 

24 

F17C: 

C160 

FOAO 

AOVERT 

MOV 

G0NE,R5 

SET  HISTO  TO  1 

25 

F180 : 

0460 

FOEE 

B 

GASENDO 

26 

F184 : 

0380 

AFINE 

RTWP 

27 

LINK 

ALPD/1 

28 

IDT 

'ALPD' 

29 

F1B4 

ALPOUT 

EQU 

>F1B4 

RIO  OF  ALPFW 

30 

FI  86 

F1C0 

AORG 

>F1C0 

31 

F1C0 : 

F1A0 

F1C4 

ALPD 

DATA 

>F1A0 ,ALPD1 

32 

F1C4 : 

C060 

11F4 

ALPD1 

MOV 

GCH2 ,R1 

33 

F1C8 : 

0821 

SRA 

Rl,2 

ADJUST  INPUT 

34 

F1CA: 

A0C1 

A 

R1  ,R3 

35 

F1CC : 

A143 

A 

R3  ,R5 

36 

F1CE: 

A1C5 

A 

R5,R7 

37 

FIDO: 

A1C6 

A 

R6  ,R7 

38 

F1D2 : 

A207 

A 

R7,R8 

39 

F1D4 : 

0828 

SRA 

R8 , 2 

ADJUST  FILTER  GAIN 

40 

F1D6 : 

C808 

DDDD 

MOV 

R8 , GDACD 

DUMMY  OUTPUT 

41 

F1DA: 

C808 

F1B4 

MOV 

R8 ,GALPOUT 

42 

FIDE: 

C207 

MOV 

R7  ,R8 

43 

F1E0 : 

C1C6 

MOV 

R6  ,R7 

44 

F1E2 : 

C185 

MOV 

R5.R6 

45 

F1E4 : 

C144 

MOV 

R4  ,R5 

46 

F1E6 : 

C103 

MOV 

R3,R4 

47 

F1E8 : 

C0C2 

MOV 

R2.R3 

48 

F1EA: 

C081 

MOV 

R1,R2 

49 

F1EC : 

0380 

RTWP 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 
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FILE:  FRONT/ 1, SLEEP 

LINK 

FRONT/ 1 

IDT 

'FRONT' 

F000 

ALPHAW 

EQU 

>F000 

•ALPHA  FILTER  WORK  SPACE 

F020 

ALPHA 

EQU 

>F020 

F080 

ADETW 

EQU 

>F080 

FOAO 

ADET 

EQU 

>FOAO 

F200 

BETAW 

EQU 

>F200 

F220 

BETA 

EQU 

>F220 

F280 

BDETW 

EQU 

>F280 

F2A0 

BDET 

EQU 

>F2A0 

FF98 

HQTIME 

EQU 

>FF98 

QUARTER  SEC  TIMER 

FF9E 

HMTIME 

EQU 

>FF9E 

MINUTE  TIMER  IN  HANDLER 

F08C 

AMESH 

EQU 

>F08C 

ALFA  DETECTOR  MESH  R6 

F28C 

BMESH 

EQU 

>F28C 

BETA  DETECTOR  MESH 

F48C 

DMESH 

EQU 

>F48C 

DELTA  MESH 

F68C 

SMESH 

EQU 

>F68C 

SIGMA  MESH 

F88C 

TMESH 

EQU 

>F88C 

THETA  MESH 

FA8C 

FMESH 

EQU 

>FA8C 

ARTIFACT  MESH 

15BC 

OMESH 

EQU 

>15BC 

OREM  MESH 

17BC 

PMESH 

EQU 

> 17BC 

PREM  MESH 

19BC 

QMESH 

EQU 

> 19BC 

QREM  MESH 
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FILE:  FRONT/ 1, SLEEP 

1 

1BBC 

RMESH 

EQU 

> 1BBC 

;REM  MESH 

2 

F492 

DISPL 

EQU 

>F492 

; DELTA  DISPLAYER ,S  TIMER 

3 

11F8 

CLOCK 

EQU 

> 11F8 

; MINUTE  CLOCK 

4 

c 

118E 

MCLOCK 

£ 

EQU 

> 118E 

; MINUTE  COUNTER 

J 

6 

1A72 

1220 

AORG 

>1220 

7 

1220:  02E0 

1200 

LWPI 

>1200 

8 

1224:  04E0 

FF98 

CLR 

0HQTIME 

1/4  SEC  HANDLER  TIMER 

9 

1228:  04E0 

FF9E 

CLR 

0HMTIME 

HANDLER  MINUTE  TIMER 

10 

122C:  04E0 

F08C 

CLR 

0AMESH 

INITIALIZE  ALFA  MESH 

11 

1230:  04E0 

F28C 

CLR 

0BMESH 

INITIALIZE  BETA  MESH 

12 

1234:  04E0 

F48C 

CLR 

0DMESH 

INITIALIZE  DELTA  MESH 

13 

1238:  04E0 

F68C 

CLR 

0SMESH 

INITIALIZE  SIGMA  MESH 

14 

123C:  04E0 

F88C 

CLR 

0TMESH 

INITIALIZE  THETA  MESH 

15 

1240:  04E0 

FA8C 

CLR 

0FMESH 

ARTIFACT  MESH 

16 

1244:  04E0 

15BC 

CLR 

0OMESH 

17 

1248:  04E0 

17BC 

CLR 

0PMESH 

18 

124C:  04E0 

19BC 

CLR 

0QMESH 

19 

1250:  04E0 

1BBC 

CLR 

0RMESH 

20 

1254:  04E0 

F492 

CLR 

0DISPL 

INITIALIZE  DELTA  TIMER 

21 

1258:  04E0 

11D2 

CLR 

0ATM 

ALPHA  TIMER 

22 

125C:  04E0 

11D4 

CLR 

0BTM 

BETA  TIMER 

23 

1260:  04E0 

11D6 

CLR 

0DTM 

DELTA  TIMER 

24 

1264:  04E0 

11D8 

CLR 

0DTN 

# OF  DELTA 

25 

1268:  04E0 

11DA 

CLR 

0STN 

# OF  SPINDLE 

26 

126C:  04E0 

11E0 

CLR 

0AT 

TOTAL  ALPHA  TIME 

27 

1270:  04E0 

11E2 

CLR 

0BT 

TOTAL  BETA  TIME 

28 

1274:  04E0 

11E4 

CLR 

0DT 

TOTAL  DELTA  TIME 

29 

1278:  04E0 

11E6 

CLR 

0DN 

TOTAL  #0F  DELTA 

30 

127C:  04E0 

11E8 

CLR 

0SN 

TOTAL  # OF  SPINDLE 

31 

1280:  04E0 

11F8 

CLR 

©CLOCK 

INITIALIZE  MINUTE  CLOCK 

32 

1284:  04E0 

118E 

CLR 

0MCLOCK 

INITIALIZE  MINUTE  COUNTER 

33 

1288:  0201 

CFFO 

LI 

R1 , >CFFO 

GET  OUTPUT  PORT  ADDRESS 

34 

128C : 04E1 

0006 

CLR 

06(1) 

SET  GAIN  1 

35 

1290:  04E1 

0008 

CLR 

08(1) 

DISABLE  MUX  AUTOINC  MODE 

36 

1294:  0201 

0100 

LI 

R1 , > 100 

37 

1298:  C801 

FFA2 

MOV 

R1 ,0>FFA2 

HANDLER  CRU  ADDRESS  SET 

38 

129C : 0201 

2000 

LI 

R1 , >2000 

MULTIPLEXED  DATA 

39 

12A0:  C801 

FF9A 

MOV 

R1 , 0>FF9A 

R8  IN  HANDLER,  STARTING  ADDR 

40 

12A4:  0201 

8300 

LI 

R1 , >8300 

STORAGE  ADDRESS 

41 

12A8:  C801 

FF90 

MOV 

R1 ,0>FF9O 

R3  IN  HANDLER,  STARTING  ADDR 

42 

12AC : 0201 

FOOO 

LI 

Rl.ALPHAW 

GET  SUBROUTINE  ALPHA  WP 

43 

12B0:  0202 

F020 

LI 

R2, ALPHA 

GET  SUBROUTINE  ALPHA  PC 

44 

12B4 : 0203 

F080 

LI 

R3 ,ADETW 

GET  ADET  WORK  SPACE  ADDRESS 

45 

12B8:  0204 

FOAO 

LI 

R4.ADET 

GET  DETECTOR  ADDRESS 

46 

12BC:  0207 

F200 

LI 

R7,BETAW 

47 

12C0 : 0208 

F220 

LI 

R8 , BETA 

48 

12C4 : 0209 

F280 

LI 

R9 , BDETW 

49 

12C8 : 020A 

F2A0 

LI 

R10.BDET 

50 

12CC : 020C 

0100 

LI 

R12 , >100 

CRU  BASE  ADDRESS  TO  PS I 

51 

12D0 : 1E00 

SBZ 

0 

INABLE  INTERRUPT 

52 

1202 : 1D03 

SBO 

3 

SET  INT3 

53 

12D4 : 0300 

0003 

LIMI 

3 

SET  INT3  ON  CPU 

54 

12D8 : 0200 

00C3 

LI 

RO, >C3 

INITIALIZE  COUNTER 

55 

12DC : 33C0 

LDCR 

RO  , 15 

START  COUNTER 

56 

12DE : 10FF 

BACK 

JMP 

3 

WAIT  FOR  INTERRUPT  1ST 
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1 

12E0 : 0407 

BLWP 

2 

12E2:  0409 

BLWP 

3 

12E4 : 0420 

1770 

BLWP 

4 

12E8:  0420 

17D0 

BLWP 

5 

12EC:  0420 

1970 

BLWP 

6 

12F0 : 0420 

19D0 

BLWP 

7 

12F4 : 0420 

1B70 

BLWP 

8 

12F8 : 0420 

1BD0 

BLWP 

9 

12FC : 10FF 

JMP 

10 

12FE:  0420 

F1C0 

BLWP 

11 

1302:  0420 

F3C0 

BLWP 

12 

1306:  0420 

F5A0 

BLWP 

13 

130A:  0420 

F4A0 

BLWP 

14 

130E:  0420 

F620 

BLWP 

15 

1312:  0420 

F6A0 

BLWP 

16 

1316:  0420 

FA20 

BLWP 

17 

131A:  0420 

FAAO 

BLWP 

18 

131E:  10FF 

JMP 

19 

1320:  0401 

BLWP 

20 

1322:  0403 

BLWP 

21 

1324:  0407 

BLWP 

22 

1326:  0409 

BLWP 

23 

1328:  0420 

F820 

BLWP 

24 

132C : 0420 

F8A0 

BLWP 

25 

1330:  0420 

1570 

BLWP 

26 

1334:  0420 

15D0 

BLWP 

27 

1338:  10FF 

JMP 

28 

133A:  0420 

F1C0 

BLWP 

29 

133E:  0420 

F3C0 

BLWP 

30 

1342:  0420 

F620 

BLWP 

31 

1346:  0420 

F6A0 

BLWP 

32 

134A:  0420 

FA20 

BLWP 

33 

134E:  0420 

FAAO 

BLWP 

34 

1352:  10FF 

JMP 

35 

1354:  0407 

BLWP 

36 

1356:  0409 

BLWP 

37 

1358:  0420 

1770 

BLWP 

38 

135C : 0420 

17D0 

BLWP 

39 

1360:  0420 

1970 

BLWP 

40 

1364:  0420 

19D0 

BLWP 

41 

1368:  0420 

1B70 

BLWP 

42 

136C : 0420 

1BD0 

BLWP 

43 

1370:  10FF 

JMP 

44 

1372:  0420 

F1C0 

BLWP 

45 

1376:  0420 

F3C0 

BLWP 

46 

137A:  0420 

F620 

BLWP 

47 

137E:  0420 

F6A0 

BLWP 

48 

1382:  0420 

FA20 

BLWP 

49 

1386:  0420 

FAAO 

BLWP 

50 

138A:  1001 

JMP 

51 

138C:  10A8 

POST 

JMP 

52 

138E : 10FF 

BRIDGE 

JMP 

53 

1390:  0401 

BLWP 

54 

1392:  0403 

BLWP 

55 

1394:  0407 

BLWP 

56 

1396:  0409 

BLWP 

PDOS  ASM  R2.4 
FILE:  FRONT/ 1, SLEEP 

7 

9 

0PLPFF 

0PREMD 

0QLPFF 

0QREMD 

0RLPFF 

0REMD 

$ ; 2ND 

0ALFD 

0LPD 

0LOW 

0DELD 

0SIG 

0SIGD 

0ARTIF 

0ARTD 

$ ; 3RD 

1 
3 
7 
9 

0THETA 

0THED 

0OLPFF 

0OREMD 

$ ;4TH 

0ALPD 

0LPD 

0SIG 

0SIGD 

0 ART IF 

0ARTD 

$ ; 5TH 

7 

9 

0PLPFF 

0PREMD 

0QLPFF 

0QREMD 

0RLPFF 

0REMD 

$ ; 6TH 

0ALPD 

0LPD 

0SIG 

0SIGD 

0 ART IF 

0ARTD 

BRIDGE 

BACK 

$ ; 7TH 

1 
3 
7 
9 
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1 

1398:  0420 

F820 

BLWP 

0THETA 

2 

139C : 0420 

F8A0 

BLWP 

0THED 

3 

13A0 : 0420 

1570 

BLWP 

0OLPFF 

4 

13A4 : 0420 

15D0 

BLWP 

0OREMD 

5 

13A8 : 10FF 

JMP 

$ 

6 

13AA:  0420 

F1C0 

BLWP 

0ALPD 

7 

13AE:  0420 

F3C0 

BLWP 

0LPD 

8 

13B2:  0420 

F5A0 

BLWP 

0LOW 

9 

13B6:  0420 

F4A0 

BLWP 

0DELD 

10 

13BA:  0420 

F620 

BLWP 

0SIG 

11 

13BE:  0420 

F6A0 

BLWP 

0SIGD 

12 

13C2 : 0420 

FA20 

BLWP 

0ARTIF 

13 

13C6 : 0420 

FAAO 

BLWP 

0ARTD 

14 

13CA:  10FF 

JMP 

$ 

15 

13CC : 0407 

BLWP 

7 

16 

13CE:  0409 

BLWP 

9 

17 

13D0 : 0420 

1770 

BLWP 

0PLPFF 

18 

13D4 : 0420 

17D0 

BLWP 

0PREMD 

19 

13D8 : 0420 

1970 

BLWP 

0QLPFF 

20 

13DC : 0420 

19D0 

BLWP 

0QREMD 

21 

13E0:  0420 

1B70 

BLWP 

0RLPFF 

22 

13E4:  0420 

1BD0 

BLWP 

0REMD 

23 

13E8:  10FF 

JMP 

$ 

24 

13EA:  0420 

F1C0 

BLWP 

0ALPD 

25 

13EE:  0420 

F3C0 

BLWP 

0LPD 

26 

13F2 : 0420 

F620 

BLWP 

0SIG 

27 

13F6 : 0420 

F6A0 

BLWP 

0SIGD 

28 

13FA:  0420 

FA20 

BLWP 

0ARTIF 

29 

13FE:  0420 

FAAO 

BLWP 

0ARTD 

30 

1402:  10FF 

JMP 

$ 

31 

1404:  0401 

BLWP 

1 

32 

1406:  0403 

BLWP 

3 

33 

1408:  0407 

BLWP 

7 

34 

140A:  0409 

BLWP 

9 

35 

140C : 0420 

F820 

BLWP 

0THETA 

36 

1410:  0420 

F8A0 

BLWP 

0THED 

37 

1414:  0420 

1570 

BLWP 

0OLPFF 

38 

1418:  0420 

15D0 

BLWP 

0OREMD 

39 

141C : 10FF 

JMP 

$ 

40 

141E:  0420 

F1C0 

BLWP 

0ALPD 

41 

1422:  0420 

F3C0 

BLWP 

0LPD 

42 

1426:  0420 

F620 

BLWP 

0SIG 

43 

142A : 0420 

F6A0 

BLWP 

0SIGD 

44 

142E:  0420 

FA20 

BLWP 

0ARTIF 

45 

1432:  0420 

FAAO 

BLWP 

0ARTD 

46 

1436:  10AA 

JMP 

POST 

47 

1438: 

END 

PDOS  ASM  R2.4 
FILE:  FRONT/ 1, SLEEP 


; 8TH 


; 9TH 


; 10TH 


; 11TH 


; 12TH 


; JUMP  LIMIT 


APPENDIX  2 

LINEAR  FIR  FILTER  DESIGN  AND  IMPLEMENTATION 
An  alpha  filter  illustrates  the  filter  design  algorithm 
used  in  the  present  system.  Other  filter  functions  are  also 
given  with  the  sampling  rates. 

LiDSJX  £IB  jLilfcsi:  -P.e_si.gn 

1)  Calculate  the  number  of  zeros  using  the  following  equa- 
tion. 

fs  = N * B 

where  fs  : sampling  frequency, 

N : required  number  of  zeros  of  the 
filter , 

B : required  passband. 

2)  Place  the  zeros  on  the  unit  circle  of  the  Z-plane. 

3)  Reduce  the  stopband  gain  by  placing  zeros  in  the  stop- 
band. 

4)  Analyze  the  filter  response  using  DINAP  [Ba78]. 

5)  If  necessary,  add  more  zeros  or  adjust  the  zero  positions 
to  make  the  filter  transfer  function  simple  and  then  go  to 
step  4. 
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Linear  Filter  Implementation 

Ex)  Alpha  Filter 

Passband  : 5.0  - 24.0  Hz 

Zeros  : 0,  2/5,  4/5,  6/5,  8/5  , 2/3,  4/3,  1 pi 

Transfer  Function  : 

-5  -1  -2 

( 1 + Z ) (-1  + Z ) ( 1 + Z/2  + Z ) 

Sampling  Rate  : 160  Hz 


Working  registers  are  used  as  shift  registers. 
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FREQUENCY  ANALYSIS  NODE  NUMBER  15 
MAONI ruoE 

8 32702  *• 


7.  49442 


0 83360 
0.  00100.  * 


144  OOOOO 


-100  OOOOO 


144  OOOOO 


-180  OOOOO 

O lOC-Ol 


O 32C  02 


0 OOE  02 


0 13E  03 


O I6E  03 
0 1 4E  03 
f-REOUENCY  < Ml  Z > 


Alpha  filter  frequency  character i st i c 


APPENDIX  3 

TOKEN  PROCESSOR  FILES 


The  following  are  the  parts  of  the  main  processor  pro- 
gram. Each  name  of  the  file  represents  the  functional  name 
of  the  token  processor.  The  'pdpro.*'  is  the  program  for 
the  translator  of  the  main  processor  and  the  others  are  for 
the  channel  descriptor.  The  token  translator,  alpha  (al- 
pha.*) and  REM  (rslow.*)  descriptor  are  given  in  order. 
'Makefile'  is  a kind  of  command  file  which  executes  a 
series  of  shell  commands  in  the  UNIX  operarting  system, 
'lexspecs',  'gram.g',  ' semaction . c ' are  the  input  files  for 
the  lexical  analyzer  and  parser  generator.  These  files  are 
used  for  the  signal  language  compiler  specification.  Most  of 
the  signal  understanding  rules  (channel  descriptor)  are  in 
these  files. 
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***  pdpro . Makefile  *** 


FILES=main.o  lex.yy.o  parser. o 

firstpass:  $(FILES) 

cc  -o  firstpass  $(FILES)  -11 
main.o:  raain.c 

cc  -c  main.c 
lex.yy.o:  lex.yy.c 

cc  -c  lex.yy.c 
lex.yy.c:  lexspecs 

lex  lexspecs 
parser. o:  parser. c 

cc  -c  parser. c 


***  pdpro. lexspecs  *** 


%{ 

^include  <stdio.h> 
int  linenurarO; 
int  datanura=0; 
extern  int  listing; 
%} 


HEXNUM 

[ A-Fa-f0-9 ] 

/* 

hex 

notation  * 

/ 

ALPHA 

[8-9A-Fa-f] 

/* 

1st 

bit 

”8" 

of 

2nd 

nibble 

*/ 

BETA 

[ 4-7  C-Fc-f ] 

/* 

2nd 

bit 

mpi 

of 

2nd 

nibble 

*/ 

DELTA 

[2367ABabEFeF] 

/* 

3rd 

bit 

»t  2 M 

of 

2nd 

nibble 

*/ 

SIGMA 

[ 1 3579BbDdFf ] 

/* 

4th 

bit 

n i ii 

of 

2nd 

nibble 

*/ 

THETA 

[8-9A-Fa-f] 

/* 

1 St 

bit 

li  8 n 

of 

3rd 

nibble 

*/ 

ARTIFACT 

[4-7  C-Fc-f] 

/* 

2nd 

bit 

ii  4 it 

of 

3rd 

nibble 

*/ 

REM 

[8-9A-Fa-f] 

/* 

1 st 

bit 

n 8 ii 

of 

4 th 

nibble 

*/ 

OREM 

[4-7C-Fc-f] 

/* 

2nd 

bit 

114  ii 

of 

4 th 

nibble 

»/ 

PREM 

[2367AaBbEeFf] 

/* 

3rd 

bit 

ii2*i 

of 

4 th 

nibble 

*/ 

QREM 

[ 1 3579BbDdFf ] 

/* 

4 th 

bit 

11 1 11 

of 

4 th 

nibble 

*/ 

%% 

A . *\n  {linenum  +=  1 ; 

if(listing)  printfC” [%05d]%s” , linenum, yy text) ; 
REJECT; 

} 

[ A-Fa-fO-9] (4 ,4 } ( r ) ; /#  address,  ignore  */ 

[01 ] [A-Fa-f0-9] [4-7C-Fc-f] [A-Fa-f0-9]  { 

/*  remove  artifact  */ 
yytexttl]  = 'O’; 
yytext [2]  = ’4  ' ; 
yytext[3]  = 'O'; 

REJECT; 

} 


[01 ] [2367ABabEFef] [A-Fa-fO-93 [A-Fa-fO-93  ( 


/*  alpha 

with  delta 

*/ 

switchCyy text[ 1 3 ) 

{ 

case  ' A * 
case  'a* 

yy text [ 1 3 

= '2 ' ; 

break ; 

case  'B' 
case  'b* 
case  'E1 

yy text[ 1 3 

= *3  * ; 

break ; 

case  ’e’ 
case  'F* 

yytext [ 1 3 

= '6  1 ; 

break; 

case  ' f * 

yy text [ 1 3 

II 

break ; 

default 

} 

REJECT; 

break ; 

} 

[01 3 [2367ABEFabef 3 [A-Fa-fO-93 [89ABCDEFabcdef 3 { 


/*  rem  with  delta  */ 
switch(yytext[33 ) 

{ 

case 

»8 

yytext[33  = 

•O’ ; 

break; 

case 

,9, 

yytext[33  = 

*1  •; 

break; 

case 

'A* 

case 

* a' 

yytext[33  = 

’ 2 ' ; 

break ; 

case 

•B' 

case 

* b* 

yytext[33  = 

’3'; 

break; 

case 

»E' 

case 

'e' 

yytext[33  = 

*6  1 ; 

break ; 

case 

»F' 

case 

i f i 

yytext[33  = 

•7*1 

break ; 

default 

} 

REJECT; 

} 

break ; 

[01 ] [A-Fa-fO-93 [A-Fa-fO-93 [89ABCDEFabcdef ] { 


/*  rem  with  slow  wave  */ 
switch(yy text [33 ) { 

case  * A' : 

case  ' a ’ : 
case  T B f : 

yytext[33  = '2'; 

break ; 

case  ’ b ’ : 
case  * E * : 

yytext[33  = *3'; 

break ; 

case  * e ' : 
case  * F ' : 

yytext[33  = ’6'; 

break ; 

case  ' f 1 : 
default  : 
REJECT; 

yytext[33  = '7’; 
break;  } 

break ; 

} 

[A-Fa-fO-93(4,4}  { 
datanurn++ ; 
return( 1 ) ; 

} 

/*  white  space, 
/*  white  space, 
/*  and  count  the 


/*  octal  number 


ii  ii 

\t 

\n 


ignore  */ 
ignore  */ 
number  of 


*/ 


columns  */ 


{ 
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printf( "Scanner  error  in  line  %d:"); 
printfC"  Illegal  Character linenura) ; 
if ( yy text[0] >=  ' ’) 

print f ( "5&c\n" , yy  text  [ 0] ) ; 


else 

} 


printfC "ascii  $03o0 ,yytext[0] ) ; 


***  pdpro. parser . c *** 


extern  char  yytext[]; 
extern  int  linenum; 
extern  int  datanura; 


parser ( ) 

{ 

int  tok; 


} 


whileC (tok=yylex() ) !=0)  { 

printf("Js  ",yytext); 
if  (datanum%8  ==  0) 
printfC "0) ; 


} 


*##  pdpro. main. c *** 


#include  <stdio.h> 
int  listing=0,  trace=0; 

FILE  *obj  f ; 

static  char  *objnarne="asra.  inp"  ; 

raainCargc , argv) 
int  argc; 
char  *argv [ ] ; 

{ 

int  i ; 

for(i=1;i<argc&&argv[i][0]==,-,; i++)  { 
switchC argv [ i] [ 1 ] ) { 
case  ’ 1 * : 

listing=1 ; 
break ; 
case  1 t ' : 

tracer  1 ; 
break ; 
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case  ’ o’ : 

objname=argv[++i] ; 
break ; 
case  • d ' : 

obj  name=0 ; 
break; 

} 

} 

if(i<argc)  { 

fclose ( stdin) ; 

if ( fopen (argvt i] , "r " ) ==NULL)  { 
printf("%s  error:  Cannot  open  JsO 
, argv [ 0] ,argv[i]) ; 
exit( 1 ) ; 

} 

} 

parser ( ) ; 
exit(0) ; 


#**  alpha. Makefile  *** 


f ilesrmain. o lex.yy.o  parser. o ptables.o  semaction.o 

alpha:  $(files) 

cc  -o  alpha  $(files)  -11 
lex.yy.o:  lex.yy.c 

cc  -c  lex.yy.c 
lex.yy.c:  lexspecs 

lex  lexspecs 
main.o:  raain.c 

cc  -c  raain.c 

parser. o:  parser. c ptables.h 
cc  -c  parser. c 

ptables.o:  ptables.c  ptables.h 
cc  -c  ptables.c 
semaction.o:  semaction.c 
cc  -c  semaction.c 
ptables.c:  gram.g 

slrgen  -u  gram.g  ptables.c 


**#  alpha. lexspecs  *** 


%{ 

//include  <stdio.h> 
int  linenum=0; 
int  datanum=0; 
int  alphatime=0; 


int  minute=0; 
extern  int  listing; 
%} 


- 136 


DIGIT 

[01] 

HEXNUM 

[0-9A-Fa-f] 

ALPHA 

[89A-Fa-f] 

NOALPHA 

[0-71 

BETA 

[4-7C-Fc-f ] 

DELTA 

[2367ABEFabef] 

SIGMA 

[13579BDFbdf] 

THETA 

[8-9A-Fa-f] 

ARTIFACT 

[4-7C-Fc-f] 

REM 

[8-9A-Fa-f] 

OREM 

[4-7C-Fc-f] 

PREM 

[2367ABEFabef ] 

QREM 

[ 1 3579BDFbdf ] 

alphadata 

[0] {ALPHA} {HEXNUM} {HEXNUM} 

notalpha 

[0] {NOALPHA} {HEXNUM} {HEXNUM} 

alphamin 

[ 1 ] {ALPHA} {HEXNUM} {HEXNUM} 

notalphamin 

%% 

[ 1 ] {NOALPHA} {HEXNUM} {HEXNUM} 

.*  \n  { 

linenura  +=  1 ; 
if ( listing) 

printf ( n [%05d]%s"f linenum , yy text ) ; 
REJECT; 

} 

{alphadata}  { 

alphatirae++ ; 
return ( 1 ) ; 

} 

{notalpha}  { 

return(2) ; 

} 


{alphamin}  { 

alphatime++ ; 
rainute++ ; 
return ( 3) ; 

} 

{ notalphamin}  { 


minute++ ; 
return ( 4 ) ; 

} 

" " ; /*  white  space,  ignore  #/ 

\t  ; /#  white  space,  ignore  */ 

\n  ; /*  and  count  the  number  of  columns 

{ 

printf ( "Scanner  error  in  line  ?d:H); 
printfC  Illegal  Character linenum) 
if (yytext[0]>= ’ •) 

printf("?c0,yytext[0]); 


*/ 
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el  se 
} 


printfC "ascii  ?03o0 ,yytext[0] ) ; 


***  alpha. gram. g *** 


^terminals 

alphadata  1 

notalpha  2 

alphamin  3 

notalphamin  4 

$endraarker 
0 

^grammar 
data  : 


datal 


data2 


null 


datal 

! 

//5 

• 

• 

data2  datal 

\ null 

alphadata 

#1 

! notalpha 

#2 

I alphamin 

#3 

i notalphamin 

//4 

***  alpha . parser . c *** 


//include  <stdio.h> 

//include  "ptables . h" 

extern  int  trace,  linenum; 
extern  char  yytext[]; 

//define  ERROR  0 
//define  SHIFT  1 
//define  REDUCE  2 
// define  ACCEPT  3 

static  gotof ( st , sym) 
register  int  st ; 
int  sym; 

{ 


register  int  ra,n; 
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for (mrgotop[sym-numterms] , 

n=gotop [ sym-numterms+1 ]-2;m<n;m+=2)  { 
if (gotov[m]==st) 

return(gotov[m+1 ] ) ; 

} 

return(gotov[n+1 ] ) ; 

} 

static  actionf ( st, sym) 
short  st; 

register  short  sym; 

{ 

register  int  m,n; 

for(mractionp[st],n=actionp[st+1]-2; 
m<=n;m+=2)  { 

if(actionv[m]==sym) 

return(actionv[m+1]); 

} 

return(actionv[n+1 ] ) ; 

} 

//define  SIZE^STATE  16000 

#def ine  SIZE^TRACE  16000 

int 
int 
int 
int 


static 
int  n ; 
{ 

} 

static 

{ 

} 

static 
int  n ; 
{ 


} 

static 

short 

{ 


st^stacktSIZE^STATE] ; 
tr^stack[SIZE_TRACE] ; 

stptr  =0;  /*state . stack  pointer*/ 

trptr  =0;  /*trace  stack  pointer*/ 


pspop(n) 


stptr  -=n; 


pstop( ) 

return(st_stack[stptr-1 ] ) ; 


pspushC  n) 


st^stackt stptr ] =n; 
stptr  += 1 ; 


shifty trace( token) 
token ; 


tr_stack[trptr ] =token; 
trptr  += 1 ; 
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} 


dump^trace ( ) ; 


static  r educe^trace (n , nonterm) 
short  n, nonterm; 

{ 

trptr  -=n; 

tr_stack[ trptr ] monterm; 

trptr  +=1  ; 

dump.  trace( ) ; 

} 

static  dump  traceO 

{ 

int  char^sum,  i; 

char^sum  =0; 

for  (i  =0;  i <trptr;  i++){ 

if  (char^sum  +strlen 
(snamesttr  stack[i]])  > =80 ) { 
printf ( ”0) ; 
char  sum  =10; 

} 

printf("  %s  ",  snames[tr^stack[i] ] ) ; 

char^sum  +=strlen(snames[tr^stack[i] ] ) +2; 

} 

printf ( "0) ; 

} 


%s") ; 

5&s0  , sn am es[  token  ] , 


static  tracelex ( token ) 
short(token) ; 

{ 

printf ("  NEW  TOKEN 

printf (»  INSTANCE 

yytext) ; 

} 

parserO  { 

register  int  token,  act; 


pspushC  0) ; 

token=tnmap[yylex()] ; 
while(l)  { 

act  = actionf(pstop() , token) ; 

switch(act%4)  { 

case  ACCEPT: 

return ( 0) ; 
case  ERROR: 

printf("#%d  syntax  error:  at"); 
printfC"  '%s * 0 , linenum , yy text ) ; 
return ( 1 ) ; 
case  SHIFT: 
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pspush(act/4) ; 

if  (trace  ==1)  shift^tr ace ( token ) ; 
token=tnmap[yylex() ] ; 
if  (trace  ==1)  tracelex ( token ) ; 
break; 

case  REDUCE: 

semaction( pr sera [ act/ 4 ] ) ; 
pspop(prlen[act/4]); 
ps push (got of (pstopO ,prlhs[act/4])) ; 
if  (trace  ==1) 

reduce_trace(prlen[act/4] , 
prlhs[act/4] ) ; 

break ; 

} 

} 

} 


***  alpha. ptables. c *#* 


extern  short  numsyms; 
extern  short  nuraprods; 
extern  short  nuraterras; 
extern  short  nuranterras; 
extern  short  nurastates; 
extern  short  tnraap[]; 
extern  short  prlhs[]; 
extern  short  prlen[]; 
extern  short  prsem[3; 
extern  char  * snames[]; 
extern  int  gotop[]; 
extern  short  gotov[]; 
extern  int  actionp[3; 
extern  short  actionvi]; 


***  alpha. main. c *** 


//include  <stdio.h> 
int  listingrO,  trace=0; 

FILE  *obj  f ; 

static  char  *objnarae="asra. inp" ; 

raain(argc , argv) 
int  argc; 
char  *argv [ ] ; 

{ 


int  i ; 
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} 


for ( i=1 ; i<argc&&argv [i][0]==T— *; i++)  { 
switch( argv [ i] [ 1 ] ) { 
case  ’ 1*  : 

listing=1 ; 
break ; 
case  ’ t ’ : 

trace= 1 ; 
break; 
case  ' o* : 

objname=argv[++i] ; 
break ; 
case  1 d T : 

obj  name=0 ; 
break ; 

} 

} 


if(i<argc)  { 

f close ( stdin) ; 

if ( fopen(argv[ i] , "r")==NULL)  { 

printf(”%s  error:  Cannot  open  5&sO); 
print f(argv[0] ,argv[i]) ; 
exit( 1 ) ; 

} 


} 


parser ( ) ; 
exit( 0) ; 


***  alpha. semaction. c *** 


#include  <stdio.h> 
int  nsec=0; 

extern  int  minute; 
extern  int  alphatirae; 

describe ( ) 

{ 

if  ( alphatime  >=  120  ) pr intf ( "H" ) ; 
else  if  ( alphatime  > 48  ) printf ( nM" ) ; 
else  if  ( alphatime  > 10)  printf ( "L" ) ; 
else  printf ( "N" ) ; 

if  ( minute  % 60  ==  0)  printfC’O); 

} 

semaction(n ) 
int  n ; 


{ 


switch(n)  { 
case  0: 

break; 

case  1 : 

if  ( nsec  < 4 ) { 

alphatime  +=  nsec; 
nsec  = 0; 

} 

break; 

case  2: 

if  ( nsec  < 4 ) nsec++; 
break ; 

case  3: 

if  ( nsec  <4  ) alphatime  +=  nsec; 
describe( ) ; 
alphatime  = 0; 
break ; 

case  4: 

describe ( ) ; 
alphatime  = 0; 
break ; 

case  5: 

printf ( "0) ; 
break ; 

default: 

printf ( "*-*-*-*  semantic  action"); 
printfC"  % d not  handled. 0 ,n) ; 
break ; 


} 

} 

***  rslow. Makefile  *** 


files=raain.o  lex.yy.o  parser. o ptables.o  semaction.o 

rslow:  $(files) 

cc  -o  rslow  $(files)  -11 
lex.yy.o:  lex.yy.c 

cc  -c  lex.yy.c 
lex.yy.c:  lexspecs 

lex  lexspecs 
main.o:  main.c 

cc  -c  main.c 

parser. o:  parser. c ptables.h 
cc  -c  parser. c 

ptables.o:  ptables.c  ptables.h 
cc  -c  ptables.c 
semaction.o:  semaction.c 
cc  -c  semaction.c 
ptables.c:  gram.g 

slrgen  -u  gram.g  ptables.c 
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***  rslow. lexspecs  *** 


%{ 

//include  <stdio.h> 
int  linenura=0; 
int  datanura=0; 
int  rslowtime=0; 
int  r slowrecord=0 ; 
int  del tarecord=0 ; 
int  rainute=0; 


extern  int  listing; 

%} 


DIGIT 

HEXNUM 

ALPHA 

BETA 

DELTA 

NODELTA 

SIGMA 

THETA 

ARTIFACT 

REM 

NOREM 

OREM 

PREM 

QREM 

rslowdata 
notr slow 
r slowmin 
notrslowmin 
delta 
notdelta 


[01] 

[0-9A-Fa-f] 

[89A-Fa-f] 

[ 4-7  C-Fc-f] 

[2367ABEFabef] 

[014589CDcd] 

[1357  9BDFbdf ] 

[8-9A-Fa-f] 

[ 4-7  C-Fc-f] 

[8-9A-Fa-f] 

[0-7] 

[4-7  C-Fc-f] 

[2367ABEFabef ] 

[ 1 3579BDFbdf ] 

[0] {HEXNUM} {HEXNUM} {REM} 

[0]  {HEXNUM} {HEXNUM} {NOREM} 

[ 1 ]  {HEXNUM} {HEXNUM} {HEXNUM} 

[ 1 ] {HEXNUM} {HEXNUM} {HEXNUM} 
[01 ]{ DELTA} {HEXNUM} {HEXNUM} 
[01  ] {NODELTA} {HEXNUM} {HEXNUM} 


%% 

* . *\n 


{delta } 


{notdelta } 


{ 

linenum  +=  1 ; 
if ( listing) 

printf ( " [?05d]%s" , linenum , yy text ) ; 
REJECT; 

} 

{ 

if  ((  rslowrecord  > 0)  &&  ( rslowtime  > 0)) 
rslowtime  - = 1 ; 
deltarecord  = 8; 

REJECT; 

} 

{ 

if  ( deltarecord  > 0)  deltarecord  -=1 ; 
REJECT; 


{ rslowdata } 


{notrslow} 


{rslowmin} 


{notrslowmin} 


ii  ti 

\t 
\n 

{ 

printfC "Scanner  error  in  line  $d:"); 
printfC "Illegal  Character:", linenum) ; 
if(yytext[0]>='  ') 

printf ( "?c0,yytext[0] ) ; 

else 

printf ( "ascii  $03o0 ,yy text[0] ) ; 

} 


*#*  rslow.gram.g  *** 


^terminals 

rslowdata  1 

notrslow  2 

rslowmin  3 

notrslowmin  4 

% endmarker 
0 

/^grammar 
data  : 

datal  #5 

! 

datal  : 

data2  datal 
i null 


rslowdata  #1 

inotrslow  #2 


} 

{ 

r slowtime++ ; 
rslowrecord  = 8; 
return( 1 ) ; 

} 

{ 

if  ( rslowrecord  > 0)  rslowrecord  -=1 
return ( 2) ; 

} 

{ 

r slowtime++ ; 
minute++ ; 
return(3) ; 

} 

{ 

minute++ ; 
return(4) ; 

} 

/*  white  space,  ignore  */ 

/*  white  space,  ignore  */ 

/*  and  count  the  number  of  columns  */ 


data2 
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irslowmin  #3 
inotrslowmin  #4 
> 

null  : 


*#*  rslow. semaction. c *** 


//include  <stdio.h> 
int  nsec=0; 


extern  int  minute; 
extern  int  rslowtime; 


describe ( ) 

{ 

if  ( rslowtime  >=  16  ) printf ( "H" ) ; 
else  if  ( rslowtime  > 8 ) printf ( "M" ) 
else  if  ( rslowtime  > 1) 

printf ("L") ; 
else  printf ( "N") ; 

if  ( minute  % 60  ==  0)  printfC'O); 

} 

semaction(n) 
int  n ; 


switch(n)  { 
case  0: 

break ; 

case  1 : 

if  ( nsec  < 3 ) { 

rslowtime  +=  nsec; 
nsec  = 0; 


} 

break ; 

case  2: 

if  ( nsec  < 
break ; 

case  3: 

if  ( nsec  < 
describe ( ) ; 
rslowtime  = 
break ; 

case  4: 

describe ( ) ; 
rslowtime  = 
break ; 

case  5: 

printf ( n0) ; 
break ; 


3 ) nsec++; 

3 ) rslowtime  += 
0; 

0; 


nsec  ; 
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default: 

printf( »#-*-#-*  semantic  action"); 
printf("  $d  not  handled. 0 ,n) ; 
break ; 


} 


} 
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